Neutral Density Variations Research Articles

SubAuroral Red Arcs Generated by Inner Magnetospheric Heat Flux and by SubAuroral Polarization Streams

AbstractSubauroral red (SAR) arcs are commonly observed ionospheric red line emissions. They are usually attributed to subauroral electron heating by inner magnetospheric heat flux (IMHF). However, the role of IMHF in changing the ionosphere‐thermosphere (IT) still remains elusive. We conduct controlled numerical experiments with the Thermosphere‐Ionosphere Electrodynamic General Circulation Model (TIEGCM). Coulomb collisional heat flux derived with the Comprehensive Inner Magnetosphere Ionosphere (CIMI) model and empirical subauroral polarization streams (SAPS) are implemented in TIEGCM. The heat flux causes electron temperature enhancement, electron density depletion, and consequently SAR arcs formed in the dusk‐to‐midnight subauroral ionosphere region. SAPS cause more substantial plasma and neutral heating and plasma density variations in a broader region. The maximum enhancement of subauroral red line emission rate is comparable to that caused by the heat flux. However, the visibility of SAR arcs also depends on the relative enhancement to the background brightness.

The Impact of Polar Vortex Strength on the Longitudinal Structure of the Noontime Mid-Latitude Ionosphere and Thermosphere

Ground-based ionospheric, CHAMP/STAR, and GOCE satellite neutral density ρ observations under deep solar minimum conditions were used to find whether there is a dependence of longitudinal variations on polar vortex strength. Ionospheric stations at fixed-dipole geomagnetic Φ ≈ 38° and geographic φ ≈ 40°N latitudes located in ‘near-pole’ and ‘far-from-pole’ longitudinal sectors were used in the analysis. No significant longitudinal NmF2 (electron concentration in the F2-layer maximum) dependence on the polar vortex strength was revealed. Geomagnetic control was shown to be responsible for the observed longitudinal NmF2 variations. Satellite-observed longitudinal variations in neutral density did not show any visible reaction to the polar vortex strength. However, the impact of sudden stratospheric warming (SSW) on the upper atmosphere is strong enough to change the neutral density longitudinal distribution. The impact of SSW shows a global occurrence and ‘works’ within 3–5 days in geographic coordinates in the vicinity of the SSW peak. Atomic oxygen values retrieved under ‘strong’ and ‘weak’ polar vortex strengths confirm the results obtained on longitudinal variations in NmF2 and ρ. In conclusion, no visible effects related to ‘strong’ or ‘weak’ polar vortex strengths have been revealed in either NmF2 or satellite neutral density longitudinal variations. Alternatively, such effects may be very small and therefore cannot be confirmed experimentally.

Temporally resolved relative krypton neutral density during breathing mode of a hall effect thruster recorded by TALIF

A key issue in the development of theory and models for plasma propulsion devices is to describe the instabilities and fluctuations of the devices. It has been widely recognized that many Hall effect thrusters (HETs) exhibit oscillations at frequencies in the range of ∼ 20 kHz. These ionization-related oscillations are generally referred to as Breathing Mode oscillations and have been the subject of considerable research. Here, for the first time, we report direct temporally resolved measurements of the ground state neutral density variation during the period of the oscillation. We used the laser-based Two-Photon Absorption Laser Induced Fluorescence (TALIF) technique to measure neutrals within the plume of a 1.5 kW HET operating on krypton (Kr). Our TALIF scheme employs a frequency-doubled, pulsed dye laser operating at ∼ 212 nm to probe ground state Kr atoms. A novel phase-binning approach is used to recover the time-dependent signal by assigning the timing of each collected TALIF signal (laser shot) relative to the phase of the discharge current. We find that the neutral density fluctuates quite strongly over the period of the oscillation, and that this fluctuation leads the current fluctuation as expected.

Large‐Scale Wave‐Driven Interactions and Plasma‐Neutral Coupling in the Low‐Latitude Ionosphere‐Thermosphere

AbstractThe plasma and neutral density variations, interactions and coupling processes within ±30° latitudes are examined concurrently by the DMSP‐F18 and Swarm‐C satellite during geomagnetically quiet years in 2020–2021. The wavenumber (WN) patterns are computed in the form of neutral and electron density for two altitudes and their latitudinal profiles are analyzed. We observe that the WN1 structure of the electron density has a significant seasonal dependence in the topside ionosphere and dominates all other structures but WN2 neutral density amplitude dominates all other structures in the middle thermosphere (∼440 km). Additionally, we analyze vertical‐temporal‐latitudinal tidal structures from the Climatological Tidal Model of the Thermosphere (CTMT) to find evidence for the modulation of the large‐scale waves (LSWs) neutral density structures. Through the examination of the in situ observational and modeling approaches, we show that the tidal contributors of WN structures obtained from CTMT can capture the influence of terrestrial sources on the WN structures of plasma‐neutral density and imprint the corresponding vertical coupling in the IT system. Correlation analysis reveals that the amplitudes of the WN1 and WN3 structures of electron density in topside ionosphere and those of neutral density in the middle thermosphere show intermittent but significant correlations with each other, unlike the WN2 and WN4 structures. This study provides new insights into the topside ionospheric response to wave driving in the lower atmosphere, which ultimately improves our capability to understand the interaction and vertical coupling of large‐scale structures, thereby advancing our predictive capabilities of space weather critical for satellite operations.

The Daytime Variations of Thermospheric Temperature and Neutral Density Over Beijing During Minor Geomagnetic Storm on 3–4 February 2022

AbstractOn 3 February 2022, 38 satellites launched by SpaceX re‐entered the atmosphere and were subsequently destroyed. An investigation found that a minor geomagnetic storm occurred on 3–4 February 2022 led to a neutral density enhancement and large atmospheric drag. To better understand the responses of the thermosphere to geomagnetic storms, the method proposed by Li et al. (2023, https://doi.org/10.1029/2022ja030988) was employed to extract exospheric temperature (Tex) from ionosonde electron density profiles (∼150–200 km) in Beijing (geolocation: 39.56°N; 116.2°E; geomagnetic location: 30.16°N; 172.08°W) station. The retrieved Tex was plugged into the NRLMSISE‐00 model to calculate the corresponding neutral density. Derived results showed a ∼2%–7% enhancement in Tex and a ∼15%–38% enhancement in neutral density at 430 km. The relative deviation in neutral density on the satellites’ orbital trajectory ranges from ∼10% (210 km) to ∼35% (500 km) on 3 February, and from ∼13% (210 km) to ∼60% (500 km) on 4 February. Furthermore, the neutral density reproduced the variations observed by the SWARM‐C satellite fairly well both on quiet and disturbed days. These results suggest that even a minor geomagnetic storm can cause significant changes in neutral temperature and neutral density at middle latitudes. Additionally, the application of our inversion method, combined with the global, long‐term and real‐time ionospheric observations from ionosondes, provides an opportunity to improve the capability of thermosphere forecasting and nowcasting.

Simulated long-term evolution of the thermosphere during the Holocene – Part 1: Neutral density and temperature

Abstract. In the previous work of Yue et al. (2022), the ionospheric evolution during the Holocene (9455 BC to 2015 AD) was comprehensively and carefully investigated for the first time using the Global Coupled Ionosphere-Thermosphere-Electrodynamics Model developed at the Institute of Geology and Geophysics, Chinese Academy of Sciences (GCITEM-IGGCAS), driven by realistic geomagnetic fields, CO2 levels, and solar activity derived from ancient media records and modern measurements. In this study, we further quantify the effects of the three drivers on thermospheric neutral density and temperature variations during the Holocene. We find that the oscillations of solar activity contribute more than 80 % of the thermospheric variability, while either CO2 or the geomagnetic field contributes less than 10 %. The effect of CO2 on the global mean neutral density and temperature is comparable to that of the geomagnetic field throughout the Holocene but is more significant after 1800 AD. In addition, thermospheric density and temperature show approximately linear variations with the dipole moment of the geomagnetic field, CO2, and F10.7, with only the linear growth rate associated with the geomagnetic field varying significantly in universal time and latitude. The increasing dipole moment and CO2 cool and contract the thermosphere, while solar activity has the opposite effect. The higher the altitude, the greater the influence of the three factors on the thermosphere. Different factors produce different seasonal variations in thermosphere changes. Furthermore, we predict that a 400 ppm increase in CO2 will result in a 50 %–70 % and 84–114 K reduction in global mean neutral density and temperature, respectively, which should directly affect the orbit and lifetime of spacecraft and space debris.

An unstable 0D model of ionization oscillations in Hall thruster plasmas

The breathing mode is an instability typical of Hall thrusters, which is characterized by oscillations of the discharge current with amplitude of the order of its mean value and frequency in the 5–30 kHz range. The strong link between this instability and the ionization processes is generally recognized. If, on one hand, 1D simulations have shown to be able to reproduce the breathing mode, on the other hand 0D models fell short in recovering self sustained oscillations, making it hard to identify the core physical mechanism governing their formation. In this work an original 0D model is presented and characterized by means of linear stability analysis and direct numerical integration. The electric field is allowed to vary in response to variations of the neutral density, acting on the ionization rate via the electron temperature and the ion dynamics. It is shown that the model is able to reproduce self-sustained oscillations with the typical characteristics of the breathing mode, even when fluctuations of the electron temperature are neglected. The stability of the model is strictly determined by the rigidity with which variations of neutral density reflect into variations of electron mobility.

Thermospheric Neutral Density Variation During the “SpaceX” Storm: Implications From Physics‐Based Whole Geospace Modeling

AbstractThe Starlink satellites launched on 3 February 2022 were lost before they fully arrived in their designated orbits. The loss was attributed to two moderate geomagnetic storms that occurred consecutively on 3–4 February. We investigate the thermospheric neutral mass density variation during these storms with the Multiscale Atmosphere‐Geospace Environment (MAGE) model, a first‐principles, fully coupled geospace model. Simulated neutral density enhancements are validated by Swarm satellite measurements at the altitude of 400–500 km. Comparison with standalone TIEGCM and empirical NRLMSIS 2.0 and DTM‐2013 models suggests better performance by MAGE in predicting the maximum density enhancement and resolving the gradual recovery process. Along the Starlink satellite orbit in the middle thermosphere (∼200 km altitude), MAGE predicts up to 150% density enhancement near the second storm peak while standalone TIEGCM, NRLMSIS 2.0, and DTM‐2013 suggest only ∼50% increase. MAGE also suggests altitudinal, longitudinal, and latitudinal variability of storm‐time percentage density enhancement due to height dependent Joule heating deposition per unit mass, thermospheric circulation changes, and traveling atmospheric disturbances. This study demonstrates that a moderate storm can cause substantial density enhancement in the middle thermosphere. Thermospheric mass density strongly depends on the strength, timing, and location of high‐latitude energy input, which cannot be fully reproduced with empirical models. A physics‐based, fully coupled geospace model that can accurately resolve the high‐latitude energy input and its variability is critical to modeling the dynamic response of thermospheric neutral density during storm time.

Importance of lower atmospheric forcing and magnetosphere-ionosphere coupling in simulating neutral density during the February 2016 geomagnetic storm

During geomagnetic storms a large amount of energy is transferred into the ionosphere-thermosphere (IT) system, leading to local and global changes in e.g., the dynamics, composition, and neutral density. The more steady energy from the lower atmosphere into the IT system is in general much smaller than the energy input from the magnetosphere, especially during geomagnetic storms, and therefore details of the lower atmosphere forcing are often neglected in storm time simulations. In this study we compare the neutral density observed by Swarm-C during the moderate geomagnetic storm of 31 January to 3 February 2016 with the Thermosphere-Ionosphere-Electrodynamics-GCM (TIEGCM) finding that the model can capture the observed large scale neutral density variations better in the southern than northern hemisphere. The importance of more realistic lower atmospheric (LB) variations as specified by the Whole Atmosphere Community Climate Model eXtended (WACCM-X) with specified dynamics (SD) is demonstrated by improving especially the northern hemisphere neutral density by up to 15% compared to using climatological LB forcing. Further analysis highlights the importance of the background atmospheric condition in facilitating hemispheric different neutral density changes in response to the LB perturbations. In comparison, employing observationally based field-aligned current (FAC) versus using an empirical model to describe magnetosphere-ionosphere (MI) coupling leads to an 7–20% improved northern hemisphere neutral density. The results highlight the importance of the lower atmospheric variations and high latitude forcing in simulating the absolute large scale neutral density especially the hemispheric differences. However, focusing on the storm time variation with respect to the quiescent time, the lower atmospheric influence is reduced to 1–1.5% improvement with respect to the total observed neutral density. The results provide some guidance on the importance of more realistic upper boundary forcing and lower atmospheric variations when modeling large scale, absolute and relative neutral density variations.

Observations of Modulation of Ion Flux in the Coma of Comet 67P/Churyumov‐Gerasimenko

On 6–8 June 2015, the Ion and Electron Sensor on board Rosetta observed keV-range water-group pickup ions arriving from the solar direction. Based on magnetic field intensification and variations, the appearance of the ions was likely to have been caused by a coronal mass ejection. During the 3-day period when Rosetta was 200 km from the comet, peak ion energy/charge (E/q) varied over a range from 50 eV to 1 keV in concert with neutral gas density variations caused by the rotation of the comet and its variable solar illumination. Thermal ion densities showed the same variations. The neutral density variations provided a unique opportunity to observe the repeated slowing of the solar wind by mass loading caused by charge exchange between energetic water-group ions and thermal water-group molecules. Such solar wind slowing was observed previously only by flyby missions that provided single events.

Proton loss of inner radiation belt during geomagnetic storm of 2018 based on CSES satellite observation* *Project supported by the Research Fund from the National Institute of Natural Hazards, Ministry of Emergency Management of China (Grant No. 2021-JBKY- 11), the National Natural Science Foundation of China (Grant Nos. 41904149 and 12173038), and the Stable Support Projects of Basic Scientific Research

The proton distribution in inner radiation belt is often affected by strong geomagnetic storm disturbance. Based on the data of the sun-synchronous CSES satellite, which carries with several high energy particle payloads and was launched in February 2018, we analyzed the extensive proton variations in the inner radiation belt in a wide energy range of 2 MeV– 220 MeV during 2018 major geomagnetic storm. The result indicates that the loss mechanism of protons was energy dependence which is consistent with some previous studies. For protons at low energy 2 MeV–20 MeV, the fluxes were decreased during main phase of the storm and did not come back quickly during the recovery phase, which is likely to be caused by Coulomb collision due to neutral atmosphere density variation. At higher energy 30 MeV–100 MeV, it was confirmed that the magnetic field line curvature scattering plays a significant role in the proton loss phenomenon during this storm. At highest energies > 100 MeV, the fluxes of protons kept a stable level and did not exhibit a significant loss during this storm.

Thermosphere modeling capabilities assessment: geomagnetic storms

The specification and prediction of density fluctuations in the thermosphere, especially during geomagnetic storms, is a key challenge for space weather observations and modeling. It is of great operational importance for tracking objects orbiting in near-Earth space. For low-Earth orbit, variations in neutral density represent the most important uncertainty for propagation and prediction of satellite orbits. An international conference in 2018 conducted under the auspices of the NASA Community Coordinated Modeling Center (CCMC) included a workshop on neutral density modeling, using both empirical and numerical methods, and resulted in the organization of an initial effort of model comparison and evaluation. Here, we present an updated metric for model assessment under geomagnetic storm conditions by dividing a storm in four phases with respect to the time of minimum Dst and then calculating the mean density ratios and standard deviations and correlations. Comparisons between three empirical (NRLMSISE-00, JB2008 and DTM2013) and two first-principles models (TIE-GCM and CTIPe) and neutral density data sets that include measurements by the CHAMP, GRACE, and GOCE satellites for 13 storms are presented. The models all show reduced performance during storms, notably much increased standard deviations, but DTM2013, JB2008 and CTIPe did not on average reveal a significant bias in the four phases of our metric. DTM2013 and TIE-GCM driven with the Weimer model achieved the best results taking the entire storm event into account, while NRLMSISE-00 systematically and significantly underestimates the storm densities. Numerical models are still catching up to empirical methods on a statistical basis, but as their drivers become more accurate and they become available at higher resolutions, they will surpass them in the foreseeable future.

Observations of Low‐Latitude Traveling Ionospheric Disturbances by a 630.0‐nm Airglow Imager and the CHAMP Satellite Over Indonesia

AbstractWe report the first comparison of ground and satellite measurements of low‐latitude traveling ionospheric disturbances (TIDs). Three TID events were simultaneously observed by a 630.0‐nm airglow imager and the CHAllenging Minisatellite Payload (CHAMP) satellite at Kototabang, Indonesia (geographic coordinates: 0.2° S, 100.3° E, geomagnetic latitude: 10.6° S). In 630.0‐nm airglow images of all three events, there are clear southward‐moving structures. Events 1 and 2 are a single pulse with horizontal scales of ∼500–1,000 km. Event 3 shows five wave fronts with a horizontal scale size of 500–1,000 km. All three TIDs are medium‐scale TIDs. Horizontal wavelengths of both airglow intensity at an average emission altitude of 250 km and CHAMP neutral density variations measured at 400 km are estimated by fitting a sinusoidal function to the observed data. The estimated horizontal wavelengths for airglow and neutral density data are 1,031 and 880 km for event 1 and 560 and 420 km for event 3, respectively. These values between airglow and CHAMP are comparable, suggesting both instruments are observing the same wave. For event 1, the CHAMP electron density mapped along the geomagnetic field line onto the airglow altitude does not show wave structure similar to the airglow variation. For events 2 and 3, the plasma density did not show wavy structures similar to the waves seen in the airglow image and CHAMP neutral density. These results suggest that the TIDs observed in airglow images are not caused by ionospheric plasma instability but by gravity waves in the thermosphere.

Quantifying the Storm Time Thermospheric Neutral Density Variations Using Model and Observations

AbstractAccurate determination of thermospheric neutral density holds crucial importance for satellite drag calculations. The problem is twofold and involves the correct estimation of the quiet time climatology and storm time variations. In this work, neutral density estimations from two empirical and three physics‐based models of the ionosphere‐thermosphere are compared with the neutral densities along the Challenging Micro‐Satellite Payload satellite track for six geomagnetic storms. Storm time variations are extracted from neutral density by (1) subtracting the mean difference between model and observation (bias), (2) setting climatological variations to zero, and (3) multiplying model data with the quiet time ratio between the model and observation. Several metrics are employed to evaluate the model performances. We find that the removal of bias or climatology reveals actual performance of the model in simulating the storm time variations. When bias is removed, depending on event and model, storm time errors in neutral density can decrease by an amount of 113% or can increase by an amount of 12% with respect to error in models with quiet time bias. It is shown that using only average and maximum values of neutral density to determine the model performances can be misleading since a model can estimate the averages fairly well but may not capture the maximum value or vice versa. Since each of the metrics used for determining model performances provides different aspects of the error, among these, we suggest employing mean absolute error, prediction efficiency, and normalized root mean square error together as a standard set of metrics for the neutral density.

Impact of nitric oxide, solar EUV and particle precipitation on thermospheric density decrease

Recent studies have revealed that the post-storm thermospheric mass density at a fixed altitude could be lower than the pre-storm mass density. This phenomenon has been attributed to an overcooling effect caused by enhanced infrared emissions from nitric oxide (NO). Here we report that thermospheric density decrease can take place under both storm and non-storm conditions. Relevant observations are presented during two superstormes, one intense storm, and one non-storm case. Thermospheric neutral density variations observed by GOCE and CHAMP satellites are investigated along with simultaneous measurements of thermospheric nitric oxide density and emissions from TIMED/GUVI and TIMED/SABER, respectively, together with auroral hemispheric power from DMSP, Joule heating rate from the Assimilative Mapping of Ionospheric Electrodynamics (AMIE) algorithm, and solar extreme ultraviolet (EUV) flux from SOHO/SEM and TIMED/SEE. It was found that the apparent storm-time density decrease or “overcooling” is a result of combined effects: an increase in NO cooling and a decrease in both the solar EUV flux and the auroral energy input. It was also found that whether or not a density decrease exists depends on whether there are auroral activities prior to the pre-storm reference day and how the reference day is chosen. The relative contributions of the different drivers (namely, enhanced NO cooling and variations in auroral hemispheric power, Joule heating, and solar EUV at the pre-storm and post-storm days) to thermospheric density decrease remain to be determined.

Space Weather Modeling Capabilities Assessment: Neutral Density for Orbit Determination at low Earth orbit

AbstractThe specification and prediction of density changes in the thermosphere is a key challenge for space weather observations and modeling, because it is one result of complex interactions between the Sun and the terrestrial atmosphere and also because it is of operational importance for tracking objects orbiting in near‐Earth space. For low Earth orbit, neutral density variation is the most important uncertainty for propagation and prediction of orbital elements. A recent international conference conducted under the auspices of the National Aeronautics and Space Administration Community Coordinated Modeling Center included a workshop on neutral density modeling, using both empirical and numerical methods, and resulted in organization of an initial effort in model comparison and evaluation. Here we report on the exploitable density data sets available, the selected years and storm events, and the metrics for complete model assessment. Comparisons between five models (three empirical and two numerical) and neutral density data sets that include measurements by the CHAMP, GRACE, and GOCE satellites are presented as examples of the assessment procedure that will be implemented at Community Coordinated Modeling Center. The models in general performed reasonably well, although seasonal errors sometimes are present, and impulsive geomagnetic storm events remain challenging. Numerical models are still catching up to empirical methods on a statistical basis, but hold great potential for describing these short‐term variations.

Correlations Between the Thermosphere's Semiannual Density Variations and Infrared Emissions Measured With the SABER Instrument.

This paper presents measurements of the amplitudes and timings of the combined, annual, and semiannual variations of thermospheric neutral density, and a comparison of these density variations with measurements of the infrared emissions from carbon dioxide and nitric oxide in the thermosphere. The density values were obtained from measurements of the atmospheric drag experienced by the Challenging Minisatellite Payload, Gravity Recovery and Climate Experiment A, Gravity field and Ocean Circulation Explorer, and three Swarm satellites, while the optical emissions were measured with the Sounding of the Atmosphere using Broadband Emission Radiometry (SABER) instrument on the Thermosphere Ionosphere Mesosphere Energetics and Dynamics satellite. These data span a time period of 16 years. A database containing global average densities that were derived from the orbits of about 5,000 objects (Emmert, 2009, https://doi.org/10.1029/2009JA014102, 2015b, https://doi.org/10.1002/2015JA021047) was employed for calibrating these density data. A comparison with the NRLMSISE‐00 model was used to derive measurements of how much the density changes over time due to these seasonal variations. It is found that the seasonal density oscillations have significant variations in amplitude and timing. In order to test the practicality of using optical emissions as a monitoring tool, the SABER data were fit to the measured variations. Even the most simple fit that used only filtered carbon dioxide emissions had good correlations with the measured oscillations. However, the density oscillations were also well predicted by a simple Fourier series, contrary to original expectations. Nevertheless, measurements of the optical emissions from the thermosphere are expected to have a role in future understanding and prediction of the semiannual variations.

Floating potential of emitting surfaces in plasmas with respect to the space potential

The potential difference between a floating emitting surface and the plasma surrounding it has been described by several sheath models, including the space-charge-limited sheath, the electron sheath with high emission current, and the inverse sheath produced by charge-exchange ion trapping. Our measurements reveal that each of these models has its own regime of validity. We determine the potential of an emissive filament relative to the plasma potential, emphasizing variations in emitted current density and neutral particle density. The potential of a filament in a diffuse plasma is first shown to vanish, consistent with the electron sheath model and increasing electron emission. In a denser plasma with ample neutral pressure, the floating filament potential is positive, as predicted by a derived ion trapping condition. Lastly, the filament floated negatively in a third plasma, where flowing ions and electrons and nonnegligible electric fields may have disrupted ion trapping. Depending on the regime chosen, emitting surfaces can float positively or negatively with respect to the plasma potential.

Design of compact dispersion interferometer with a high efficiency nonlinear crystal and a low power CO2 laser

When the electron density of a plasma generated in high pressure environment is measured by a conventional interferometer, the phase shifts due to changes of the neutral gas density cause significant measurement errors. A dispersion interferometer, which measures the phase shift that arises from dispersion of medium between the fundamental and the second harmonic wavelengths of laser light, can suppress the measured phase shift due to the variations of neutral gas density. In recent years, the CO2 laser dispersion interferometer has been applied to the atmospheric pressure plasmas and its feasibility has been demonstrated. By combining a low power laser and a high efficiency nonlinear crystal for the second harmonic component generation, a compact dispersion interferometer can be designed. The optical design and preliminary experiments are conducted.

Investigation of neutral particle dynamics in Aditya tokamak plasma with DEGAS2 code

Neutral particle behavior in Aditya tokamak, which has a circular poloidal ring limiter at one particular toroidal location, has been investigated using DEGAS2 code. The code is based on the calculation using Monte Carlo algorithms and is mainly used in tokamaks with divertor configuration. This code has been successfully implemented in Aditya tokamak with limiter configuration. The penetration of neutral hydrogen atom is studied with various atomic and molecular contributions and it is found that the maximum contribution comes from the dissociation processes. For the same, Hα spectrum is also simulated and matched with the experimental one. The dominant contribution around 64% comes from molecular dissociation processes and neutral particle is generated by those processes have energy of ~2.0 eV. Furthermore, the variation of neutral hydrogen density and Hα emissivity profile are analysed for various edge temperature profiles and found that there is not much changes in Hα emission at the plasma edge with the variation of edge temperature (7–40 eV).