#include <deformednucleus.h>

DeformedNucleus: Child of nucleus for deformed nuclei.

All options from the nucleus will still apply. The deformed nucleus adds new or updated features which are outlined below.

Definition at line 36 of file deformednucleus.h.

Inheritance diagram for smash::DeformedNucleus:

[legend]

Collaboration diagram for smash::DeformedNucleus:

[legend]

Public Member Functions
	DeformedNucleus (const std::map< PdgCode, int > &particle_list, int nTest)
	Constructor for DeformedNucles which takes a particle list and the number of testparticles. More...

	DeformedNucleus (Configuration &config, int nTest, bool auto_deformation)
	Constructor for DeformedNucleus, that needs the configuration parameters from the inputfile and the number of testparticles. More...

ThreeVector	distribute_nucleon () override
	Deformed Woods-Saxon sampling routine. More...

void	set_deformation_parameters_automatic ()
	Sets the deformation parameters of the radius according to the current mass number. More...

void	set_deformation_parameters_from_config (Configuration &config)
	Set parameters for spherical deformation of the nucleus from the values specified in the configuration file. More...

void	set_orientation_from_config (Configuration &orientation_config)
	Set angles for orientation of nucleus from config file. More...

void	rotate () override
	Rotates the nucleus according to members nucleus_polar_angle_ and nucleus_azimuthal_angle_ and updates nucleon positions. More...

double	get_saturation_density () const

double	nucleon_density (double r, double cosx) const override
	Return the deformed Woods-Saxon probability density for the given position. More...

double	nucleon_density_unnormalized (double r, double cosx) const override
	Return the unnormalized deformed Woods-Saxon distribution for the given position. More...

void	set_beta_2 (double b2)
	Set deformation coefficient for Y_2_0. More...

void	set_beta_4 (double b4)
	Set deformation coefficient for Y_4_0. More...

void	set_polar_angle (double theta)
	Set the nucleus polar angle. More...

void	set_azimuthal_angle (double phi)
	Set the nucleus azimuthal angle. More...

double	get_beta2 ()
	return the beta2 value. More...

double	get_beta4 ()
	return the beta4 value. More...

Public Member Functions inherited from smash::Nucleus
	Nucleus ()=default
	default constructor More...

	Nucleus (Configuration &config, int nTest)
	Constructor for Nucleus, that needs the configuration parameters from the inputfile and the number of testparticles. More...

	Nucleus (const std::map< PdgCode, int > &particle_list, int nTest)
	Constructor which directly initializes the Nucleus with particles and respective counts. More...

virtual	~Nucleus ()=default

double	mass () const

double	woods_saxon (double x)
	Woods-Saxon distribution. More...

virtual void	arrange_nucleons ()
	Sets the positions of the nucleons inside a nucleus. More...

virtual void	set_parameters_automatic ()
	Sets the deformation parameters of the Woods-Saxon distribution according to the current mass number. More...

virtual void	set_parameters_from_config (Configuration &config)
	Sets the parameters of the Woods-Saxon according to manually added values in the configuration file. More...

virtual void	generate_fermi_momenta ()
	Generates momenta according to Fermi motion for the nucleons. More...

void	boost (double beta_scalar)
	Boosts the nuclei into the computational frame, such that the nucleons have the appropriate momentum and the nuclei are lorentz-contracted. More...

void	fill_from_list (const std::map< PdgCode, int > &particle_list, int testparticles)
	Adds particles from a map PDG code => Number_of_particles_with_that_PDG_code to the nucleus. More...

void	shift (double z_offset, double x_offset, double simulation_time)
	Shifts the nucleus to correct impact parameter and z displacement. More...

void	copy_particles (Particles *particles)
	Copies the particles from this nucleus into the particle list. More...

size_t	size () const
	Number of numerical (=test-)particles in the nucleus: More...

size_t	number_of_particles () const
	Number of physical particles in the nucleus: More...

size_t	number_of_protons () const
	Number of physical protons in the nucleus: More...

FourVector	center () const
	Calculate geometrical center of the nucleus. More...

void	set_label (BelongsTo label)
	Sets target / projectile labels on nucleons. More...

void	align_center ()
	Shifts the nucleus so that its center is at (0,0,0) More...

virtual double	calculate_saturation_density () const

virtual void	set_saturation_density (double density)
	Sets the saturation density of the nucleus. More...

std::vector< ParticleData >::iterator	begin ()
	For iterators over the particle list: More...

std::vector< ParticleData >::iterator	end ()
	For iterators over the particle list: More...

std::vector< ParticleData >::const_iterator	cbegin () const
	For const iterators over the particle list: More...

std::vector< ParticleData >::const_iterator	cend () const
	For const iterators over the particle list: More...

void	set_diffusiveness (double diffuse)
	Sets the diffusiveness of the nucleus. More...

double	get_diffusiveness () const

double	get_saturation_density () const

double	default_nuclear_radius ()
	Default nuclear radius calculated as: More...

void	set_nuclear_radius (double rad)
	Sets the nuclear radius. More...

double	get_nuclear_radius () const

Private Attributes
double	beta2_ = 0.0
	Deformation parameter for angular momentum l=2. More...

double	beta4_ = 0.0
	Deformation parameter for angular momentum l=4. More...

Angles	nuclear_orientation_
	Nucleus orientation (initial profile in xz plane) in terms of a pair of angles (theta, phi) More...

bool	random_rotation_ = false
	Whether the nuclei should be rotated randomly. More...

Additional Inherited Members
Protected Member Functions inherited from smash::Nucleus
void	random_euler_angles ()
	Randomly generate Euler angles. More...

Protected Attributes inherited from smash::Nucleus
std::vector< ParticleData >	particles_
	Particles associated with this nucleus. More...

double	saturation_density_ = nuclear_density
	Saturation density of this nucleus. More...

double	euler_phi_
	Euler angel phi. More...

double	euler_theta_
	Euler angel theta. More...

double	euler_psi_
	Euler angel psi. More...

Constructor & Destructor Documentation

◆ DeformedNucleus() [1/2]

smash::DeformedNucleus::DeformedNucleus	(	const std::map< PdgCode, int > &	particle_list,
		int	nTest
	)

Constructor for DeformedNucles which takes a particle list and the number of testparticles.

This constructor is only used for testing purposes.

Parameters

[in]	particle_list	Map with PDGCode and number of particles which make up the nucleus
[in]	nTest	number of testparticles

Definition at line 93 of file deformednucleus.cc.

95 : Nucleus(particle_list, nTest) {}

smash::Nucleus::Nucleus

Nucleus()=default

default constructor

◆ DeformedNucleus() [2/2]

smash::DeformedNucleus::DeformedNucleus	(	Configuration &	config,
		int	nTest,
		bool	auto_deformation
	)

Constructor for DeformedNucleus, that needs the configuration parameters from the inputfile and the number of testparticles.

Parameters

[in]	config	contains the parameters from the inputfile on the numbers of particles with a certain PDG code
[in]	nTest	number of testparticles
[in]	auto_deformation	whether or not deformation parameters should be set automatically

Definition at line 97 of file deformednucleus.cc.

     : Nucleus(config, nTest) {
   if (auto_deformation) {
     set_deformation_parameters_automatic();
     set_saturation_density(calculate_saturation_density());
   } else {
     set_deformation_parameters_from_config(config);
   }
   if (config.has_value({"Deformed", "Orientation"})) {
     Configuration subconfig = config["Deformed"]["Orientation"];
     set_orientation_from_config(subconfig);
   }
 }

Here is the call graph for this function:

Member Function Documentation

◆ distribute_nucleon()

ThreeVector smash::DeformedNucleus::distribute_nucleon ( )

overridevirtual

Deformed Woods-Saxon sampling routine.

Returns: Spatial position from uniformly sampling the deformed woods-saxon distribution

Reimplemented from smash::Nucleus.

Definition at line 112 of file deformednucleus.cc.

                                                 {
   double a_radius;
   Angles a_direction;
   // Set a sensible maximum bound for radial sampling.
   double radius_max =
       Nucleus::get_nuclear_radius() / Nucleus::get_diffusiveness() +
       Nucleus::get_nuclear_radius() * Nucleus::get_diffusiveness();
  
   // Sample the distribution.
   do {
     a_direction.distribute_isotropically();
     // sample r**2 dr
     a_radius = radius_max * std::cbrt(random::canonical());
   } while (random::canonical() >
            nucleon_density(a_radius, a_direction.costheta()) /
                Nucleus::get_saturation_density());
  
   // Update (x, y, z) positions.
   return a_direction.threevec() * a_radius;
 }

Here is the call graph for this function:

◆ set_deformation_parameters_automatic()

void smash::DeformedNucleus::set_deformation_parameters_automatic ( )

Sets the deformation parameters of the radius according to the current mass number.

The deformation parameters are taken from Moller:1993ed [35]. Corrections to the deformation parameter beta2 in Uranium come from Kuhlman:2005ts [29]. For finite nucleon size corrections to the nuclear density and radius for copper and gold, see Hirano:2009ah [22], and Hirano:2010jg [24] for uranium.

Definition at line 133 of file deformednucleus.cc.

                                                            {
   // Set the deformation parameters
   // reference for U, Pb, Au, Cu: \iref{Moller:1993ed}
   // reference for Zr and Ru: \iref{Schenke:2019ruo}
   // reference for Xe: \iref{Moller:2015fba}
   bool listed = 0;
   const std::map<int, std::string> A_map = {{238, "Uranium"},
                                             {208, "Lead"},
                                             {197, "Gold"},
                                             {63, "Copper"},
                                             {129, "Xenon"}};
   const std::map<std::string, std::string> Z_map = {{"Uranium", "92"},
                                                     {"Lead", "82"},
                                                     {"Gold", "79"},
                                                     {"Copper", "29"},
                                                     {"Xenon", "54"}};
   int A = Nucleus::number_of_particles();
   int Z = Nucleus::number_of_protons();
   switch (A) {
     case 238:  // Uranium
       if (Z == 92) {
         set_beta_2(0.28);
         set_beta_4(0.093);
       } else {
         listed = true;
       }
       break;
     case 208:  // Lead
       if (Z == 82) {
         set_beta_2(0.0);
         set_beta_4(0.0);
       } else {
         listed = true;
       }
       break;
     case 197:  // Gold
       if (Z == 79) {
         set_beta_2(-0.131);
         set_beta_4(-0.031);
       } else {
         listed = true;
       }
       break;
     case 129:  // Xenon
       if (Z == 54) {
         set_beta_2(0.162);
         set_beta_4(-0.003);
       } else {
         listed = true;
       }
       break;
     case 63:  // Copper
       if (Z == 29) {
         set_beta_2(0.162);
         set_beta_4(-0.006);
       } else {
         listed = true;
       }
       break;
     case 96:
       if (Z == 40) {  // Zirconium
         set_beta_2(0.0);
         set_beta_4(0.0);
       } else if (Z == 44) {  // Ruthenium
         set_beta_2(0.158);
         set_beta_4(0.0);
       } else {
         throw std::domain_error(
             "Number of protons for nuclei with mass number A = 96 does not "
             "match that of Zirconium or Ruthenium. The deformation parameters "
             "for additional isobars are currently not implemented."
             " Please specify at least \"Beta_2\" and \"Beta_4\" "
             "manually and set \"Automatic: False.\" ");
       }
       break;
     default:
       throw std::domain_error(
           "Mass number not listed for automatically setting deformation "
           "parameters. Please specify at least \"Beta_2\" and \"Beta_4\" "
           "manually and set \"Automatic: False.\" ");
   }
   if (listed) {
     throw std::domain_error("Mass number is listed under " + A_map.at(A) +
                             " but the proton "
                             "number of " +
                             std::to_string(Z) +
                             " does not match "
                             "its " +
                             Z_map.at(A_map.at(A)) +
                             " protons."
                             "Please specify at least \"Beta_2\" and \"Beta_4\" "
                             "manually and set \"Automatic: False.\" ");
   }
 }

Here is the call graph for this function:

Here is the caller graph for this function:

◆ set_deformation_parameters_from_config()

void smash::DeformedNucleus::set_deformation_parameters_from_config ( Configuration & config )

Set parameters for spherical deformation of the nucleus from the values specified in the configuration file.

Parameters

config The configuration for the deformation of this nucleus (projectile or target).

Definition at line 228 of file deformednucleus.cc.

                            {
   // Deformation parameters.
   if (config.has_value({"Deformed", "Beta_2"})) {
     set_beta_2(static_cast<double>(config.take({"Deformed", "Beta_2"})));
   }
   if (config.has_value({"Deformed", "Beta_4"})) {
     set_beta_4(static_cast<double>(config.take({"Deformed", "Beta_4"})));
   }
 }

Here is the call graph for this function:

Here is the caller graph for this function:

◆ set_orientation_from_config()

void smash::DeformedNucleus::set_orientation_from_config ( Configuration & orientation_config )

Set angles for orientation of nucleus from config file.

Parameters

[in] orientation_config The configuration for the deformation of this nucleus (projectile or target).

Definition at line 239 of file deformednucleus.cc.

                                        {
   // Read in orientation if provided, otherwise, the defaults are
   // theta = pi/2, phi = 0, as declared in the angles class
  
   if (orientation_config.has_value({"Theta"})) {
     if (orientation_config.has_value({"Random_Rotation"}) &&
         orientation_config.take({"Random_Rotation"})) {
       throw std::domain_error(
           "Random rotation of nuclei is activated although"
           " theta is provided. Please specify only either of them. ");
     } else {
       set_polar_angle(static_cast<double>(orientation_config.take({"Theta"})));
     }
   }
  
   if (orientation_config.has_value({"Phi"})) {
     if (orientation_config.has_value({"Random_Rotation"}) &&
         orientation_config.take({"Random_Rotation"})) {
       throw std::domain_error(
           "Random rotation of nuclei is activated although"
           " phi is provided. Please specify only either of them. ");
     } else {
       set_azimuthal_angle(
           static_cast<double>(orientation_config.take({"Phi"})));
     }
   }
  
   if (orientation_config.take({"Random_Rotation"}, false)) {
     random_rotation_ = true;
   }
 }

Here is the call graph for this function:

Here is the caller graph for this function:

◆ rotate()

void smash::DeformedNucleus::rotate ( )

overridevirtual

Rotates the nucleus according to members nucleus_polar_angle_ and nucleus_azimuthal_angle_ and updates nucleon positions.

Reimplemented from smash::Nucleus.

Definition at line 272 of file deformednucleus.cc.

                              {
   if (random_rotation_) {
     // Randomly generate euler angles for theta and phi. Psi needs not be
     // assigned, as the nucleus objects are symmetric with respect to psi.
     Nucleus::random_euler_angles();
     set_azimuthal_angle(euler_phi_);
     set_polar_angle(euler_theta_);
   }
   for (auto &particle : *this) {
     /* Rotate every vector by the nuclear azimuth phi and polar angle
      * theta (the Euler angles). This means applying the matrix for a
      * rotation of phi about z, followed by the matrix for a rotation
      * theta about the rotated x axis. The third angle psi is 0 by symmetry.*/
     ThreeVector three_pos = particle.position().threevec();
     three_pos.rotate(nuclear_orientation_.phi(), nuclear_orientation_.theta(),
                      0.);
     particle.set_3position(three_pos);
   }
 }

Here is the call graph for this function:

◆ get_saturation_density()

double smash::DeformedNucleus::get_saturation_density ( ) const

inline

Returns: the saturation density of the deformed_nucleus

See also: saturation_density_

Definition at line 100 of file deformednucleus.h.

100 { return saturation_density_; }

smash::Nucleus::saturation_density_

double saturation_density_

Saturation density of this nucleus.

Definition: nucleus.h:292

◆ nucleon_density()

double smash::DeformedNucleus::nucleon_density	(	double	r,
		double	cosx
	)		const

overridevirtual

Return the deformed Woods-Saxon probability density for the given position.

This corresponds to the nuclear density at the very same position.

Parameters

[in]	r	The radius at which to sample
[in]	cosx	The cosine of the polar angle at which to sample

Returns: The Woods-Saxon density

Reimplemented from smash::Nucleus.

Definition at line 304 of file deformednucleus.cc.

                                                                    {
   return Nucleus::get_saturation_density() /
          (1 + std::exp((r - Nucleus::get_nuclear_radius() *
                                 (1 + beta2_ * y_l_0(2, cosx) +
                                  beta4_ * y_l_0(4, cosx))) /
                        Nucleus::get_diffusiveness()));
 }

Here is the call graph for this function:

Here is the caller graph for this function:

◆ nucleon_density_unnormalized()

double smash::DeformedNucleus::nucleon_density_unnormalized	(	double	r,
		double	cosx
	)		const

overridevirtual

Return the unnormalized deformed Woods-Saxon distribution for the given position.

Parameters

[in]	r	The radius
[in]	cosx	The cosine of the polar angle

Returns: The unnormalized Woods-Saxon distribution

Reimplemented from smash::Nucleus.

Definition at line 312 of file deformednucleus.cc.

                                                                         {
   return 1.0 / (1 + std::exp((r - Nucleus::get_nuclear_radius() *
                                       (1 + beta2_ * y_l_0(2, cosx) +
                                        beta4_ * y_l_0(4, cosx))) /
                              Nucleus::get_diffusiveness()));
 }

Here is the call graph for this function:

◆ set_beta_2()

void smash::DeformedNucleus::set_beta_2 ( double b2 )

inline

Set deformation coefficient for Y_2_0.

Parameters

[in] b2 deformation coefficient for l=2

Definition at line 123 of file deformednucleus.h.

123 { beta2_ = b2; }

Here is the caller graph for this function:

◆ set_beta_4()

void smash::DeformedNucleus::set_beta_4 ( double b4 )

inline

Set deformation coefficient for Y_4_0.

Parameters

[in] b4 deformation coefficient for l=4

Definition at line 128 of file deformednucleus.h.

128 { beta4_ = b4; }

Here is the caller graph for this function:

◆ set_polar_angle()

void smash::DeformedNucleus::set_polar_angle ( double theta )

inline

Set the nucleus polar angle.

Parameters

[in] theta Polar angle of position inside nucleus

Definition at line 133 of file deformednucleus.h.

                                             {
     nuclear_orientation_.set_theta(theta);
   }

Here is the call graph for this function:

Here is the caller graph for this function:

◆ set_azimuthal_angle()

void smash::DeformedNucleus::set_azimuthal_angle ( double phi )

inline

Set the nucleus azimuthal angle.

Parameters

[in] phi Azimuthal angle of position inside nucleus

Definition at line 140 of file deformednucleus.h.

                                               {
     nuclear_orientation_.set_phi(phi);
   }

Here is the call graph for this function:

Here is the caller graph for this function:

◆ get_beta2()

double smash::DeformedNucleus::get_beta2 ( )

inline

return the beta2 value.

Definition at line 146 of file deformednucleus.h.

146 { return beta2_; }

◆ get_beta4()

double smash::DeformedNucleus::get_beta4 ( )

inline

return the beta4 value.

Definition at line 150 of file deformednucleus.h.

150 { return beta4_; }

Member Data Documentation

◆ beta2_

double smash::DeformedNucleus::beta2_ = 0.0

private

Deformation parameter for angular momentum l=2.

Definition at line 154 of file deformednucleus.h.

◆ beta4_

double smash::DeformedNucleus::beta4_ = 0.0

private

Deformation parameter for angular momentum l=4.

Definition at line 156 of file deformednucleus.h.

◆ nuclear_orientation_

Angles smash::DeformedNucleus::nuclear_orientation_

private

Nucleus orientation (initial profile in xz plane) in terms of a pair of angles (theta, phi)

Definition at line 161 of file deformednucleus.h.

◆ random_rotation_

bool smash::DeformedNucleus::random_rotation_ = false

private

Whether the nuclei should be rotated randomly.

Definition at line 165 of file deformednucleus.h.

The documentation for this class was generated from the following files:

src/include/smash/deformednucleus.h
src/deformednucleus.cc

Public Member Functions

Private Attributes

Additional Inherited Members

Constructor & Destructor Documentation

◆ DeformedNucleus() [1/2]

◆ DeformedNucleus() [2/2]

Member Function Documentation

◆ distribute_nucleon()

◆ set_deformation_parameters_automatic()

◆ set_deformation_parameters_from_config()

◆ set_orientation_from_config()

◆ rotate()

◆ get_saturation_density()

◆ nucleon_density()

◆ nucleon_density_unnormalized()

◆ set_beta_2()

◆ set_beta_4()

◆ set_polar_angle()

◆ set_azimuthal_angle()

◆ get_beta2()

◆ get_beta4()

Member Data Documentation

◆ beta2_

◆ beta4_

◆ nuclear_orientation_

◆ random_rotation_