Asymmetric Backward Peaking Radiation Pattern From a Relativistic Particle Accelerated by Lightning Leader Tip Electric Field

Terrestrial Gamma ray Flashes exhibit slopes of ionizing radiation associated with bremsstrahlung. Bremsstrahlung has a continuous spectrum of radiation from radio waves to ionizing radiation. The Poynting vector of the emitted radiation, that is, the radiation pattern around a single particle under the external lightning electric field during interaction with other particles or atoms, is not quite well known. The overall radiation pattern arises from the combination of radiation of parallel and perpendicular motions of a particle caused by the acceleration from the lightning electric field and the bremsstrahlung. The calculations and displays of radiation patterns are generally limited to a low‐frequency approximation for radio waves and separate parallel and perpendicular motions. Here, we report the radiation patterns of combined parallel and perpendicular motions from accelerated relativistic particles at low and high frequencies of the bremsstrahlung process with an external lightning electric field. The primary outcome is that radiation patterns have four relative maxima with two forward peaking and two backward peaking lobes. The asymmetry of the radiation pattern, that is, the different intensities of forward and backward peaking lobes, are caused by the Doppler effect. A novel outcome is that bremsstrahlung has an asymmetry of the four maxima around the velocity vector caused by the curvature of the particle's trajectory as it emits radiation. This mathematical modeling helps to better understand the physical processes of a single particle's radiation pattern, which might assist the interpretation of observations with networks of radio receivers and arrays of γ‐ray detectors.


Introduction
It was recently suggested that high-frequency radiation emissions observed in the atmosphere could originate from muons interacting with electric fields inside thunderclouds.This novel idea is based on a reduction of the muon detection during thunderstorm occurrences by the ground based telescope GRAPES-3 located in Ooty, India (Hariharan et al., 2019).Gamma-Ray Bursts (GRBs) are commonly thought to result from the interaction of neutron stars in outer space or comet collisions.GRBs emit photons in the energy the atmosphere (Celestin, 2016).Another working hypothesis is that bremsstrahlung radiation is emitted by thermal runaway electrons accelerated by intra-cloud lightning leader tips (Xu et al., 2015).Bremsstrahlung has a continuous electromagnetic spectrum.Low-frequency radio and optical emissions could also be due to fluorescence, where high-frequency TGFs are absorbed by air molecules (Xu et al., 2015).Numerical Monte Carlo simulations demonstrated the significance of the bremsstrahlung process as the primary process behind high-frequency emissions (Dwyer et al., 2012).Bremsstrahlung electrons emit radiation in forward peaking radiation patterns with an angle that scales with the inverse of the Lorentz factor of the relativistic electrons (Koch & Motz, 1959).
Asymmetric signal of γ-ray bursts measured by the GRB Monitor on the Fermi Gamma ray Space Telescope reveal the lightning leader charge structure.Asymmetric γ-ray pulses indicate the lightning leader charge flux, which exhibits a fast rise and slow decay of the leader tip electric field (Foley et al., 2014).The asymmetries in γ-ray pulses are thought to be caused by Compton scattering (Xu et al., 2019).The rise to decay time ratio of single γ-ray pulses was measured to be approximately 0.67 (Nemiroff et al., 1994).Data from the Burst and Transient Source Experiment reveal two different types of spectra of γ-ray bursts known as bright and dim GRBs.It was found that dim GRBs have less photon energy than bright GRBs (Norris et al., 1994).It was observed that as time passes, overall γ-ray photons transit from bright to dim photons as a photon bunch due to a time delay of approximately 100 μs between the peaks arising from hard and soft photons (Grefenstette et al., 2008).
Experimental measurements of ionizing radiation and optical emissions by the Atmosphere Space Interactions Monitor on the International Space Station recently reported the detection of 217 TGFs from June 2, 2018, to April 1, 2019 (Østgaard et al., 2019), some associated with radio emissions from charged particles that are observed on the ground.All these measurements reveal the properties of γ-ray bursts.After the combination of the measurements from ground-based radio receivers and spacecraft, it was found that TGFs are produced at the very beginning of the lightning discharge process.It is well known that the observed γ-rays originate from the bremsstrahlung process (Xu et al., 2015).There are approximately 10 17 -10 19 GRBs emitted during the bremsstrahlung process.It is well known that the initially emitted ionizing radiation is not the same in terms of energy and direction compared to the radiation measured by sensors.This difference is because the emitted radiation loses energy by back-scattering and interacting with other air molecules.The interaction causes an ionization and releases more electrons, which can explain why 10 17 -10 19 γ-rays are emitted (Dwyer, 2008).Another theory explains γ-ray bursts to originate from the large electric fields of leader tip streamers producing ∼10 12 electrons which then increase the number of electrons within the relativistic runaway electron avalanche (RREA) process that emits γ-ray photons (Babich et al., 2014(Babich et al., , 2015;;Celestin & Pasko, 2011;Chanrion & Neubert, 2010;Moss et al., 2006;Skeltved et al., 2017).This contribution reports the modeling of an asymmetric forward peaking radiation pattern and an asymmetric backward peaking radiation pattern of a single particle bremsstrahlung process.The asymmetry occurs around the horizontal axis parallel to the direction of motion of the charged particle, and it is unique to the bremsstrahlung process as the particle continuously follows a curved trajectory of an anticlockwise rotation.Radiation patterns are calculated for both relativistic and non-relativistic velocities.The main asymmetry with four radiation peaks is unique to the bremsstrahlung process and occurs when the particle radiation transits from a dipole toward forward and backward peaking radiation patterns.
The particle trajectory is mostly determined by the Coulomb field and modified by the presence of the external lightning leader tip electric field.The change of the particle trajectory caused by the external leader tip electric field is accounted for by only adjusting the radius of the curvature of the particle's trajectory over time.This is possible because, when the leader tip electric field is resolved into its vector components, one vector component is tangent to the trajectory, whereas, the other vector component is perpendicular to it.The perpendicular component points either in the same or opposite direction to the centripetal acceleration caused by the target particle depending on the polarity of the leader tip electric field.The tangent component of the leader tip electric field contributes to determining the amount of angular rotation that can be covered by the particle during Coulomb interaction.The perpendicular component of the external leader tip electric field contributes to the radius of curvature of the curved trajectory caused by Coulomb interaction.
The particle trajectory defined by the Coulomb field and the particle velocity arising from the leader tip electric field are two independently derived equations that serve as ingredients that contribute to the bremsstrahlung radiation pattern (Figure 1a).This independence is important because, under a repulsive or attractive Coulomb force, the target particle mostly determines the trajectory of the incoming particle.Also, parameters defining particle trajectory by the target particles Coulomb field can account for the trajectory changes caused by the external leader tip electric field.Hence, when the radius of curvature modification caused by the external leader tip electric field is excluded from the particle trajectory, the temporal derivative of the particle trajectory is the velocity that an incoming particle experiences from the Coulomb force.The overall particle velocity is the superposition of the velocity arising from both the leader tip and the target particle's electric field.The particle velocity is considered to be dominated by the external lightning leader tip electric field.

Aims and Objectives
The primary aim of this contribution is to construct a generic mathematical model of electromagnetic radiation emitted by a single relativistic charged particle via the bremsstrahlung process that applies to high and low-frequency radiation.Moreover, a generic mathematical model that is also capable of explaining a particle motion whose velocity and the acceleration vectors have both parallel (particle following linear path) and perpendicular (particle following circular trajectory) orientations with each other, at the same time.The secondary aim is to use the model to investigate the details of the actual particle motion and ra-YÜCEMÖZ AND FÜLLEKRUG 10.1029/2020JD033204 3 of 23 diation mechanisms during the bremsstrahlung process.Initial assumptions for the particle position vector and lightning leader tip electric field were made to achieve the stated aims.In addition, all the necessary physical details of the interaction, such as form factors and collision cross-sections, are included and investigated.They both play a crucial role in explaining high-frequency radiation.Moreover, the Doppler effect was added to transform the model from a static to a dynamic model.Transformation enables an explanation of the asymmetry of the radiation in forward and backward direction relative to the axis perpendicular to the particle's direction of motion.Finally, the second-order differential equation (Equation 13) was solved as an initial step to establish a generic model to explain radiation patterns using algebra, calculus, geometry, and the table of integral transforms.

Particle Position Vector
Starting with defining a curved path for a bremsstrahlung electron.
The position vector is formulated for a particle trajectory that is an anti-clockwize rotating spiral as a function of the retarded time characteristic for bremsstrahlung radiation.
The position vector in Equation 1 defines a spiral trajectory for an incoming particle, that is, an electron, induced by the Coulomb force of the target particle that causes the emission of bremsstrahlung radiation (Figure 2a).The spiral trajectory in Figure 2a and mathematically defined in Equation 1 is realistic even though the mean free path is quite short, e.g., nm-μm in the atmosphere with a high recombination rate.For example, a circle with a radius of 1 m could also have a radius of 2 μm, depending on the medium and YÜCEMÖZ AND FÜLLEKRUG   1).The tangential component of the electron acceleration and corresponding velocity vectors are displayed with black arrows and are tangential to the spiral trajectory (red line) and perpendicular to the position vector r(t).R(t r ) is the distance between the accelerated electron and the observer, which is a function of both retarded and chronological time (t r , t).P is the position of an observer.
 S is the Poynting vector, or radiant energy flux, which determines the direction of the energy flow per area of an emitted electromagnetic wave.The dimensionless unit vector n points in the direction of the Poynting vector  S. (b) The position vector r(t), velocity vector v(t) and the unit radiation vector n form a rectangular triangle.(c) Co-ordinate system and all the associated vectors are drawn from the target particle which is positioned at the origin, O.As the target particle is not a rigid body and displaces under the incoming particle's Coulomb field, the moving co-ordinate system is constructed (Brock, 2019).Displacement of the moving co-ordinate system is described with the collision cross-section and form factors (Supplemental information, Text S2).Position vector, r(t) connects both target and incoming particle together.Time derivative of this position vector is the tangential velocity of the incoming particle and it is perpendicular to the position vector in the moving co-ordinate system.Without radius of trajectory changes caused by the external leader tip electric field, time derivative of the position vector is the tangential component of the Coulomb velocity.With radius of trajectory changes caused by the external leader tip electric field, time derivative of the position vector is the tangential component of the overall velocity.
the recombination rate.It is still a circle but a scaled microscopic version of the initial macroscopic circle.Preserving geometry at different scales is also true in the formulated spiral trajectory (Equation 1).The decision on a specific particle trajectory considers the ratio of particle size to a curved trajectory radius.If the particle size is larger than the curvature radius, the particle trajectory is approximately a straight line.Therefore, a spiral particle trajectory is realistic because an electron has a size of <2.8 × 10 −19 m as measured by the Hadron-Electron Ring Accelerator in Hamburg, Germany at the Deutsches Electronen Synchrotron facility (Bourilkov, 2000).
Overall, acceleration caused by the Coulomb field can be resolved into its orthogonal components as centripetal and tangent Coulomb accelerations.The centripetal acceleration caused by the target particle defines the curved trajectory of the incoming particle and it is this centripetal acceleration that forces the incoming particle to stay on its curved trajectory.The tangent Coulomb velocity and acceleration components caused by the target particle due to the Coulomb field are perpendicular to the position and centripetal acceleration vectors.Hence, tangent to the particle's curved trajectory.
Similar to Coulomb field, leader tip electric field can also be resolved into its orthogonal vector components.As shown in Figure 1b, tangent to the particle trajectory component of the leader tip electric field increases the incoming particle's tangential Coulomb velocity, hence acceleration.The remaining acceleration component of the external lightning electric field acts in the opposite direction to centripetal acceleration and it only changes the radius of curvature of the incoming particle's trajectory which can be accounted with mathematical variable "b" in Equation 1 representing particles radius of curvature.
where r(t) is the position vector as a function of time t in s, R is the dimensionless bremsstrahlung asymmetry index, τ is the mean free time in s.Also, ω′ is the angular frequency of the emitted electromagnetic wave in the frame of reference of the particle in rads/s, c is the speed of light, θ n,r(t) is the angle between the emitted radiation unit vector n and the particle's position vector.The time range of the position vector is −∞ < t < +∞.In addition, the factor b in m describes the interaction distance between the incoming particle and the target particle, which is the radius of the time-dependent position vector.The radius of the position vector is directly proportional to the parameter b.This radius of curvature is related to the mean free path because the curvature increases with time, contributing to the overall arc length, i.e., the mean free path of the accelerated particle.In other words, b ∝ λ v/c , where λ v/c is the mean free path of a particle at the velocity v which is given as a percentage of the speed of light v/c.The relation between a and λ v/c will be derived in Section 3. The factor a in m is an arbitrary adjustment parameter.It is introduced to correct the radius of the curvature of a particle during the bremsstrahlung process.This correction is required because the trajectory of a relativistic particle shrinks in size over time, and a propagation close to the speed of light introduces significant changes in the mean free path.Finally, Equation 1 is a function of time t.Hence, the particle will only cover some segment, or arc length, of the spiral, or the complete arc length of a particle's spiral trajectory when t = τ.

Interaction Gap
Defining the mean free path, λ v/c and the particle velocity arising only from the Coulomb interaction, dr dt .
The arc length, or the mean free path, of a curve in polar coordinates is given by where λ v/c is the mean free path, r is the position vector given in Equation 1, and the derivative with respect to the polar angle is  dr d .Equation 2 is a geometric equation used to calculate the arc length of any curve in polar coordinates.As the particle covers its mean free path, it follows the curved spiral trajectory in small segments with time t.In this case, Equation 2 has to be formulated with the parameter t rather than θ.To satisfy the dimensional accuracy, the Jacobian is used to carry out a variable transform from θ to t to preserve the dimension m 2 in the square root of the integrand.In geometry, the unit of angle in rad is considered to be dimensionless.An example of this is the arc length formula of a circle (s = rθ).Hence,    t .
The mean free time τ is a function of time due to its dependence on particle velocity that changes over time.
Therefore, the quotient rule has to be used to find the derivative where τ′ is the first derivative of the mean free time with respect to time t.The mean free time τ of the two identical particles is defined as , where n is the particle or molecule number density per unit volume in m −3 and d is the diameter in m of both an incoming and a target particle or molecule.The basic mean free time is defined between two identical particles with same diameter d where effective collision area is defined as a function of the two identical particle's diameter as A c = πd 2 .Finally, and the limits of the integration are      0 1 f t rad.The maximum limit of the integral is 1, which means that the ratio of t to τ as time progresses should not exceed allowed mean free time τ for the particle to radiate.
Coulomb velocity of an incoming particle arising from electrostatic interaction between a target particle can be written as Second terms, t 2R 2R′ln(t) and τ 2R 2R′ln(τ) inside the square brackets of Equation 3 represents a trajectory of varying radius of curvature, that is, spiral.Changing radius of curvature with time is directly linked to the changing bremsstrahlung asymmetry, R with time, R′.In the case of R′ = 0, particle follows a constant radius trajectory where the bremsstrahlung asymmetry, R remains at a constant value.Such trajectory could be a circle.
This leads to Equation 4 means that by knowing the mean free path of an electron in the atmosphere, the parameters a and b of a position vector r(t) in equation one can be calculated.

Lightning Leader Tip Electric Field
In this section, the lightning leader tip electric field is formulated to define the particle velocity for the bremsstrahlung radiation.The velocity defined by an integration of the electric field can be equated to the formulated velocity (Equation 3) from the particle position vector to calculate the unknown parameters a and b.Moreover, this is an important step in defining particle velocity parameters such as the mean free path and time as the leader tip electric fields can be measured experimentally.
YÜCEMÖZ AND FÜLLEKRUG 10.1029/2020JD033204 6 of 23 The measured lightning discharge electric field has a characteristic behavior of a sharp rise and a slow decay.This characteristic information was used to mathematically model the lightning leader tip electric field, E in Equation 5.The electric field is mathematically modeled using the ansatz of a sharp rise and a slow decay where ], n = 1.95 is a dimensionless constant, and m = 1.5 is a dimensionless constant.The scaling factor s f = 1 scales the lightning leader tip electric field magnitude, i.e., it is a dimensionless constant that can be used with different constant values to adjust the leader tip electric field to the required strength, s ft = 1 is a dimensionless constant scaling factor of time, and the time t/t 0 is the relative time with respect to t 0 = 1 s where t[s] is an independent time variable.The time range of the electric field is s ft 0.39694 μs ≤ t ≤ ∞.Values closer to 0.397 μs describe the initial state and values closer to 72.2 μs describe the final state of the lightning leader tip electric field shown in Figure 3.The minimum value of time was determined such that it is the minimum value of the electric field before it goes to −∞ at t = 0 s.The scaling factor determines the acceleration of the particle.As a result, the magnitude of the emitted radiation intensity scales in an indirect way by scaling the electric field which scales the force and hence the acceleration.
The lightning leader tip electric field in Equation 5is derived such that it goes to zero at infinite time, as shown in Figure 3b.Any remaining residual electric field integrated over infinite time would accelerate a particle to velocities exceeding the speed of light.
The strength and duration of the electric field affect the magnitude and pattern of the emitted radiation.
To adjust the electric field's strength and duration conveniently without affecting the graphical shape of Figure 3a, the scaling factor s f for magnitude and s ft for the duration are introduced.

Velocity Function
The velocity of the particle accelerated by the leader tip electric field is found by integrating the derived electric field (Figure 3, Equation 5).The upper and lower limits of the integration are chosen such that the result of the integral gives a velocity equation as a function of time.Next, the solution of integral is scaled to accelerate the particle up to 94% of the speed of light, i.e., 280,000 km/s, and thereby prevents the particle from exceeding the speed of light in vacuum, c.Finally, the particle velocity function should decay to the thermal velocity.Afterward, the particle should either remain constant at a thermal velocity or slowly decay to 0 m/s and remain stationary at an infinite time for the particle to stop radiating.A stationary or constant velocity at infinite time is important because Equation 13 requires an integral with limits between ∓∞, and the particle cannot radiate for an infinite time.The derived velocity function (Equation 9) presented in Figure 4 accelerates the particle from a thermal velocity (92.493 km/s) up to 94% of c within ∼12.427 ns and then decelerates back to a thermal velocity.The whole process takes ∼24.854 ns.
The velocity function is the integral of the electric field function with respect to time, multiplied with the constant charged particle properties, which are elementary charge e, charge number z, and the inverse of accelerated electron mass Where t o = 1 s.The upper limit was set to be infinite in order to prevent the lower boundary of the integral from exceeding the upper limit.The lower limit involves a time parameter to prevent the indefinite integral from being a definite integral in order to create a velocity function.The absolute time is used here to create a symmetric velocity behavior in retarded time to explain the particle motion and also to satisfy the integral limits between ∓∞ in Equation 13.The factorial was incorporated in order to prevent the electric field from becoming singular at t = 0 s.
YÜCEMÖZ AND FÜLLEKRUG 10.1029/2020JD033204 8 of 23 Lightning discharges exhibit a sharp rise and slowly decaying electric field.The peak value of the electric field was chosen to be ∼3 MV/m, which is the approximate electric field for the conventional dielectric breakdown of air.For a RREA to develop, the electrons need to be relativistic such that bremsstrahlung occurs.In this case, the leader tip electric field is on the order of ∼26 MV/m at 101 kPa of atmospheric pressure (Babich et al., 2015).

(b)
The derived electric field from 0 s to ∞ is plotted on logarithmic time scale.It is important that the leader tip electric field is 0 V/m at infinity.Any function used to describe the leader tip electric field that never goes to zero (i.e., an exponential) integrated up to infinity would cause a particle to exceed the speed of light.
Separating integral (Equation 6) into the two terms of the electric field function (Equation 5) gives The factorial of time only works for an integer number of time values.However, this is only a problem when the velocity function is used with non-integer values.To approximate velocities with a factorial of non-integer time values, Stirling's approximation can be used.The velocity function can be written with the substitution of all the constant values that do not require to be changed in order to scale.Also, to preserve the physical meaning of the characteristic lightning leader tip electric field feature with a sharp rise and a slow decay taking place in the atmosphere.9) of the derived electric field.(a) There is no negative chronological time.However, negative time has a meaning in physics that it describes an off the record process of an experimentally observed electromagnetic pulse.Negative time is the retarded time in Liénard-Wiechert potentials that describes the radiation emission by the acceleration of a single particle.Retarded time in Liénard-Wiechert potentials is also the reason why Equation 13 has integral limits between −∞ and ∞.Negative time describes the duration of the process until the pulse detection time of the observer.Therefore, starting from ∼−10 ns up to to ∼10 ns, the particle gains velocity through an acceleration by the external leader tip electric field, and as it accelerates, it radiates electromagnetic radiation.When the particle velocity is peaking around t = 0 s, the observer starts receiving a signal and can measure the pulse as the process of particle acceleration, and the propagation of the electromagnetic wave from the source to an observer requires some time -known as the retarded, or negative, time.(b) The electron's velocity on an infinite time scale displays the end of a particle acceleration outside the range of ∼∓11 ns until ∓∞.The peak velocity remains the same at 94% of the speed of light.
As the particle has a short mean free path with a short mean free time ∼ns, the time scale in the velocity function ( 8) was scaled accordingly to ns with s ftv = 1 × 10 9 to reflect atmospheric electrodynamic conditions.Previous to the scaling, as the time t is integrated up to ∞, the time scaling in the velocity function ( 8) has grown and it took ∼12.427s for the electron to reach 94% of the speed of light from its thermal velocity By knowing the particle velocity, the particle velocity equation (Equation 9) can be scaled directly by using the corresponding scaling factors for the required velocity.The mean free path (Equation 4) of a particle is a function of the total particle velocity.Although the derivative of the position vector (Equation 3) alone determines the velocity arising from the Coulomb force, Equation 3 is the total velocity when used for the mean free path (Equation 4).Hence, if the particle velocity is known, it can be equated to the derivative of the position vector on a spiral trajectory (Equation 3) to determine the parameters a (the relativistic mean free path correction parameter) and b (the interaction distance between an incoming and a target particle).This algebraic calculation transforms the velocity arising from the Coulomb field, that is, the derivative of the position vector in Equation 1, to an overall velocity, that is, the sum of the velocity arising from the Coulomb field and the velocity originating from the external lighting leader tip electric field for a particle on a spiral trajectory.
As can be seen from Figure 4, the above assumptions help to construct a velocity function (Equation 9).The complete radiation model is based on the fact that the particle reaches approximately 94% of the speed of light.Then the particle decelerates back to a thermal velocity as bremsstrahlung is the braking radiation due to obstacle particles or atoms, and as the source electric field from the leader tip starts to decay.
The relative particle velocity is .
where t is the time in s and the range of the time changes in the velocity function (Equation 9), hence in β (Equation 11) function.The time range in velocity function differs from the electric field function given in Equation 5 as a result of the integration process.The time range of the velocity function is −∞ ≤ t ≤ ∞.An important information to note about the indicated range of time t values separately in both, the scaled velocity function (Equation 9) and the electric field function (Equation 5), is that they are specific for this scaled velocity and non-scaled electric field function separately.There is no time parameter in any other equation nor in the final equation due to the nature of the definite integral of time.However, as indicated in Equation 13, the definite integral between the limits of −∞ and ∞ of time results in a time t parameter to vanish in the final equation (Equation 36).Therefore, the range of time t values presented in Sections 4 and 5 can be neglected for the final Equation 36.
Before moving to the next section, the concepts explained so far make up the base of the complete theory.The next section uses the information of Sections 2, 4, and 5 to construct the core of the theory and to derive the novel generic radiation intensity pattern of a single particle unique to the bremsstrahlung process.Figure 5 shows how this theory fits into practical observations.YÜCEMÖZ AND FÜLLEKRUG  This section introduces collision cross-sections as a part of explaining complete radiation pattern.Further information on collision cross-section and mathematical definitions consisting of form factors can be found in supplemental information file in the form of Text S2 (Punjabi & Perdrisat, 2014).Collision cross-section and form factors are very important in defining the final form of the radiation patterns at high frequency (Perdrisat et al., 2007).Form factors are not used in the final equation and in simulations of radiation patterns.For high-frequency emissions, the details of collision and interaction are crucial.They are accounted for by adding a cross-section of a particle of interest.
Hence, it is assumed that Dimensionally, χ = I × σ, where, σ in m 2 is the Coulomb scattering cross section, χ in Js −1 m 2 is the radiation cross section, and I in Js −1 the radiation intensity.As a result of Liénard-Wiechert potentials, the radiation emitted by a single charged particle (Jackson, 1999, p. 675) is given by where In addition, z is the dimensionless charge number, c is the speed of light, ϵ 0 is the permittivity of free space, ω′ is the angular frequency of the emitted wave in rad/s in the particle's frame of reference and t is the time in s.
The radiation unit vector n, the position vector r(t) in Equation 1, the particle motion, the position of an observer, and the overall theoretical concept about the bremsstrahlung radiation emission modeling that the equations are built upon are shown in Figure 2.
As n is a unit vector in the direction of the emitted radiation, the integrand of the integral in Equation 13can be re-written in the sinusoidal form of the cross-product as or in the sinusoidal form of dot product as  shows how the particle's velocity changes as a function of time.The particle is still allowed to propagate to positive time values.However, the starting time is from negative time values (retarded time).Negative time is required to reflect the reality that when the signal is measured by an observer at t = 0 s, it has to travel some time from the source to the observer.The acceleration process that causes the emission also requires some time to take place.The observer in chronological time does not observe all of these processes.Hence, they occur in negative retarded time.(c) sketches when the observer would receive the signal.It is plotted using the Dirac Delta function to sketch when the observer would receive the signal, which does not reflect the actual timing.The first signal could also be received at t = 0 s, however not in retarded time as it would then contradict the observations where there is no negative time. .

ds (17)
The exponential term can be simplified such that Further simplification of the common parameters results in where e −iωs can be neglected as there is no 1 c term to reduce the rate of exponential decay.Therefore, this term goes to zero quickly with time such that Finally, the velocity Equation 9 converging to 0 m/s at ∓∞ (see Figure 4b) enables the integral in Equation 21 to be evaluated.Convergence is due to the nature of derived electric field function, which goes to zero at infinite time and does not accelerate the particle any more.Furthermore, the use of functions such as Planck's radiation curve, Heidler current (Heidler et al., 2013), or Poisson distribution function with a variable transform to represent electric field would not be beneficial in deriving a velocity function as they diverge due to their exponential nature.The solution to the integral is approximated by the use of a formula in the book of Fourier transforms of exponential functions (Bateman et al., 1954, p. 121, Equation 23).The solution of the integral requires index ν in the variable (ix) ν that multiplies the exponential function in Exponential Fourier Transform of function f(x) = f(s) to be ν > −1 (Bateman et al., 1954, p. 121, Equation 23).The derivation of the velocity (Equation 9) results in the beta function (Equation 11) satisfying this requirement.
As the integrand of the integral in Equation 14multiplies with the β function consisting of 2 terms given in Equation 11, the integral can be separated and solved individually for each term such that all terms can be added together to give the final result.
Hence, the approximated solution of the integral ζ for the first term of the scaled velocity function ( 11) The approximated solution of the integral ζ for the second term of the scaled velocity function ( 11) where .
Therefore, the final radiation pattern for both high and low-frequency emission is the sum of the two terms where R is the bremsstrahlung asymmetry control parameter, which plays a crucial role in distinguishing the parameter y from α, which is a requirement in order to be able to approximate the solution to an integral (dimensionless).D v (z) is the parabolic cylinder function.
The parabolic cylinder function is given by (Whittaker & Watson, 1927, p. 347 where, W v/2+1/4,1/4 (1/2z 2 ) is a Whittaker function (Whittaker & Watson, 1927, p. 346) and where, M κ,μ (1/2z 2 ) is another Whittaker function (Kiyosi Ito and The Mathematical Society of Japan, 1993;Whittaker & Watson, 1927, p. 347) and YÜCEMÖZ AND FÜLLEKRUG 10.1029/2020JD033204 For the second term of Equation 27with Whittaker M function of negative μ where 1 F 1 (1/2 + μ − κ; 2μ + 1; z) is a confluent hypergeometric function of the first kind (Abramowitz & Stegun, 1972) and where is the confluent hypergeometric function of the second kind, named Kummer's U function such that By direct comparison of Equation 26with 25 and Equation 26with 27, the constant quantities z, κ and μ are defined as Finally, in order to complete Equation 25, the observation angles (θ n,β and θ n,r(t) ) need to be reduced from two to one as they are related to each other.
Currently, the radiation pattern is observed from cos(θ n,r(t) ) and sin(θ n,β ).In addition, cos(θ n,r(t) ) coming from substitute parameter α and y in Equation 25and sin(θ n,β ) being a variable parameter in Equation 25defining the radiation around the particle, also known as the solid angle.In addition, by looking at Figure 2b, it can be seen that the velocity vector, which is the derivative of the position vector r(t) is always perpendicular to the position vector where radiation unit vector pointing in the direction of the emitted radiation makes up the hypotenuse of the formed rectangular triangle.Hence, the angles of the rectangular triangle can be equated to each other.
The velocity function defined in Equation 9provides the particle speed attained by the acceleration due to the external leader tip electric field, and it is scalar.The derivative of the position vector r(t) (Equation 1) provides the velocity caused by the Coulomb interaction of an incoming and target particle, which also determines the curved trajectory of the bremsstrahlung given by the position vector (Equation 1).Hence, the particle's overall velocity moving perpendicular to the position vector r(t) is the sum of the Coulomb velocity arising from the derivative of the position vector r(t) plus the scalar speed attained by the external leader tip electric field given in Equation 9. Equation 9 is not the velocity, but the speed.It only provides magnitude but gives no information about the direction of the particle's motion.The direction of the particle's motion is already provided by the position vector r(t) given by Equation 1.As a result and n n r t sin cos (34)

Dynamic Bremsstrahlung -Special Relativity and Doppler Effect
Now the Doppler shift is added to the emitted radiation.
Figure 6 displays a stationary laboratory frame ω where the relativistic effects and the Doppler shift is observed, particle's moving frame of reference ω′ and the relative velocity cβS SpecialR between the observer and the two different frames that determines what an observer would detect.In addition, S SpecialR (dimensionless) is a scaling factor of the relative velocity between the two different frames of references.In particular, S SpecialR = 1 means the observer is stationary, and the particle is propagating at its velocity.As the S SpecialR approaches zero, the relative velocity between the observer and the particle falls, meaning that the observer is getting closer to the particle's velocity and would start to observe or detect a non-relativistic dipole radiation pattern.
The equation for translating the angular frequency of the emitted wave into the laboratory frame for the Doppler shift is given by (Jackson, 1999, p. 720, Equation 15.40). .
As shown in Figure 2b, θ n,β is the angle between the emitted radiation unit vector and the particle velocity vector.
Substituting Equations 34 and 35 into the Equation 25gives YÜCEMÖZ AND FÜLLEKRUG 10.1029/2020JD033204 16 of 23 .In addition, S SpecialR = 0 meaning the observer is moving with the particle at the same velocity at all times, hence would see no radiation emission as the particle is not accelerating.As S SpecialR → 0, the observer would see dipole radiation pattern.As S SpecialR → 1, the observer would approach rest and would see the particle propagating at its own velocity defined by Equation 9with the corresponding radiation pattern.
where the Lorentz factor γ is γ =   2 1 1 (dimensionless), and ω is the received angular frequency in the laboratory (stationary) frame of reference in rad/s.In addition, the new definition of , and the new definition of y is The scaling factor s f in Equation 36 is carried over from the definition of the electric field Equation 5.As a result of the integration process, s f scales the magnitude of the emitted radiation, which is directly linked with its previous purpose of scaling the electric field, which defines the magnitude of the emitted radiation.The scaling factor s f can be approximated at high frequency with a scalar equation that gives the magnitude of the emitted high-frequency radiation from a single particle in an instantaneous circular motion (Jackson, 1999, p. 679, Equation 14.84). Hence, and the critical angular frequency (Jackson, 1999, p. 679, Equation 14.81) is where p is the radius of curvature of instantaneous circular motion in m, ω c is the critical angular boundary frequency.Beyond critical angular boundary frequency, the emitted radiation would have a minimal value in all directions such that it can be neglected.

Results
This section presents predictions of radiation patterns of a single electron accelerated under an external lightning leader tip electric field using the derived Equation 36.
Radiation patterns change as a result of changing particle velocity, and acceleration as a particle get relativistic.As the acceleration of a charged particle increases, it emits higher intensity, hence higher frequency radiation.Therefore, the input parameter, emitted radiation frequency represents particle acceleration.To simulate relativistic and non-relativistic radiation patterns, the range of input velocity and frequency values were chosen and tabulated in Table 1.
The new radiation patterns reveal that as well as forward peaking, there are also backward peaking lobes as the particle gains speed.The magnitude of the acceleration determines the magnitude of the radiation intensity.The horizontal axis is the radiation intensity per solid angle per emitted frequency range in eight different radiation patterns displayed in Figure 7.Although there is overall more energy in the forward peaking lobes, the peak radiation intensity in the backward direction is slightly higher in patterns 5, 6, and 7 as they are more horizontally projected.However, the total radiation intensity is higher in the forward direction as the peaking lobes have larger beam areas due to the Doppler effect's presence.
The peak intensity of the backward lobes starts to decay as a result of the Doppler shift from pattern seven onwards.The Doppler shift effect can be observed from the consideration of both radiation length in the horizontal direction (as the horizontal axis is a measure of the radiation intensity per solid angle, i.e., more horizontally projected) and the total area of the emitted radiation pattern.The larger the area of the radiation pattern, the more the beam extends in the horizontal direction or, when it is more projected toward the horizontal axis, the more energy the emitted wave has.An increasing particle velocity translates into increasing radiation frequency.Therefore an increase in energy in the forward direction of the particle lowers the frequency and hence lowers the energy at the backward direction of the particle due to the Doppler effect.For example, in the eighth radiation pattern, which is at ∼94% of the speed of light and has 1 GHz frequency, increasing the frequency of the observed radiation leads to an increase in the Doppler effect that reduces the frequency, hence the energy (E = hf) of the beam at the back of the particle.As a consequence, backward radiation lobes start to shrink such that overall radiation pattern looks like a basic forward peaking radiation pattern (the pattern 8).
Radiation patterns on the polar plot have to start from 2π rad until 0 rad to show the particle's complete radiation.Progression toward the zero is because a particle starts from the retarded time and radiates until the chronological time t = 0 s and beyond, as shown in Figures 2 and 5.If the polar plot were made with time as a representation of the solid angle, the plot would have started from ∼−10 ns (when the particle starts accelerating, see Figures 2 and 5) until 0 s.In addition, this is the same as starting from positive time and ending at t = 0 s.Therefore, the final result is independent of positive or negative start time as long as the time flow is toward zero seconds.The independence of the sign is caused by the absolute time and symmetry of the velocity function.Progression toward zero seconds is caused by the initial start time being retarded negative time.
The derived mathematical model (Equation 36) is specifically for the bremsstrahlung process.Unlike many of the other radiation emission processes of charged particles (i.e., from linear acceleration or crossing the boundary between two different dielectric media), the bremsstrahlung process affects the shape of the emitted radiation.This bremsstrahlung effect causes an asymmetry of the emitted radiation about the particle's velocity vector or, in other terms, the direction of motion.To understand this effect, we can compare an electron to a car that travels in the dark with the headlight turned on.The headlights are the emitted electromagnetic radiation by the car.When the car gets into the bend, like the bremsstrahlung process of an electron, an observer outside the car can immediately tell the radiation shape of the headlight would be asymmetric by looking at the reflections on the road compared to the case when the car follows a straight path.This novel effect is clearly shown in supporting information in Figure S1, which displays a real-life example of visible asymmetry of the car headlights (electromagnetic radiation) on a bend.
YÜCEMÖZ AND FÜLLEKRUG Note.Emitted radiation frequency by a single particle is an input parameter to the final Equation 36.It plays a crucial role in delivering valuable information and determining the particle's total acceleration, hence the total velocity is indirect as there is no input acceleration parameter in the final Equation 36.In addition, β function determines the particle's total velocity arising significantly from the external lightning leader tip electric field.The remaining contribution to the overall velocity comes from the Coulomb force between an incoming and target particle, where the target particle causes an incoming particle to follow and covers some arc length of the spiral trajectory described in the position vector (Equation 1) during its mean free time of the collision.  .(g) A particle reaching relativistic speed attains a focused beam in both forward and backward direction.(h) A particle reaching to ultra-relativistic speed experiences the take over of the Doppler effect that minimizes the lowfrequency radiation in the backward direction and emits almost fully in the forward direction.

Discussion
First, the created model describing radiation patterns specifically for the bremsstrahlung process is modeled to be generic and applicable to all bremsstrahlung process on earth, outer space, or other planets.particle or ion via Coulomb field.In addition, the incoming particle is assumed to follow a spiral particle trajectory.This assumption is based on the in-balance of centripetal (Coulomb field between two particles + Earth's external magnetic field) and centrifugal (external leader tip electric field) forces acting on the incoming particle when it interacts with a target particle via Coulomb field.Any in-balance in centripetal and centrifugal force acting on an incoming particle that follows a curved trajectory would result in a spiral particle trajectory.Hence, the choice of spiral trajectory considers the effects of the Earth's external magnetic field, B on an incoming particle trajectory in an indirect way.
Particle speed defined in Equation 9 relies on the ansatz of a sharp rise and slow decay form of leader tip electric field.This ansatz originates and is backed up by the asymmetric γ-ray bursts (GRBs) measurement that reveals the internal charge structure of the leader tip described in the introduction section, in the second paragraph.
Although acceleration is modeled from the electric field that assumes the same behavior as the leader tip electric field observed in the atmosphere, it only provides a magnitude component of the electric field acceleration.Independent input parameters such as particle trajectory, the bremsstrahlung asymmetry, R, and the Doppler effect are fundamental physics where the particle trajectory, and the bremsstrahlung asymmetry, R are specific to the bremsstrahlung process.These input parameters enable the created model to apply to any bremsstrahlung processes.
The bremsstrahlung process has a continuous electromagnetic spectrum meaning it is capable of explaining both ionizing and low-frequency radiations.The bremsstrahlung radiation patterns presented in Figure 7 are unique to a bremsstrahlung process.They are always the same in geometry with forward-backward radiation patterns at relativistic particle speeds and all the asymmetries.However, when the model is applied to any other bremsstrahlung processes other than TGFs, only radiation intensities are subject to change but not the radiation shapes.
Overall, the eight different radiation patterns plotted in Figure 7 are not the complete radiation patterns, especially for high-frequency radiation.As mentioned in the previous section, a high-frequency radiation pattern is sensitive to the interaction process described by scattering cross-sections that an incoming bremsstrahlung electron experiences against target particles.Moreover, plotted radiation patterns are the second-order derivative term of Equation 12, and it misses the information of interaction described by Coulomb scattering cross-section.In future work, a complete radiation pattern will be calculated using the complete Equations 12 and 36, and Equation 1 provided in the supplemental information, S2 with appropriate form factors, and a dimensional analysis to describe high-frequency emissions.Furthermore, the first part of the derived equation demonstrates a success as Figure 7a (patterns 1) and 7h (pattern 8), display the common expected radiation patterns such as low-frequency dipole radiation pattern and relativistic forward peaking radiation patterns (Jackson, 1999, p. 669, Figure 14.4).
The dipole radiation pattern 1 in Figure 7a displays perfect symmetry about both the velocity axis and the axis perpendicular to the velocity vector.In the absence of a Doppler and bremsstrahlung asymmetric radiation effect, the four maxima of the relativistic radiation pattern are symmetric about both the velocity axis and the axis perpendicular to the velocity.Therefore, following the statement of Noether's theorem, conservation of energy arising from the symmetry in the radiation lobes about the two axes mentioned above, causes the dipole radiation pattern 4 and 5 (in Figures 7d and 7e) to collapse and form 4 symmetric maxima.Initially, there were symmetries about two axes, and as the particle gains up speed, the number of symmetric axes should still be preserved.Hence, if there is only a forward peaking radiation with none in the backward direction, one of the symmetry is broken (about axes perpendicular to the velocity).The idea of preserving the initial number of symmetric axes is quite powerful in understanding why there are four maxima (two forward peakings and two backward peakings) and why the dipole collapses to form 4 maxima of the radiation pattern.All the radiation patterns in Figure 7 incorporate the Doppler effect and the bremsstrahlung trajectory asymmetry.These two effects cause the radiation pattern to be asymmetric.However, the emitted total radiation energy is still conserved compared between the presence and absence of radiation pattern asymmetry causing physical effects (The bremsstrahlung trajectory and the Doppler effect) at any given particle velocity.In other words, the total energy radiated at a given velocity is the same regardless of the Doppler and the bremsstrahlung trajectory effect.While these two effects increase the YÜCEMÖZ AND FÜLLEKRUG 10.1029/2020JD033204 radiation beam intensity in one direction, they also decrease the radiation intensity in other directions by the same amount such that the overall energy radiated by the particle at a given speed remains the same.
On the other hand, the reported radiation asymmetry of the bremsstrahlung (R-parameter) about the velocity vector due to the particle's curved trajectory was found to be existing at relativistic speeds when forward and backward peaking radiation patterns are present.This bremsstrahlung radiation asymmetry (R-parameter) we get is a novel finding.At the non-relativistic speeds, when radiation has a dipole pattern, asymmetry of the bremsstrahlung was found to be absent.This absence indicates that asymmetry is associated with being a physical effect rather than a mathematical artifact.The absence of the bremsstrahlung radiation asymmetry at a non-relativistic dipole radiation pattern of a particle following a curved trajectory can again be compared to a car having its headlights on each side of the door radiating perpendicular with respect to the velocity vector of the car.In this case, one can see that two emitted radiation beams would not be asymmetric as two beams are in the opposite direction to each other and are independent of the curvature of the trajectory as headlights radiate perpendicular to the trajectory.
Asymmetries are important because they alter the radiation intensity of the emitted radiation.Asymmetries increase the radiation intensity in one direction and decrease it in the other direction.The bremsstrahlung radiation pattern of a single particle with involved asymmetries can be used with a Monte Carlo simulation for multiple particles to predict radiation intensity map across an area where lightning discharge takes place.The radiation intensity map is the expected radiation power at each different observation location around the area of lightning discharge location.The predicted radiation intensity map allows improving the lightning detection and location systems.It helps to find the location of peak radiation intensity, hence contributes to the development of better radiation protection systems.

Summary
In summary, the radiation patterns were found to be peaking in a backward direction as well as an already known forward direction.In addition, four maxima (2 in the forward direction and 2 in the backward direction) were found to be due to the conservation of symmetry axes arising from the initial dipole radiation pattern.With the introduction of the Doppler effect, peaking lobes of radiation in forward and backward directions were found to be asymmetric about an axis perpendicular to the particle's velocity vector.Moreover, the novel second asymmetry of peaking lobes with respect to particle's velocity vector was found to be unique to the bremsstrahlung due to the particle following a curved trajectory.Finally, it is also found that low-frequency radiation peaks in the backward direction, whereas high-frequency peaks in the forward direction.M. Yücemöz would like to thank to his supervisor M. Füllekrug for all the opportunities, suggestions, guidance, reviews, and support throughout his PhD study.The Engineering and Physical Sciences Research Council (EPSRC) and MetOffice sponsor M. Yücemöz's PhD project under contract numbers EG-EE1239 and EG-EE1077.M. Füllekrug acknowledges support from the Natural Environment Research Council (NERC) for under grants NE/L012669/1 and NE/H024921/1.M. Yücemöz wishes to thank Dr. Adrian Hill for mathematical support in helping to solve integral with a divergence problem and his family for their support and good wishes.

Figure 1 .
Figure 1.(a) The flow chart shows the general structure of the theory as it is developed throughout the text.(b) Explains how particle's curved trajectory caused by Coulomb field changes as a result of the external non-uniform leader tip electric field and how it is accounted for in the theoretical model.The yellow lines are the external electric field, black arrows are the tangential Coulomb velocity, thick blue arrows are the centripetal acceleration.Finally, green lines are the acceleration caused by the external electric field.As can be seen in region B, one component of the leader tip electric field acceleration is in the direction of tangential Coulomb velocity.Whereas, the other component is in the opposite direction of the centripetal acceleration.Leader tip electric field acceleration in the opposite direction of the centripetal acceleration of the Coulomb field only changes the radius of curvature of the trajectory.Hence, it can be accounted with parameter "b" in the trajectory Equation1, defining the radius of curvature.The acceleration of the leader tip electric field tangent to the trajectory is accounted with the leader tip electric field velocity Equation8.Finally, region A explains how the bremsstrahlung process ends, when the leader tip electric field acceleration component opposite to centripetal acceleration (contributing to the centrifugal acceleration) is larger than the centripetal acceleration of the target particle.The condition at region A allows an incoming particle to escape by overcoming the Coulomb force.

Figure 2 .
Figure 2. (a) Trajectory of the bremsstrahlung electrons given by Equation 1 in a polar coordinate system and radiation emissions by the change in velocity over time by a Coulomb force of other charges.O represents the target particle that defines the electron's spiral trajectory due to the Coulomb force (Equation1).The tangential component of the electron acceleration and corresponding velocity vectors are displayed with black arrows and are tangential to the spiral trajectory (red line) and perpendicular to the position vector r(t).R(t r ) is the distance between the accelerated electron and the observer, which is a function of both retarded and chronological time (t r , t).P is the position of an observer.S is the Poynting vector, or radiant energy flux, which determines the direction of the energy flow per area of an emitted electromagnetic wave.The dimensionless unit vector n points in the direction of the Poynting vector 

Figure 3 .
Figure 3. Derived electric field in Equation 5. (a) The leader tip electric field is modeled to fit observed electric field characteristics of lightning discharges.Lightning discharges exhibit a sharp rise and slowly decaying electric field.The peak value of the electric field was chosen to be ∼3 MV/m, which is the approximate electric field for the conventional dielectric breakdown of air.For a RREA to develop, the electrons need to be relativistic such that bremsstrahlung occurs.In this case, the leader tip electric field is on the order of ∼26 MV/m at 101 kPa of atmospheric pressure(Babich et al., 2015).(b) The derived electric field from 0 s to ∞ is plotted on logarithmic time scale.It is important that the leader tip electric field is 0 V/m at infinity.Any function used to describe the leader tip electric field that never goes to zero (i.e., an exponential) integrated up to infinity would cause a particle to exceed the speed of light.

Figure 4 .
Figure 4. Velocity function (Equation9) of the derived electric field.(a) There is no negative chronological time.However, negative time has a meaning in physics that it describes an off the record process of an experimentally observed electromagnetic pulse.Negative time is the retarded time in Liénard-Wiechert potentials that describes the radiation emission by the acceleration of a single particle.Retarded time in Liénard-Wiechert potentials is also the reason why Equation 13 has integral limits between −∞ and ∞.Negative time describes the duration of the process until the pulse detection time of the observer.Therefore, starting from ∼−10 ns up to to ∼10 ns, the particle gains velocity through an acceleration by the external leader tip electric field, and as it accelerates, it radiates electromagnetic radiation.When the particle velocity is peaking around t = 0 s, the observer starts receiving a signal and can measure the pulse as the process of particle acceleration, and the propagation of the electromagnetic wave from the source to an observer requires some time -known as the retarded, or negative, time.(b) The electron's velocity on an infinite time scale displays the end of a particle acceleration outside the range of ∼∓11 ns until ∓∞.The peak velocity remains the same at 94% of the speed of light.
493 km/s at an example ambient air temperature of 188.15 K, when considering the mesosphere at an altitude of 80 km.Similarly, as the time t is integrated up to ∞, the particle attained a large velocity with the multiplication of the constant term e ez m in front of the integral.This is solved with a new velocity scaling factor s fv = 8.19 × 10 −11 .

Figure 5 .
Figure5.Complete order of events from initial particle acceleration until pulse detection by an observer.(a) The yellow line indicates the particle's starting position.It starts accelerating from −10 ns until t = 0 s as a result of an external leader tip electric field and Coulomb electric field (Bremsstrahlung) of the target particle in the atmosphere.(b) shows how the particle's velocity changes as a function of time.The particle is still allowed to propagate to positive time values.However, the starting time is from negative time values (retarded time).Negative time is required to reflect the reality that when the signal is measured by an observer at t = 0 s, it has to travel some time from the source to the observer.The acceleration process that causes the emission also requires some time to take place.The observer in chronological time does not observe all of these processes.Hence, they occur in negative retarded time.(c) sketches when the observer would receive the signal.It is plotted using the Dirac Delta function to sketch when the observer would receive the signal, which does not reflect the actual timing.The first signal could also be received at t = 0 s, however not in retarded time as it would then contradict the observations where there is no negative time.
from t to s requires the Jacobian.Assuming that s = −t it follows that 

(
20)Bringing the whole integral in Equation 16 together with the simplified exponential gives

Figure 6 .
Figure 6.Particle radiation at two different frames of references.Laboratory (stationary) frame of reference ω and moving frame of reference ω′ with the velocity ( ) SpecialR v t c S c.In addition, S SpecialR = 0 meaning the observer is as ratio with the speed of light β and emitted radiation frequency(Hz) Figure8displays the effects of the bremsstrahlung radiation asymmetry control quantity R. The bremsstrahlung asymmetry R depends on the radius of the curvature of the incoming particle trajectory undergoing the bremsstrahlung process.In addition, the bremsstrahlung asymmetry R increases as the radius of curvature of the bremsstrahlung trajectory decreases.

Figure 7 .
Figure 7.The Radiation patterns emitted by anti-clockwize rotating charged particle -bremsstrahlung process.Plot is in Polar co-ordinates.Horizontal axis gives the radiation intensity per Solid angle, Ω, per emitted angular radiation frequency, ω.In addition, angle of the Polar plot is the Solid angle, Ω. ((a)-(f)) The radiation pattern connecting dipole with a forward-backward peaking radiation pattern is called the transition pattern.The radiation patterns and formula explaining the radiation patterns is not complete if it does not demonstrate the transition patterns.The transition patterns explain how the particle attains relativistic forward-backward peaking lobes.The radiation pattern of a particle during the bremsstrahlung, starting from a low velocity up until a relativistic velocity starts from the dipole radiation pattern and exhibits the forward and backward peaking radiation pattern.The transition is demonstrated in radiation patterns 4 and 5, where the dipole collapses to form four maxima.The values used for plotting are: mean free time τ = 30 μs, number of charges z = 1, a = 100 μm, b = 1 nm (a and b are related to mean free path), s ft = 1, s f = 1, S SpecialR = 1, velocity-time scaling factor s ftv = 1 × 10 9 and velocity scaling factor s fv = 8.19 × 10 −11 .Finally, the bremsstrahlung asymmetry is R = 1/8.In addition,   1 1 9 3 R .(g) A particle reaching relativistic speed attains a focused beam in As mentioned in the third paragraph of the introduction, TGFs are believed to originate as a result of the bremsstrahlung process of RREA driven by external leader tip electric field.Hence, the theoretical model assumes a single particle accelerated under an external leader tip electric field interacting with a