Class JF12FieldSolenoidal

Defined in File JF12FieldSolenoidal.h

Inheritance Relationships

Base Type

public JF12Field (Class JF12Field)

Class Documentation

class JF12FieldSolenoidal : public JF12Field 

JF12FieldSolenoidal galactic magnetic field model.

Implements a modified JF2012 Galactic magnetic field model. This implementation inherits most methods from the initial JF12Field, just the regular disk field and the poloidal halo “X field” are altered.

A solenoidal transition for the disk field (which redirects the magnetic flux of the magnetic spiral field at the boundaries and lets the field strength tend to 0 in a continuous way) and an improved version of the JF12 X field (with parabolic field lines at z = 0 instead of sharp kinks) were added in order to smooth the field and avoid violations of magnetic flux conservation, see: Kleimann et al, 2018, arXiv:1809.07528, Solenoidal improvements for the JF12 Galactic magnetic field model.

The turbulent field component is the exact same as in the initial JF12 implementation and should be used with care since the new regular and old turbulent field do not match in the transition regions of the spiral field.

Public Functions

JF12FieldSolenoidal(double delta = 3 * kpc, double zs = 0.5 * kpc)

Constructor

Parameters:

delta – Transition width for the disk field such that the magnetic flux of the spiral field lines is redirected between r = 5 kpc and r = 5 kpc + delta as well as r = 20 kpc - delta and r = 20 kpc. The input parameter delta should be non-negative and smaller than 7.5 kpc. The default value is 3 kpc.
zs – Scale height for the X-shaped poloidal halo field such that the straight field lines are replaced by parabolas for abs(z) < zs. The input parameter zs should be non-negative. The default value is 500 pc.

virtual void setUseStriatedField(bool use)

virtual void setUseTurbulentField(bool use)

void setDiskTransitionWidth(double delta)

Adjust the transition width of the disk field.

Parameters:: delta – The new transition width for the disk field with field strength transitions and flux redirection between r = 5 kpc and r = 5 kpc + delta as well as r = 20 kpc - delta and r = 20 kpc. Should be non-negative and smaller than 7.5 kpc.
Returns:: Void

double getDiskTransitionWidth() const

void setXScaleHeight(double zs)

Adjust the scale heigth of the X field.

Parameters:: zs – The new scale height of the X field with parabolic field lines for abs(z) < zs. Should be non-negative.
Returns:: Void

double getXScaleHeight() const

virtual Vector3d getXField(const double &r, const double &z, const double &sinPhi, const double &cosPhi) const

virtual Vector3d getDiskField(const double &r, const double &z, const double &phi, const double &sinPhi, const double &cosPhi) const

void deactivateOuterTransition()

Disable the transition of the spiral field strength to 0 at the outer boundary such that only the magnetic flux at the 5 kpc ring is redirected. Thus, the spiral field lines are continued to r = 20 kpc as in the initial JF12 field. You can reactivate the outer transition afterwards via setDiskTransitionWidth which sets both transition widths at the inner and outer boundary.

Returns:: Void

double getDiskTransitionPolynomial(const double &r) const

Evaluate the polynomial which provides the transition of the spiral field strength to zero in the transition regions. This transition is differentiable at the boundaries of the unaltered spiral field and contuinuous at the outer boundaries of the spiral field where the field strength goes to zero.

Parameters:: r – Distance of the current position to the z axis in the usual galactocentric cylindrical coordinates. Should be non-negative.
Returns:: The value of the transition polynomial at r if r is inside one of the transition regions. Otherwise, return (5 kpc)/r, i.e. the scaling of the spiral field, if r is inside the region where the spiral field remains unaltered.

double getDiskTransitionPolynomialDerivative(const double &r) const

Evaluate the derivative of the polynomial which provides the transition of the spiral field strength to zero in the transition regions. The derivative is needed to restore solenoidality in the transition region.

Parameters:: r – Distance of the current position to the z axis in the usual galactocentric cylindrical coordinates. Should be non-negative.
Returns:: The value of the transition polynomial derivative at r if r is inside one of the transition regions. Otherwise, return 0.

double getHPhiIntegral(const double &r, const double &phi) const

Evaluate an angular azimuth integral over the piecewise constant spiral field strengths at r_1 = 5 kpc. This integral is needed to restore the solenoidality of the spiral field in the transition regions.

Parameters:

r – Distance of the current position to the z axis in the usual galactocentric cylindrical coordinates. Should be non-negative.
phi – Azimuth angle of the current position in galactocentric cylindrical coordinates. Can be any double.

Returns:

The value of the azimuth integral over the spiral field strengths at r_1 = 5 kpc from at fixed phi0 to the current phi which is mapped back to r_1 = 5 kpc along the spiral field line passing through (r,phi).

double getSpiralFieldStrengthConstant(const double &r, const double &phi) const

Find the correct magnetic spiral arm for the current position and return its field strength at r_1 = 5 kpc.

Parameters:

r – Distance of the current position to the z axis in the usual galactocentric cylindrical coordinates. Should be non-negative.
phi – Azimuth angle of the current position in galactocentric cylindrical coordinates. Can be any double.

Returns:

The value of the spiral field strength b_j

void randomStriated(int seed = 0)

void randomTurbulent(int seed = 0)

void setStriatedGrid(ref_ptr<Grid1f> grid)

Set a striated grid and activate the striated field component

Parameters:: grid – scalar grid containing random +1/-1 values, 100 parsec grid spacing

void setTurbulentGrid(ref_ptr<Grid3f> grid)

Set a turbulent grid and activate the turbulent field component

Parameters:: grid – vector grid containing a random field of Brms = 1

ref_ptr<Grid1f> getStriatedGrid()

ref_ptr<Grid3f> getTurbulentGrid()

void setUseRegularField(bool use)

void setUseDiskField(bool use)

void setUseToroidalHaloField(bool use)

void setUseXField(bool use)

bool isUsingRegularField()

bool isUsingStriatedField()

bool isUsingTurbulentField()

bool isUsingDiskField()

bool isUsingToroidalHaloField()

bool isUsingXField()

double logisticFunction(const double &x, const double &x0, const double &w) const

Vector3d getRegularField(const Vector3d &pos) const

Vector3d getToroidalHaloField(const double &r, const double &z, const double &sinPhi, const double &cosPhi) const

Vector3d getStriatedField(const Vector3d &pos) const

double getTurbulentStrength(const Vector3d &pos) const

Vector3d getTurbulentField(const Vector3d &pos) const

virtual Vector3d getField(const Vector3d &pos) const

inline virtual Vector3d getField(const Vector3d &position, double z) const

inline size_t addReference() const

inline size_t removeReference() const

inline int removeReferenceNoDelete() const

inline size_t getReferenceCount() const

Protected Attributes

bool useRegularField

bool useStriatedField

bool useTurbulentField

bool useDiskField

bool useToroidalHaloField

bool useXField

double rArms[8]

double pitch

double sinPitch

double cosPitch

double tanPitch

double cotPitch

double tan90MinusPitch

double bDisk[11]

double bRing

double hDisk

double wDisk

double bNorth

double bSouth

double rNorth

double rSouth

double wHalo

double z0

double bX

double thetaX0

double sinThetaX0

double cosThetaX0

double tanThetaX0

double cotThetaX0

double rXc

double rX

double sqrtbeta

ref_ptr<Grid1f> striatedGrid

ref_ptr<Grid3f> turbulentGrid

double bDiskTurb[8]

double bDiskTurb5

double zDiskTurb

double bHaloTurb

double rHaloTurb

double zHaloTurb

mutable size_t _referenceCount