Small Scale Research Articles

Density Fluctuations in the Intracluster Medium: An Attempt to Constrain Viscosity with Cosmological Simulations

The impact of viscosity in the intracluster medium (ICM) is still an open question in astrophysics. To address this problem, we have run a set of cosmological simulations of three galaxy clusters with a mass larger than M Vir > 1015 M ⊙ at z = 0 using the smoothed particle magnetohydrodynamics-code OpenGadget3. We aim to quantify the influence of viscosity and constrain its value in the ICM. Our results show significant morphological differences at small scales, temperature variations, and density fluctuations induced by viscosity. We observe a suppression of instabilities at small scales, resulting in a more filamentary structure and a larger amount of small structures due to the lack of mixing with the medium. The conversion of kinetic to internal energy leads to an increase of the virial temperature of the cluster of ∼5%–10%, while the denser regions remain cold. The amplitude of density and velocity fluctuations are found to increase with viscosity. However, comparison with observational data indicates that the simulations, regardless of the viscosity, match the observed slope of the amplitude of density fluctuations, challenging the direct constraint of viscosity solely through density fluctuations. Furthermore, the ratio of density to velocity fluctuations remains close to 1 regardless of the amount of viscosity, in agreement with the theoretical expectations. Our results show for the first time in a cosmological simulation of a galaxy cluster the effect of viscosity in the ICM, a study that is currently missing in the literature.

Meso-inflationary Peccei–Quinn Symmetry Breaking with Nonminimal Coupling

We study a realization of the inflationary scenario where the Peccei–Quinn (PQ) symmetry is spontaneously broken during inflation, facilitated by its nonminimal coupling to gravity. This results in, effectively, a two-field inflation: the early stage is driven by an inflaton field with the PQ symmetry intact, and the later stage is driven by the PQ scalar after its effective mass becomes tachyonic, causing destabilization from the origin. The nonminimal coupling serves the dual purpose of restoring the PQ symmetry during early inflation and flattening the PQ potential posttachyonic shift, allowing for continued slow-roll. We analyze the inflationary background solutions and scalar perturbations, which are amplified at small scales via significant isocurvature perturbations generated near the symmetry-breaking epoch. These perturbations lead to second-order gravitational waves, detectable by next-generation space-based experiments.

Microlevel analysis of how itinerant immigrant entrepreneurs underwrite local economic development in an African city

ABSTRACT Itinerant merchants are actively involved in various economic activities within and around African cities, albeit on a small scale. However, there remains a gap in understanding their contribution to host cities’ local economic development (LED). This research aims to fill this gap by examining how itinerant immigrant retailers support LED in host cities, using Wa in Ghana. Recognising their role is essential for informing location-specific development plans to boost local economies. Employing a qualitative research approach, we interviewed 52 participants, including itinerant immigrants, city residents, and officials from relevant state agencies. Thematic analysis using NVivo version 11 was employed to analyse the responses. Itinerant merchants play a crucial role in local economies by creating employment opportunities, providing income sources for residents, ensuring convenient access to goods and services, and contributing to market diversification. We discuss how these insights could inform LED policies and planning in sub-Saharan Africa.

Scale-up of Sodium Persulfate Mediated, Nitroxide Catalyzed Oxidative Functionalization Reactions.

Oxidation is a valuable tool in preparative organic chemistry. Oxoammonium salts and nitroxides have proven valuable as reagents and catalysts in this endeavor. The objective of this study is to scale up the oxidative amidation, ester formation, and nitrile formation using nitroxide as an organocatalyst. Oxidative functionalization reactions were scaled from the 1 mmol to the 1 mole level. Sodium persulfate was used as the primary oxidant, and a nitroxide was employed as a catalyst. The products of the reactions were isolated in analytically pure form by extraction with no need for column chromatography. The oxidative amidation and esterification of aldehydes can be scaled up from 1 mmol to 1 mole effectively, with comparable product yields being obtained at each increment. This work shows that conditions developed on a small scale can be transferred to a larger scale without reoptimization. The oxidative functionalization of aldehydes to prepare nitriles is not amenable to direct scale-up due to the concomitant formation of significant quantities of the corresponding carboxylic acid, thereby compromising the product yield. Two of the three oxidative transformations studied here can be scaled up successfully from the 1 mmol to the 1 mole level.

Slow-fast dynamics in non-linear enzyme cascades gives rise to spatial multiscaling.

Slow-fast dynamics in a well-mixed biological system is known to be associated with functional modularity with temporal hierarchy of processes. Yet, the significance of slow-fast dynamics in spatially distributed reaction-diffusion systems remains to be clarified.We assessed such systems with slow-fast dynamics, such as blood coagulation and predator-prey models to identify spatial patterns with independent parametric regulation, and performed their reaction-rate based sensitivity analysis with the help of computational one-dimensional reaction-diffusion models.In blood coagulation, spatial distribution of fibrin had two scales: a larger region with a high and almost constant fibrin concentration (corresponding to solid blood clot) and a narrow transition zone, where fibrin concentration decayed to zero (semi-liquid clot). Slow variables (e.g., the concentration of the coagulation factor IXa) regulated the clot size (large scale), while fast variables (e.g., concentrations of thrombin and factor Xa) regulated the thickness of the gel-to-liquid transition zone (small scale). Analysis of the simplified model of blood coagulation, composed only of factor X and IX activation, and of the predator-prey model confirmed this finding to be sufficiently general, and revealed that the limitation of fast enzyme production (due to precursor depletion or to kinetics with saturation) was crucial for this multi-scaling. Independent regulation of spatial scales also required species with slow kinetics to be upstream from species with fast kinetics.The shared slow-fast dynamics of the examined systems endow them with the ability to form spatial patterns, representing a new mechanism complementing the Turing-type and the Positional Information type pattern formation.

Small scale laboratory monotonic and cyclic pull out testing on grout and resin encapsulated cable bolts

Axial studies on cable bolts can be conducted using various scale testing apparatuses. Large scale testing, while providing a powerful platform for testing, is expensive and time consuming. This study presents details of a small scale pull out testing campaign on cable bolts and investigates the results achieved. Six popular types of cable bolts were studied using an anti rotation apparatus while encapsulated in cementitious grout and resin. The resin samples were tested under both monotonic and cyclic loading patterns. The results showed that grouted bulbed cables require higher displacement to reach their maximum load capacity which is lost at failure, while plain cables tend to hold lower loads for a longer time. Resin samples provided strain softening behaviour with low capacities, particularly in absence of cable indentation or bulbs. Cyclic loading tended to adversely affect the post peak behaviour of the resin samples, especially in the bulbed cables. Failed samples inspected after the testing suggested a non-uniform damage profile along the cable with extensive damage at the exit point transitioning into almost no damage at the entry point.

Fluid wetting and penetration characteristics in T-shaped microchannels

A thorough understanding of media tightness in automotive electronics is crucial for ensuring more reliable and compact product designs, ultimately improving product quality. Concerning the fundamental characteristics of fluid leakage issues, the dynamic wetting and penetration behavior on small scales is of special interest and importance. In this work, four T-shaped microchannels with one inlet and two outlets are experimentally investigated in terms of contact angle dynamics and interface movement over time, generating novel insight into the wetting mechanisms and fluid distribution. With a main channel width of 1 mm, a crevice width of w=0.3mm,0.4mm\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$w = {0.3}\\,\\hbox {mm}, {0.4}\\,\\hbox {mm}$$\\end{document} and a rounding edge radius of r=0.1mm,0.2mm\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$r = {0.1}\\,\\hbox {mm}, {0.2}\\,\\hbox {mm}$$\\end{document}, the geometrical effects on the fluid penetration depth in the crevice and the interface edge pinning effect are analyzed quantitatively using an automated image processing procedure. It is found that the measured dynamic contact angles in all parts can be well described by molecular kinetic theory using local contact line velocities, even with local surface effects and abrupt geometry changes. Moreover, a smaller crevice width, a sharper edge and a larger flow velocity tend to enhance the interface pinning effect and prevent fluid penetration into the crevice. The rounding radius has a more significant effect on the interface pinning compared with crevice width. The experimental data and image processing algorithm are made publicly available.

Investigating the Kinematics of Central and Satellite Galaxies Using Normalizing Flows

Galaxy clustering contains information on cosmology, galaxy evolution, and the relationship between galaxies and their dark matter hosts. On small scales, the detailed kinematics of galaxies within their host halos determines the galaxy clustering. In this paper, we investigate the dependence of the central and satellite galaxy kinematics on θ , the intrinsic host halo properties (mass, spin, concentration), cosmology (Ωm, σ 8), and baryonic feedback from active galactic nuclei and supernovae (A AGN1, A AGN2, A SN1, A SN2). We utilize 2000 hydrodynamic simulations in CAMELS run using IllustrisTNG and SIMBA galaxy formation models. Focusing on central and satellite galaxies with M * > 109 M ⊙, we apply neural density estimation (NDE) with normalizing flows to estimate their p(Δr ∣ θ ) and p(Δv ∣ θ ), where Δr and Δv are the magnitudes of the halocentric spatial and velocity offsets. With NDE, we accurately capture the dependence of galaxy kinematics on each component of θ . For central galaxies, we identify significant spatial and velocity biases dependent on halo mass, concentration, and spin. For satellite distributions, we find significant deviations from a Navarro–Frenk–White profile and evidence that they consist of distinct orbiting and infalling populations. However, we find no strong dependence on θ besides a weak dependence on host halo spin. For both central and satellite galaxies, there is no notable dependence on cosmological parameters and baryonic feedback. These results provide key insights for improving the current halo occupation distribution (HOD) models. This work is the first in a series that will reexamine and develop HOD frameworks for improved modeling of galaxy clustering at smaller scales.

New insights into impact-induced removal of the deposited droplet

This paper presents a comprehensive investigation into the collision dynamics of equal and unequal-sized nanodroplets on a flat surface using molecular dynamics simulations, revealing new insights into scaling laws and energy dissipation mechanisms. The simulations, conducted with the Large-Scale Atomic/Molecular Massively Parallel Simulator software, involved an initially stationary droplet on the surface and a suspended droplet with varying diameter ratios (λ) and impact velocities. The results show that at low Weber numbers (We < 24.15), the droplets tend to deposit after impact, while at higher Weber numbers (We ≥ 24.15), they undergo spreading and retraction, ultimately rebounding. The study reveals that the dimensionless contact time (t*) and maximum spreading factor (βmax*) in collisions between droplets of different sizes do not follow the same scaling relationship observed in single nanodroplet impacts. By redefining the Weber and Reynolds numbers (Re), the new scaling relationships t* ∼ We2/3Re−1/3λ−1/3 and βmax* ∼ We2/3Re−1/3λ−1/3 are proposed and validated. This work represents a further in-depth study of previous research on single nanodroplet impact, introducing for the first time the diameter ratio in unequal droplet impacts into the variation patterns of contact time and maximum spreading diameter. Moreover, these findings highlight the importance of revisiting and potentially revising classical theories to accommodate the unique physical phenomena that emerge at smaller scales.

Assessing Sustainable Urbanization in Small Cities: A Life Satisfaction Perspective Using the Delphi Method and Analytic Hierarchy Process (AHP)

In recent years, small cities have experienced rapid development, and the living standards of their residents have increasingly approximated those of larger urban centers. Due to their smaller scale, these cities have inherent advantages in transitioning to sustainable development models. The post-pandemic era, characterized by job reductions in large cities and the economic appeal of lower living costs, has seen a migration flow from large to small cities. The Chinese government has traditionally emphasized urban growth. While recent policies have started to incorporate environmental considerations, there remains a limited focus on the well-being and living conditions of city inhabitants. This study aims to evaluate the sustainable urbanization of small cities from a life satisfaction perspective, emphasizing indicators closely related to residents' well-being. Employing the Delphi method and Analytic Hierarchy Process (AHP), this research identifies persuasive evaluation metrics tailored to the unique contexts of small cities. This work addresses the research gap in sustainable development indicators at the small city level. By emphasizing the perspective of residents' quality of life, it advocates for local governments to prioritize the well-being and experiences of their inhabitants.

Ectoparasite infection levels differ between fish from upwelling-exposed and sheltered rocky reefs areas in Brazil

Ubiquitous and abundant within marine ecosystems, parasites play essential ecological roles such as shaping host population dynamics, altering competition between species, and influencing energy flows through communities. Their diversity and population dynamics are demonstrably shaped by both seasonal and geographic variations. These variations have been often explored at broad spatial scale. However, parasite communities can exhibit significant disparities even at small spatial scales, driven by factors such as wave exposure, temperature fluctuations, and benthic habitat composition. We investigated how crustacean parasites of fish – caligids and gnathiids - differed between two distinct habitats, which are separated solely by a few kilometres, at Arraial do Cabo, Brazil. These two habitats were characterised by a sheltered embayment (hereinafter referred to as “inside”) or an exposed upwelling habitat (hereinafter referred to as “outside”). Individual fish from four species were examined in both habitats. We found that the infestation rate of caligids varied among fish species and, gnathiids varied between the two sampling sites. Gnathiids were absent from fish outside, while they were present on fish inside the embayment. This disparity suggests a critical role of local environmental factors in shaping gnathiid distribution. Potential drivers include temperature fluctuations, substrate composition, and wave exposure, which differed markedly between the two sites. Conversely, caligid parasites infected fish in both locations. While environmental factors may also influence caligid abundance, they appear to exhibit greater tolerance compared to gnathiids. These findings emphasize the importance of incorporating fine-scale environmental heterogeneity when investigating parasite distribution patterns.

She’s Got Her Mother’s Hair: Unveiling the Origin of Black Hole Magnetic Fields through Stellar to Collapsar Simulations

Abstract Relativistic jets from a Kerr black hole (BH) following the core collapse of a massive star (“collapsar”) is a leading model for gamma-ray bursts (GRBs). However, the two key ingredients for a Blandford–Znajek-powered jet—rapid rotation and a strong magnetic field—seem mutually exclusive. Strong fields in the progenitor star’s core transport angular momentum outward more quickly, slowing down the core before collapse. Through innovative multidisciplinary modeling, we first use MESA stellar evolution models followed to core collapse to explicitly show that the small length scale of the instabilities—likely responsible for angular momentum transport in the core (e.g., Tayler–Spruit)—results in a low net magnetic flux fed to the BH horizon, far too small to power GRB jets. Instead, we propose a novel scenario in which collapsar BHs acquire their magnetic “hair” from their progenitor proto–neutron star (PNS), which is likely highly magnetized from an internal dynamo. We evaluate the conditions for the BH accretion disk to pin the PNS magnetosphere to its horizon immediately after the collapse. Our results show that the PNS spin-down energy released before collapse matches the kinetic energy of Type Ic-BL supernovae, while the nascent BH’s spin and magnetic flux produce jets consistent with observed GRB characteristics. We map our MESA models to 3D general-relativistic magnetohydrodynamic simulations and confirm that accretion disks confine the strong magnetic flux initiated near a rotating BH, enabling the launch of successful GRB jets, whereas a slower-spinning BH or one without a disk fails to do so.

Diagenetic History and Biosignature Preservation Potential of Fine‐Grained Rocks at Hogwallow Flats, Jezero Crater, Mars

AbstractThe Mars 2020 Perseverance rover discovered fine‐grained clastic sedimentary rocks in the “Hogwallow Flats” member of the “Shenandoah” formation at Jezero crater, Mars. The Hogwallow Flats member shows evidence of multiple phases of diagenesis including Fe/Mg‐sulfate‐rich (20–30 wt. %) outcrop transitioning downward into red‐purple‐gray mottled outcrop, Fe/Mg clay minerals and oxides, putative concretions, occasional Ca sulfate‐filled fractures, and variable redox state over small (cm) spatial scales. This work uses Mastcam‐Z and SuperCam instrument data to characterize and interpret the sedimentary facies, mineralogy and diagenetic features of the Hogwallow Flats member. The lateral continuity of bedrock similar in tone and morphology to Hogwallow Flats that occurs over several km within the western Jezero sedimentary fan suggests widespread deposition in a lacustrine or alluvial floodplain setting. Following deposition, sediments interacted with multiple fluids of variable redox state and salinity under habitable conditions. Three drilled sample cores were collected from this interval of the Shenandoah formation as part of the Mars Sample Return campaign. These samples have very high potential to preserve organic compounds and biosignatures. Drill cores may partially include dark‐toned mottled outcrop that lies directly below light‐toned, sulfate‐cemented outcrop. This facies may represent some of the least oxidized material observed at this interval of the Shenandoah formation. This work reconstructs the diagenetic history of the Hogwallow Flats member and discusses implications for biosignature preservation in rock samples for possible return to Earth.

Higher order method based on the coupling of local discontinuous Galerkin method with multistep implicit-explicit time-marching for the micropolar Navier–Stokes equations

In this paper, the spatial local discontinuous Galerkin (LDG) approximation coupled with the temporal multistep implicit-explicit (IMEX) evolution for the micropolar Navier–Stokes equations (MNSE) is adopted to construct an efficient fully discrete method. First, the multistep IMEX-LDG methods are constructed up to the third order, which have small computational scale and facilitate implementation to higher orders. Second, the unconditional stability of the fully discrete method on Cartesian grids is derived theoretically, meaning that the temporal step size τ is upper bounded by a positive constant independent of the spatial mesh size h. Last, numerical results for nonlinear MNSE with different boundary conditions are presented, confirming the theoretical validity and effectiveness of the method.

Insights Into Magma Reservoir Dynamics From a Global Comparison of Volcanic and Plutonic Zircon Trace Element Variability in Individual Hand Samples

AbstractTrace element compositional trends in zircons separated from single hand samples have been used to infer dynamic processes in magma reservoirs. Here, we compile published zircon trace element chemistry to quantify any systematic difference between the range of compositions observed in zircon from individual volcanic and plutonic hand samples and compare these results with geochemical modeling to derive implications for magma reservoir dynamics. We find that both rock types span a wide range of hand‐sample scale variability (i.e., wide range of coefficients of variation), but there is no systematic difference in the average variability between plutonic and volcanic samples (i.e., no difference in the mean coefficient of variation). This indicates that dynamic processes related to eruption are not necessarily required as a fundamental process to create hand sample‐scale compositional heterogeneity beyond what is present due to dynamic processes in the reservoir recorded in plutonic samples. Modeling of felsic systems (>68.5 wt.% SiO2) indicates that the similar average variability in felsic volcanic and plutonic hand samples cannot be reproduced by closed‐system crystallization of compositionally distinct melts locally within a magma reservoir (i.e., isolated melt pockets in a crystal mush) but requires mixing of at least two felsic melt compositions at a small spatial scale. This study provides a framework for focused studies on individual volcanic‐plutonic systems exploring how plutonic and volcanic zircon compositional variability records the time and length scales of magma reservoir processes.

Pangenome and pantranscriptome as the new reference for gene family characterisation – a case study of basic helix-loop-helix (bHLH) genes in barley

Genome-wide identification and comparative gene family analyses have been commonly performed to investigate species-specific evolution linked to various traits and molecular pathways. However, most previous studies were limited to gene screening in a single reference genome, failing to account for the gene presence/absence variations (gPAVs) in a species. Here, we propose an innovative pangenome-based approach of gene family analyses based on orthologous gene groups (OGGs). Using the basic helix-loop-helix (bHLH) transcription factor family in barley as an example, we identified 161 ∼ 176 bHLHs in 20 barley genomes, which could be classified into 201 OGGs. These 201 OGGs were further classified into 140 core, 12 soft-core, 29 shell, and 20 line-specific/cloud bHLHs, revealing a complete profile of bHLH in barley. Using a genome-scan approach, we overcome the genome annotation bias and identified on average 1.5 un-annotated core bHLHs per barley genome. We found that all core bHLHs belong to whole genome/segmental duplicates whilst dispensable bHLHs were more likely to result from small scale duplication events. Interestingly, we noticed that the dispensable bHLHs tended to enrich in specific subfamilies SF13, SF27, and SF28, implying the potential biased expansion of specific bHLHs in barley. We found that 50% of the bHLHs contain at least one intact transposon element within the 2kb upstream-to-downstream region. bHLHs with CNV have 1.48 TEs on average, significantly higher than 1.36 for core bHLH without CNV, supporting TEs’ potential role in bHLH expansion. Selection pressure analyses showed that dispensable bHLHs had experienced clear relaxed selection compared to core bHLHs, consistent with their conservation patterns. We further integrate pangenome with recently available barley pantranscriptome data in 5 tissues and discovered apparent transcriptional divergence within and across bHLH subfamilies. We conclude that pangenome-based gene family analyses can better describe the genuine evolution status of bHLHs untapped before and provided novel insights into bHLH evolution in barley. We expect this study will inspire similar analyses in many other gene families and species.

Mechanical properties of laser powder bed fusion processed Inconel alloy IN718 in different heat treatment conditions through small scale specimen testing

Mechanical properties of laser powder bed fusion processed Inconel alloy IN718 in different heat treatment conditions through small scale specimen testing

HBIM for conservation and valorization of structural heritage: The Stylite Tower at Umm ar-Rasas, Jordan

The Heritage Building Information Modeling (HBIM) is a consistent approach to support the multidisciplinary process of knowledge, conservation, and valorization of the historical building heritage, exploiting the combination of advanced 3D modelling and heterogeneous informative data. In the present paper, a highly customized HBIM workflow is developed for the enhancement and conservation of the byzantine Stylite Tower in the Jordan archaeological site of Umm ar-Rasas. The case study is considered a benchmark to test the effectiveness and the scalability of the proposed methodology for the broader set of ancient structures located in archaeological sites. The proposed HBIM process has a dual aim: the organization of heterogeneous and fragmented sources concerning the history and the actual state maintenance of the Tower and the development of a flexible tool to perform fast qualitative and quantitative analysis oriented towards conservation and management projects. The small scale of the structure and an innovative semi-automatic modelling process based on visual coding allowed for a stone-by-stone segmentation of the Tower. The model procedure is based on an information management workflow based on CIDOC-CRM ontology and the interchange format Industry Foundation Classes to ensure interoperability, developing an open-access interactive visualization of the model of the Tower, embedding organized document sources. The proposed HBIM model can also be a valuable tool for enhanced monument fruition and dissemination activities of the structure.

Radio–FIR correlation: A probe into cosmic ray propagation in the nearby galaxy IC 342

Resolved studies of the correlation between the radio and far-infrared (FIR) emission from galaxies at different frequencies can unveil the interplay between star formation and the relativistic interstellar medium (ISM). Thanks to the LOFAR LoTSS observations combined with VLA, Herschel, and WISE data, we study the role of cosmic rays and magnetic fields in the radio–FIR correlation on scales of ≳200 pc in the nearby galaxy IC 342. The thermal emission traced by the 22 μm emission, constitutes about 6%, 13%, and 30% of the observed radio emission at 0.14, 1.4, and 4.8 GHz, respectively, in star-forming regions and less in other parts. The nonthermal spectral index becomes flatter at frequencies lower than 1.4 GHz (αn = −0.51 ± 0.09, Sν ∝ ναn) than between 1.4 and 4.8 GHz (αn = −1.06 ± 0.19) on average, and this flattening occurs not only in star-forming regions but also in the diffuse ISM. The radio–FIR correlation holds at all radio frequencies; however, it is tighter at higher radio frequencies. A multi-scale analysis shows that this correlation cannot be maintained on small scales due to diffusion of cosmic ray electrons (CREs). The correlation breaks at a larger scale (≃320 pc) at 0.14 GHz than at 1.4 GHz (≃200 pc), indicating that the CREs traced at lower frequencies have diffused a longer path in the ISM. We find that the energy index of CREs becomes flatter in star-forming regions, in agreement with previous studies. Cooling of CREs due to the magnetic field is evident globally only after compensating for the effect of star formation activity that both accelerates CREs and amplifies magnetic fields. Compared with other nearby galaxies, we show that the smallest scale of the radio–FIR correlation is proportional to the propagation length of the CREs on which the ordered magnetic field has an important effect.

Verification of electromagnetic simulation capabilities in global gyrokinetic particle-in-cell code GTS

Recently, the numerical scheme presented by Mishchenko et al. [Phys. Plasmas 21, 052113 (2014); 21, 092110 (2014)] enabled explicit gyrokinetic simulations of low-frequency electromagnetic instabilities in tokamaks at experimentally relevant values of plasma β. This scheme resolved the long-standing cancellation problem that previously hindered gyrokinetic particle-in-cell code simulations of magnetohydrodynamic phenomena with inherently small parallel electric fields. Moreover, the scheme did not employ approximations that eliminate critical tearing-type instabilities. Here, we report on the implementation of this numerical scheme in the global gyrokinetic particle-in-cell code GTS. This implementation allows for a more complete and accurate picture of interaction between small scale turbulence and MHD modes in tokamaks. Additionally, we present a comprehensive set of verification simulations of numerous electromagnetic instabilities relevant to present-day tokamaks. These simulations encompass the kinetic ballooning mode, the internal kink mode, the tearing mode, the micro-tearing mode, and the toroidal Alfven eigenmode destabilized by energetic ions, which are all instrumental in understanding tokamak physics. We will also showcase the preliminary nonlinear simulations of kinetic ballooning instabilities and (2,1) island formation due to tearing mode instability. These simulations validate the accuracy of the scheme implementation and pave the way for studying how these instabilities affect plasma confinement and performance.