Neutral Density Research Articles

Investigation of Si incorporation in (010) β-Ga2O3 films grown by plasma-assisted MBE using diluted disilane as Si source and suboxide Ga2O precursor

Traditionally, elemental Ga and Si have been used to supply Ga and Si, respectively, in molecular beam epitaxy (MBE) to grow Si-doped β-Ga2O3. In this work, we investigated the feasibility of enhancing the β-Ga2O3 growth rate by using a Ga-suboxide precursor in a plasma-assisted MBE. Additionally, Si doping of β-Ga2O3 using diluted disilane and Ga-suboxide as the Si and Ga precursors, respectively, was studied. The growth rate and film quality under different suboxide fluxes were inspected. We found that Si concentration has an inverse relationship with Ga2O flux due to atom competition. A room-temperature mobility of 115 cm2/V s was measured for an electron concentration of 1.2 × 1017 cm−3 on the sample grown using a Ga2O beam equivalent pressure of 1.1 × 10−7 Torr and a disilane flow rate of 0.006 sccm. Temperature-dependent Hall characterization was performed on this sample, revealing compensating acceptor and neutral impurity densities of 2.70 × 1015 and 8.23 × 1017 cm−3, respectively.

One-dimensional simulation and validation of divertor detachment induced through nitrogen seeding on HL-2A

Divertor detachment is a promising method to solve the power exhaust problem in tokamak devices or even in future magnetic fusion reactors. In this work, a detachment experiment (HL-2A shot #38008) with mixed gas seeding (60% nitrogen and 40% deuterium) is simulated using the SD1D module in BOUT++. In the process from attachment to detachment, the target electron temperature and the target ion saturation current in simulations are found to be consistent with the experimental results measured by Langmuir probes on the target plate. In order to understand the underlying detachment mechanism on HL-2A, this work analyses the role of different particle species in the cases with different seeding rates. It shows that the plasma density varies little and the density of neutrals ( D and D2 ) slightly decreases with the increase of nitrogen seeding rate, such that plasma–neutral interactions cannot effectively reduce plasma energy and plasma momentum in the divertor. The case with a high seeding rate predicts that increasing seeding rate cannot reach a target temperature lower than 2.5eV , which is the required temperature for strong plasma–molecule interactions. Thus, the plasma–molecule interactions may not be important in the divertor during nitrogen seeding. This work also studies the parallel forces including (1) force due to the parallel electric field, (2) friction forces, (3) ion- and electron-thermal forces, and (4) collisional reactions (e.g. charge exchange recombination and ionisation). It is found that the friction force between nitrogen ions and other particle species (primarily D+ ) is the dominant force towards the target, while ion- and electron-thermal forces are the dominant force pushing nitrogen ions back to upstream. Parallel forces determine the parallel distribution of nitrogen impurities, and therefore decide the region of nitrogen radiation.

Multi-machine benchmark of the self-consistent 1D scrape-off layer model DIV1D from stagnation point to target with SOLPS-ITER

This paper extends a 1D dynamic physics-based model of the scrape-off layer (SOL) plasma, DIV1D, to include the core SOL and possibly a second target. The extended model is benchmarked on 1D mapped SOLPS-ITER simulations to find input settings for DIV1D that allow it to describe SOL plasmas from upstream to target—calibrating it on a scenario and device basis. The benchmark shows a quantitative match between DIV1D and 1D mapped SOLPS-ITER profiles for the heat flux, electron temperature, and electron density within roughly 50% on: (1) the Tokamak Configuration Variable (TCV) for a gas puff scan; (2) a single SOLPS-ITER simulation of the Upgraded Mega Ampere Spherical Tokamak; and (3) the Upgraded Axially Symmetric Divertor EXperiment in Garching Tokamak (AUG) for a simultaneous scan in heating power and gas puff. Once calibrated, DIV1D self-consistently describes dependencies of the SOL solution on core fluxes and external neutral gas densities for a density scan on TCV whereas a varying SOL width is used in DIV1D for AUG to match a simultaneous change in power and density. The ability to calibrate DIV1D on a scenario and device basis is enabled by accounting for cross field transport with an effective flux expansion factor and by allowing neutrals to be exchanged between SOL and adjacent domains.

The Physics of Indirect Estimators of Lyman Continuum Escape and their Application to High-Redshift JWST Galaxies

ABSTRACT Reliable indirect diagnostics of LyC photon escape from galaxies are required to understand which sources were the dominant contributors to reionization. While multiple LyC escape fraction (fesc) indicators have been proposed to trace favourable conditions for LyC leakage from the interstellar medium of low-redshift ‘analogue’ galaxies, it remains unclear whether these are applicable at high redshifts where LyC emission cannot be directly observed. Using a library of 14 120 mock spectra of star-forming galaxies with redshifts 4.64 ≤ z ≤ 10 from the SPHINX20 cosmological radiation hydrodynamics simulation, we develop a framework for the physics that leads to high fesc. We investigate LyC leakage from our galaxies based on the criteria that successful LyC escape diagnostics must (i) track a high-specific star formation rate, (ii) be sensitive to stellar population age in the range 3.5–10 Myr representing the times when supernova first explode to when LyC production significantly drops, and (iii) include a proxy for neutral gas content and gas density in the interstellar medium. O32, ΣSFR, MUV, and H β equivalent width select for one or fewer of our criteria, rendering them either necessary but insufficient or generally poor diagnostics. In contrast, UV slope (β), and E(B − V) match two or more of our criteria, rendering them good fesc diagnostics (albeit with significant scatter). Using our library, we build a quantitative model for predicting fesc based on direct observables. When applied to bright z > 6 Ly α emitters observed with JWST, we find that the majority of them have $f_{\rm esc} \lesssim 10~{{\ \rm per\ cent}}$.

Investigating the role of edge neutrals in exciting tearing mode activity and achieving flat temperature profiles in LTX-β

We present observations, numerical simulations, and analysis from experiments in the Lithium Tokamak Experiment-Beta (LTX-β) in which the electron temperature profile (Te (r)) shifts from flat to peaked and a tearing mode is also destabilized when the average density (ne ave ) exceeds ∼1019 m−3. Flat Te (r) is obtained routinely in LTX-β, with a lithium coated, low-recycling first wall, once the external fueling is stopped and density decays [Boyle et al 2023 Nucl. Fusion 63 056020]. In the present experiment, flat Te profiles can be sustained while maintaining constant ne ave below a line averaged density threshold (ne ave th ) of ∼1019 m−3. Above ne ave th, Te (r) shifts from flat to peaked and a tearing mode is destabilized. Due to low recycling, the achieved ne ave can be controlled precisely by external fueling and hence, a certain threshold of the edge neutral inventory from the external fueling is experimentally manifested through ne ave th . The goal of the present work is to investigate the role of edge neutrals in determining Te (r) and MHD stability in the unique low-recycling regime of LTX-β. Our hypothesis is that the peaking of Te (r) beyond ne ave th is due ultimately to the edge cooling by the cold neutrals beyond a critical fueling flux. At lower fueling flux, flat Te (r) results in broader pressure profile and lower resistivity, which in turn stabilizes the tearing mode. This hypothesis is supported by edge neutral density estimation by DEGAS 2 code. Mode analysis by singular value decomposition confirms the tearing mode structure to be m/n = 2/1 (m and n being the poloidal and toroidal mode numbers). Linear tearing stability analysis with M3D-C1 predicts that plasmas with ne ave > 1019 are highly susceptible to a n = 1 tearing mode. ORBIT simulations, however, confirmed that the tearing modes do not contribute to the loss of fast ions from neutral beam injection. This study shows for the first time that the neutral inventory at the edge could be one of the deciding factors for the achievability of the unique operation regime of flat Te (r) and the excitation of tearing activity that could be disruptive for the plasmas.

Modelling a thrust imparted by a highly ionized magnetic nozzle rf plasma thruster

Influence of the local-ionization-induced neutral depletion on the thrust imparted by the magnetic nozzle plasma thruster is discussed by simply considering reduction of the neutral density due to the ionization in the thruster model combining the global source model and the one-dimensional magnetic nozzle model. When increasing the rf power, it is shown that the increase rate of the plasma density is reduced, while the electron temperature continues to increase due to a decrease in the neutral density. Since the major components of the thrust are originated from the electron pressures in the source and in the magnetic nozzle, the increase in the electron temperature contributes to the increase in the thrust in addition to the gradual density increase by the rf power. The model qualitatively predicts the reduction of the thruster efficiency by the neutral depletion for the high-power condition, compared with the constant neutral density model.

Three-dimensional Venusian ionosphere model

To study the Venusian ionosphere, a 3D ionospheric model is included in the Venus Planetary Climate Model (Venus PCM), which we present here with a number of recent extensions and improvements. Our ionospheric model of Venus consists of 15 charged species (electron and 14 ions), 13 photo-ionization for 8 species, an ion-neutral chemistry (61 reactions) and an ambipolar diffusion scheme for the ion vertical transport. The electron temperature is assumed independent of the solar activity. Simulation results are compared with observation from the Pioneer Venus (PV) and Venus Express (VEX) for high and low/intermediate solar activity respectively. The model shows that ambipolar diffusion dominates photochemical equilibrium above 180 km and above 130 km altitude on the dayside and the nightside respectively. On the dayside, the electron density predicted by Venus PCM and its variation with solar zenith angle (SZA) are in good agreement with PV observations at high solar activity, even if the secondary electron peak is underestimated by the model, due to the absence of (photo-)electron impact ionization process, in Venus PCM. At higher solar activity, the predicted horizontal and vertical variations in electron and ion densities are both significantly underestimated on the nightside, probably due to the currently predicted weak day-to-night transport and the absence of photo-electron impact ionization process which is the main source of ion production on the nightside. On the dayside at low solar activity, the electron density predicted by Venus PCM are over-estimated by a factor 1.25 at the altitude of the main ionospheric peak and by a factor 2–3 at 250 km, compared to VEX observations. We suggest that this difference between high and low solar activity is linked to the overestimation of the neutral density at low solar activity and the independence of electron temperature with solar activity used in Venus PCM, which should have a significant effect on ion chemistry and the scale height of ion species.

Simulating the Effects of Partial Neural Conduction Delays in the Visual Evoked Potential.

The purpose of this study was to understand the double peaks or broadening of P100 observed in some cases of optic neuritis by inducing conduction delays in healthy eyes through stimulus luminance manipulation in analogy to the perceptual Pulfrich effect. Checkerboard pattern reversal visual evoked potentials (VEPs) with check sizes of 0.8degrees, 0.4degrees, and 0.2degrees were recorded in healthy participants using two experiment variants. Variant (1) involved binocular stimulation with inter-ocular luminance difference achieved by a 1.8 neutral density (ND) filter, along with monocular control conditions. Variant (2) included monocular stimulation with hemifields having a luminance difference (half of monitor with ND filter), along with single-hemifield control conditions. In both variants, VEP curves under mixed stimulation were compared to synthesized VEPs computed from offline summation of curves from the relevant control conditions, followed by assessing P100 characteristics. Despite considerable variability between participants, the binocular variant demonstrated marked differences between VEPs from mixed recordings and synthesized curves, whereas in the hemifield variant, agreement was strong. The anticipated double peak or broadened deflection pattern was observed to varying extents in participants, often contingent on check size, with nominal peak time frequently failing to indicate partial conduction delays. The present findings corroborate the hypothesis that nominal peak time does not always reflect conduction delays if only a subset of fiber bundles is affected. Peak shape might provide additional diagnostic evidence of a partial conduction delay. Enhancing the understanding of VEP waveform changes associated with partial conduction delays could offer diagnostic insights for optic neuritis.

The role of isotope mass on neutral fueling and density pedestal structure in the DIII-D tokamak

Experimental measurements on DIII-D of hydrogen neutral penetration lengths ( λn0 ) on the high field side (HFS) are longer by a factor of 2 than for deuterium consistent with the thermal velocity ratio for neutrals at the same temperature (vthH/vthD=2) . This ratio is constant for both low and high pedestal electron density. At low pedestal density (ne∼4×1019 m −3) , the neutral penetration length is greater than the density pedestal width for both isotopes, and the additional 40% increase of neutral penetration in hydrogen widens the pedestal by the same amount. As the density pedestal height increases (ne∼6×1019 m −3) , the neutral penetration lengths drop below the density pedestal widths for both isotopes, and the increased penetration of hydrogen has no increased effect on the pedestal width compared to deuterium. Extrapolating to future reactor-relevant high electron density pedestals, the isotope-mass change in neutral fueling on the HFS from the deepest neutral penetration of hydrogen, to the shortest neutral penetration of tritium will be negligible (0.2–0.4 cm) in comparison to estimates of the density pedestal width (6–8.5 cm).

Analysis of anode surface roughness influence on heaterless hollow cathode discharge

This research delves into the influence of cylindrical and planar anode surface roughness on heaterless hollow cathode discharge characteristics. Three surface roughness levels, Ra 1.6 μm, Ra 3.2 μm, and Ra 6.4 μm, along the cylindrical anode’s azimuthal direction and the planar anode’s radial direction, have been selected and modified by the machining process. A central finding is the correlation between the discharge voltage and anode surface roughness. As the roughness increases, cylindrical and planar anodes require less discharge voltage to sustain the primary discharge. This reduction is likely due to the enhanced surface area from the roughness, which aids in electron current collection. The discharge voltage oscillations appear inversely related to the surface roughness of planar anodes, which may be associated with changes in the neutral gas density gradient between the planar anode and the cathode. In contrast, cylindrical anodes show fewer effects from their surface roughness, likely due to the distinct neutral gas flow dynamics. These findings offer insights into standardizing hollow cathode testing and allow future research to explore these interactions more deeply.

On the role of minor neutrals in determining the characteristic features of the Venus ionosphere at low altitudes

ABSTRACT Electron density profiles obtained from radio occultation experiments reveal significant ionization occurring at altitudes below 120 km in the Venusian ionosphere. The presence of such ionization remains unexplained when considering only major neutrals and the resulting ion chemistry. In this study, we investigate the role of minor neutrals, including NO, O2, C, Ar, H2, and H, in the lower altitudes of the Venusian ionosphere using a one-dimensional photochemical model and Akatsuki radio science measurements. We derive the density of minor neutrals from the profiles of major neutrals. Our analysis demonstrates that NO+ and O$_2^+$ ions exhibit a notable presence below 140 km. When we incorporate the minor neutrals, electron density increases below 115 km but remains unchanged above this altitude. Specifically, the density of NO+ increases below 125 km, with corresponding reductions in O$_2^+$. The increase in NO+ above 115 km can be attributed to the charge-exchange reaction between O$_2^+$ and NO. However, below 115 km, the photoionization of NO results in a net increase in NO+ and electron density. These findings highlight the significant role played by minor neutrals, particularly below 115 km, in shaping the Venusian ionosphere at low altitudes.

Ultra‐Fast Kelvin Wave Packets in Mars' Atmosphere and Their Interactions With Tides as Viewed by MAVEN/NGIMS and MRO/MCS

AbstractA key element of successful aerobraking operations at Mars is accurate thermospheric density predictions. Evidence suggests that much of the longitude variability in Mars' aerobraking region is associated with atmospheric tides, and the day‐to‐day variability is connected with tidal modulation by longer‐period global‐scale waves. Specifically, ultra‐fast Kelvin waves (UFKWs) and their modulation of the tidal spectrum play a key role in coupling Mars' lower (<∼80 km) and middle (∼80–100 km) atmosphere with the aerobraking region above. In this study, over 5 years of Mars Atmosphere and Volatile Evolution Neutral Gas and Ion Mass Spectrometer CO2 density observations are employed to reveal prominent, frequent, and persistent 2.5‐ to 4.5‐day UFKW packets in the whole Martian middle and upper thermosphere (ca. 150–200 km), and large secondary waves arising from their nonlinear interactions with the tidal spectrum. Detailed analyses focusing on a prominent ∼2.5‐day UFKW event in late 2015 demonstrate primary and secondary wave amplitudes growing twofold with altitude from ∼7%–14% near 150 km to ∼12%–25% near 200 km and their combined effects to account for ∼60%–80% of the altitude‐longitudinal variability of Mars' thermospheric density. Concurrent temperature measurements from Mars Reconnaissance Orbiter Mars Climate Sounder reveal consistent wave signatures near 80 km altitude suggesting propagation of both primary and secondary waves from the lower atmosphere. This study demonstrates that UFKWs and secondary waves from UFKW‐tide interactions are sources of significant altitude‐longitude variability in the Mars' aerobraking region that should be accounted for when analyzing satellite observations and nonlinear models.

Low-pressure inductively coupled plasmas in hydrogen: impact of gas heating on the spatial distribution of atomic hydrogen and vibrationally excited states

Non-equilibrium inductively coupled plasmas (ICPs) operating in hydrogen are of significant interest for applications including large-area materials processing. Increasing control of spatial gas heating, which drives the formation of neutral species density gradients and the rate of gas-temperature-dependent reactions, is critical. In this study, we use 2D fluid-kinetic simulations with the Hybrid Plasma Equipment Model to investigate the spatially resolved production of atomic hydrogen in a low-pressure planar ICP operating in pure hydrogen (10–20 Pa or 0.075–0.15 Torr, 300 W). The reaction set incorporates self-consistent calculation of the spatially resolved gas temperature and 14 vibrationally excited states. We find that the formation of neutral-gas density gradients, which result from spatially non-uniform electrical power deposition at constant pressure, can drive significant variations in the vibrational distribution function and density of atomic hydrogen when gas heating is spatially resolved. This highlights the significance of spatial gas heating on the production of reactive species in relatively high-power-density plasma processing sources.

Retrieval of Thermospheric O and N2 Densities From ICON EUV

AbstractAs activity in Earth orbit continues to grow, it is important to characterize the environment of near‐Earth space. One means of remotely sensing lower thermospheric neutrals is by measurement of O and N2 density through the observation of far‐ultraviolet (FUV) airglow of atomic oxygen at 135.6 nm and the N2 Lyman‐Birge‐Hopfield (LBH) bands (~130–180 nm), as has been done on the Ionospheric Connection Explorer (ICON), Global‐scale Observations of the Limb and Disk (GOLD), and Thermosphere Ionosphere Mesosphere Energetics and Dynamics (TIMED) missions. This technique is not without limitations, however, as the FUV measurements suffer from contamination by ionospheric emissions at low latitudes and auroral emissions excited by precipitating energetic electrons and protons at high latitudes. Previous work has shown the potential for making measurements of O and N2 density in the lower‐middle thermosphere using observations of extreme‐ultraviolet (EUV) airglow. This measurement approach has a potential advantage in that it does not have an inherent ionospheric emission that must be accounted for. Additionally, these emissions are primarily excited directly by solar UV rather than electron impact and thus have the potential to enable expansion of neutral density observations into the auroral zone and polar cap where the FUV measurement cannot be applied. This article demonstrates a new approach and algorithm designed to retrieve thermospheric O and N2 density from 150 to 400 km using measurements from the ICON EUV instrument. The retrieval results throughout 2020 are summarized and compared to measurements from ICON FUV, GOLD, and SWARM.

Including the parallel mass flow in calculating the steady-state solutions and stability of the momentum balance equations for a quasisymmetric stellarator

The Helically Symmetric Experiment (HSX) is a quasisymmetric stellarator with minimal parallel viscous damping in a helical direction. The parallel flow (Vǁ) along the magnetic field is similarly weakly damped by viscosity. In this paper, the self-consistent steady-state parallel and poloidal momentum balance equations are used to show that a large Vǁ on the order of the ion thermal velocity can increase the ion resonant radial electric field (Er) beyond the value calculated using the typical approximation that Vǁ is zero. By altering the damping of Vǁ, either by degrading the quasisymmetry or varying the neutral density, the ion resonant Er can shift in a controllable fashion. It is shown explicitly that there exist stable and unstable steady-state solutions in the two-dimensional space of Vǁ and Er. A stability analysis of each solution is performed by calculating the eigenvalues and eigenvectors of the Jacobian. The unstable solution corresponds to a saddle point in which the eigenvalues have opposite signs. The analysis leads to the conclusion that unstable solutions occur when the derivative of the total poloidal damping with respect to Er is positive. A hysteresis in Er and Vǁ is observed when the radial current density is linearly increased to a maximum and then decreased back to zero. Jumps in the radial electric field and the parallel flow are observed as the radial current density drives the evolution from one stable point to the next. This result is similar to experimental data observed on several devices.

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.

Lyα Absorption in a “Croissant-like” Heliosphere

Lyα absorption profiles have been used to detect astrospheres and heliospheric absorption from the hydrogen wall and heliotail. Using magnetohydrodynamic models of the heliosphere, we can compare simulated to observed Lyα profiles to probe the neutral hydrogen within and near the heliosphere. There is an ongoing controversy whether the heliosphere has a long “comet-like” tail or a short “croissant-like” tail. Here we present the first Lyα absorption investigation using a croissant-like heliosphere. With identical boundary conditions we compare the BU model, which presents a croissant-like tail, and the Moscow model, which presents a comet-like tail. The BU and Moscow models present nearly identical Lyα profiles toward nose targets (α Cen and 36 Oph). Differences in Lyα profiles are shown toward the tail target (HD 35296). Despite the shortened heliotail of the croissant model, significant downwind heliosheath absorption is seen, just 5% shallower and shifted by 4 km s−1. This implies that an extended tail model is not required to reproduce the heliosheath Lyα absorption observations. Finer observation gratings may be able to resolve this shift. Additionally, when using higher interstellar medium (ISM) neutral and plasma densities and lower magnetic field (∣B LISM∣ = 3.2 μG, α BV ≈ 40°) than in the Moscow model, we find better agreement with observed Lyα profiles. None of the models presented show agreement in all directions simultaneously. Furthermore, we show that for the ISM conditions with the least certainty (n p,LISM, n H,LISM, T LISM, B LISM), B LISM has the most significant effect on the structure of the hydrogen wall and Lyα profiles.

Magnetic Signatures of the Interaction Between Europa and Jupiter's Magnetosphere During the Juno Flyby

AbstractBased on a hybrid model of Europa's magnetospheric interaction, we provide context for the magnetic field perturbations observed by the Juno spacecraft during its only close flyby of the moon in September 2022. By systematically varying the incident flow conditions and the density profile of Europa's atmosphere, we demonstrate that the observed, large‐scale signatures of magnetic field draping are consistent with a dawn‐dusk asymmetry in the moon's neutral envelope. During the flyby, such an asymmetry would have enhanced the magnetic perturbations in Europa's anti‐Jovian hemisphere, explaining why the spacecraft already detected strong field line draping while still several moon radii away. Conversely, a reduced neutral density in the sub‐Jovian hemisphere can explain why the perturbations in the flow‐aligned field component remained nearly constant as Juno approached Europa. While a dawn‐dusk asymmetry in Europa's atmosphere has been predicted by theoretical work, our results provide the first in situ hints of its presence.

Evaluation of TIEGCM based on GOCE neutral density

The Thermosphere Ionosphere Electrodynamic General Circulation Model (TIEGCM), as one of the most advanced physical models of the Earth’s thermosphere and ionosphere, is not only widely used in scientific research, but also has essential reference value in aerospace operations. In this study, we use Gravity field and steady-state Ocean Circulation Explorer (GOCE) neutral density to evaluate the accuracy of the TIEGCM. The assessment is performed on both time and spatial scales. The time scales are conducted annually, monthly, and daily, while the spatial scales are carried out in terms of altitude, latitude, and local time. On the time scales, the performance of the TIEGCM on the monthly time scale is better than that on the annual time scale. Also, the performance on the daily time scale is better than that on the monthly time scale. The relative deviation shows a significant seasonal variation, that is, larger in winter and summer and smaller in spring and autumn. In addition, the relative deviation shows a negative correlation with F10.7 and Ap. On the spatial scale, with the increase in altitude, the average relative deviation of the model becomes larger in general. The relative deviation is usually larger at middle latitudes in the Northern Hemisphere and high latitudes in the Southern Hemisphere. Finally, on the scale of local time, the relative deviation changes more dramatically in local morning than at dusk.

Electrohydrodynamically printed microlens arrays with controllable curvature based on surface functionalization

Precisely controlling the optical characteristics of the microlens array (MLA) is a key issue for the MLA's fabrication and its applications. A cost-effective method for fabricating MLA with different curvatures and diameters by surface functionalization and electrohydrodynamic inkjet (E-jet) printing is proposed. The surface with self-assembled monolayers has low surface energy, which greatly improves the contact angle of the microlens. The UV-ozone treatment and neutral density filters with different optical densities are used to control the wettability of the substrates. Based on this method, the curvatures of the microlens can be controlled easily and the numerical aperture (NA) of the microlens can be adjusted from 0.18 to 0.49. Compared with traditional fabrication techniques, the method we demonstrated is flexible, low-cost, rapid, and capable of fabricating MLA with large area, high packing density, desirable NA, and diameters.