Skip to main content
Login | Suomeksi | På svenska | In English

Browsing by master's degree program "Magisterprogrammet i elementarpartikelfysik och astrofysikaliska vetenskaper"

Sort by: Order: Results:

  • Leppälä, Ari (2023)
    The thesis presents results for simulating an orbit determination process with least-squares fitting (LSF) of synthetic observations for non-cooperative satellites. Five satellites with different masses, sizes, and orbital inclinations are orbiting in Low Earth Orbit (LEO). Three different inclinations (53°, 85°, and 98°) are simulated using Orekit (Orbit Extrapolation Kit), an open-source astro- dynamical software library. Satellites are observed with four hypothetical radars. The thesis is a feasibility study. It addresses four main research questions: How reliable are predictions of the satellite position for 24 h, 48 h, 72 h, and 96 h after the first detection when using the observations from four consecutive overpasses and the LSF-estimated orbit? How accurately the Keplerian or- bital elements be modelled with LSF? How does the satellite’s mass affect the results’ reliability, and can the satellite’s drag, reflection, and absorption coefficients be estimated with Orekit? The simulation utilizes the Orekit software library in a Python environment, and the simulation incorporates several perturbing forces: Earth’s gravitational potential model EGM2008, an atmo- sphere model with space weather data NRLMSISE-00, tidal forces, and point-like masses of the Sun and the Moon. Synthetic observations were simulated by assuming Gaussian distribution for the uncertainties of the measurements. The initial estimate of the orbit is done with Gooding’s method, and the Levenberg–Marquardt algorithm is used for the LSF. The propagation of the actual satellite for 96 h was compared with the 96-h propagation of estimated orbital elements and with a satellite with estimated size and mass. Simulations were conducted with overestimated and underestimated initial masses. Each combination of a satellite, an observing radar, and an orbit was repeated 100 times. The results for reliability were promising. The likelihood of the sky-plane projected error of the position staying under the 3.0-km for up to 3 days was found to be 90% for two of the simulated radars, but only in cases of orbits having an inclination of 53°. Two other orbits had less successful results, but simulations revealed a potential increase in reliability by adjusting the observational strategy. There was a small bias in the results hinting that underestimating the satellite’s mass could result in less accurate orbit estimates than overestimating it. The satellite’s drag, absorption, and reflection coefficients could not be obtained due to the short observational time intervals simulated.
  • Hietala, Hilppa (2020)
    The aim of this thesis is to explore applications of machine learning to the study of asteroid spectra, and as such, its research question can be summarized as: How can asteroid spectra be analyzed using machine learning? The question is explored through evaluation of the obtained solutions to two tasks: the optimal locations of spectrophotometric filters for asteroid classification success and the formation of an asteroid taxonomy through unsupervised clustering. First, background theory for asteroids and particularly spectroscopy of asteroids is presented. Next, the theory of machine learning is briefly discussed, including a focus on the method utilized to solve the first task: neural networks. The first task is executed by developing an optimization algorithm that has access to a neural network that can determine the classification success rate of data samples that would be obtained using spectrophotometric filters at specific locations within the possible wavelength range. The second task, on the other hand, is evaluated through determining the optimal number of clusters for the given dataset and then developing taxonomies with the clustering algorithm k-means. The obtained results for the first task involving the optimal locations of filters for spectrophotometry seem reliable, and correlate relatively well with well-known mineralogical features on asteroid surfaces. The taxonomic systems developed by the unsupervised clustering also succeeded rather well, as many of the formed clusters seem to be meaningful and follow the trends in other asteroid taxonomies. Therefore, it seems that based on the two investigated tasks, machine learning can be applied well to asteroid spectroscopy. For future studies, larger datasets would be required for improving the overall reliability of the results.
  • Yu, Sicheng (2024)
    Aims. In this thesis, we review the history in studying the evolution path of magnetic white dwarfs and explain the longstanding questions and debates over their magnetic origin. We intend to find magnetic white dwarfs in various forms (isolated or with companions) from the spectral database of LAMOST DR7 to complement the current magnetic white dwarf catalogue. We then move on to compare our results with some commonly accepted scenarios regarding their magnetic origin. Methods. Low-resolution spectra is the main source in this study, we intend to locate signs of Zeeman-splitting in spectra of isolated white dwarfs, measure separations of substructures due to Zeeman-splitting and estimate magnetic field strength. Magnetic white dwarfs in binary or multiple systems are found by seeking signs of cyclotron radiation due to mass transfer and particle movement in magnetic fields. Photometric survey from Transiting Exoplanet Survey Satellite (TESS) was used to fold periodic light curves for targets of interest, in order to further study the nature of our candidates, especially the ones that are believed to coexist with companions. Results. We identified 31 isolated magnetic white dwarfs in the LAMOST DR7 database by the discovery of Zeeman-splitting components. Their estimated magnetic field strength ranging from below 1 mega gauss (MG) to a scale ten times larger. Two Polars/Intermediate Polars were found with both Zeeman-splitting components and broad Balmer emissions usually seen in cataclysmic variables. We also discovered two candidates of detached magnetic binaries. These systems are believed to be the progenitors of polars or intermediate polars. Despite their rarity, these candidates serve as vital hints in clarifying the ongoing debates concerning the magnetic field origins in white dwarfs.
  • Khalil, Hossam (2024)
    Understanding the baryonic physics on the galaxy group level is a prerequisite for cosmological studies of large-scale structures. While the majority of baryons in galaxy groups are located in their intragroup medium (IGrM), one poorly understood aspect of galaxy groups is their hot intragroup X-ray emission. In this thesis, a new all-sky catalogue of X-ray detected groups (AXES-2MRS) is presented, based on the identification of large X-ray sources discovered in the ROSAT All-Sky Survey (RASS) with the 2MRS Bayesian Group Catalogue. In addition to X-ray luminosity coming from the shallow survey data of RASS, detailed X-ray properties of the groups have been obtained by matching the catalogue to archival X-ray observations conducted by XMM-Newton. The relationship between X-ray and optical properties of AXES-2MRS is explored through scaling relations, namely $\sigma_{v}-L_{X}$, $\sigma_{v}-kT$, $\sigma_{v}-M$, and $kT-L_{X}$ which denote (velocity dispersion vs. X-ray luminosity), (velocity dispersion vs. X-ray temperature), (velocity dispersion vs. hydrostatic mass), (X-ray temperature vs. X-ray luminosity), respectively. The scaling relations reveal similarities between our low-redshift catalogue and high-redshift studies implying that our knowledge about galaxy groups is redshift-invariant. This study enhances the representation of the underexplored low-z, low-luminosity galaxy groups, particularly in low-mass systems ($< 10^{14} M_{\odot}$). This enhances the completeness of galaxy group catalogs, addressing the persistent issue of missing faint, low-mass systems. Moreover, previous catalogues, based on detecting the peak of the X-ray emission preferentially sample the high dark matter (DM) halo-concentration groups, while AXES-2MRS includes many low DM halo-concentration groups.
  • Nincă, Ilona Ştefana (2020)
    Cadmium Telluride (CdTe) has a high quantum efficiency and a bandgap of 1.44 eV. As a consequence, it is being used to efficiently detect gamma rays. The aim of this thesis is to explore the properties of the CdTe pixelated detector and the procedures conducted in order to fine-tune the electronic readout system. A fully functional CdTe detector would be useful in medical imaging techniques such as Boron Neutron Capture Therapy (BNCT). BNCT requires a detector with a good energy resolution, a good timing resolution and a good stopping power. Although the CdTe crystal is a promising material, its growing process is difficult due to the fact that different types of defects appear inside the crystal. The quality assurance process has to be thorough in order for suitable crystals to be found. An aluminum oxide layer (Al2O3) was passivated onto the surface of the crystal. The contacts for both sides were created using Titanium Tungsten (TiW) and gold (Au) sputtering deposition, followed by an electroless nickel growth. I tested the CdTe pixelated detector with different radioactive sources such as Am-241, Ba-133, Co-57, Cs-137 and X-ray quality series in order to study the sensitivity of the device and its capacity to detect gamma and X-rays.
  • Arvo, Jukka (2023)
    Tutkielmani käsittelee kullan pysyvän isotoopin Au-197 tuottamista neutroniaktivaatiolla luonnollisesta elohopeanäytteestä. Kokeen kannalta pääasiallinen reaktiomenetelmä oli Hg-196 neutronikaappaus. Kyseinen transmutaatio suoritettiin myös kokeellisesti. Elohopeaa sisältävänä näytteenä käytettiin Ardentin valmistamia Futura Standard -hammasamalgaamikapseleita. Turvallisuussyistä kapselit olivat koejärjestelyssä alkuperäisessä purkissaan. Kapseleita oli kaikkiaan 50 kappaletta, ja jokaisessa oli 400 mg elohopeaa. Yhteensä näytteessä oli siis 20 grammaa elohopeaa. Näytettä säteilytettiin STUKin tiloissa 14 vuorokauden ajan kolmen AmBe-neutronilähteen avulla. Valmistajan ilmoittamat neutronituotot käytetyille lähteille ovat 2.0E+7 n/s, 2.1E+6 n/s ja 6.7E+5 n/s. Lähteiden ilmoitetut aktiivisuudet ovat vastaavasti 333 GBq, 37.0 GBq ja 11.1 GBq. Neutronien hidastamiseen käytettiin HDPE-tankoa. Säteilytyksen jälkeen näytettä mitattiin STUKin gammaspektrometrian laboratorion B6 p-HPGe BE5030 -germaniumilmaisimella, ja kullan synty voitiin todentaa spektristä löytyvien karakterististen gamma- ja röntgenpiikkien avulla. Koejärjestely onnistui, ja työni osoittaa, että kullan pysyvän isotoopin Au-197 valmistaminen luonnollisesta elohopeasta havaittavissa määrin on mahdollista toteuttaa melko yksinkertaisella koejärjestelyllä käyttäen varsin pienitehoisia neutronilähteitä. Varsinaisen kokeen lisäksi käsittelen työssäni myös kullanteon historiaa sekä aiheeseen liittyvää teoriaa.
  • Vuori, Mikko (2023)
    A method for deriving the complex refractive index of a mm-sized single particle in a specific wavelength using laboratory measurements is presented. Laboratory measurements were done using the 4π scatterometer, which measures Mueller matrix elements of a particle suspended in air using acoustic levitation as a function of scattering angle. To obtain the complex refractive index of the particle, measurements were compared to simulations from a newly developed SIRIS4 Fixed Orientation (SIRIS4 FO) geometric optics simulation. The 4π scatterometer is a unique instrument which measures Mueller matrix elements from a particle using linear polarizers and a detector rotating about the particle on a rotational stage. The scatterometer uses an acoustic levitator as a sample holder which provides nondestructive measurements and full orientation control of the sample. To compare the measurement results to simulations, SIRIS4 single-particle geometric optics code was modified to handle particles in a fixed orientation. The original code is able to calculate the Mueller matrix elements for a given 3D model, but averages the results over the orientation of the particle. The modified SIRIS4 FO calculates the Mueller matrix elements over the full solid angle as functions of the two scattering angles, which give the direction of observation of the scattered wave compared to the direction of the incident wave. A 3D model of the shape of the measured particle was constructed using X-ray microtomography, and was translated to SIRIS4 FO. The complex refractive index was obtained with a nonlinear least squares analysis by minimizing the sum of squared residuals between the measurements and simulations with varying refractive index values. Finally, confidence regions were constrained for the results, by estimating the computed residuals between simulations and measurements as the random errors in the nonlinear model.
  • Blomqvist, Sofia (2024)
    The matter in neutron stars exist under extreme conditions, and the cores of these stars harbour densities unreachable in any laboratory setting. Therefore, this unique environment provides an exceptional opportunity to investigate high-density matter, described by the theory of Quantum Chromodynamics (QCD). This thesis centers on the exploration of twin stars, hypothetical compact objects that extend beyond the neutron star sequence. Originating from a first-order phase transition between hadronic matter and quark matter, our focus is on understanding the constraints on these phase transitions and their effect on the observable properties of twin stars. In our investigation of twin stars, we construct a large ensemble of possible equations of state featuring a strong first-order phase transition. We approximate the low- and high-density regions with polytropic form and connect them to chiral effective field theory results at nuclear densities and extrapolated perturbative QCD at high densities. The resulting equations of state are then subjected to astrophysical constraints obtained from high-mass pulsars and gravitational wave detections to verify their compatibility with observations. Within our simple study, we identify two distinct types of twin stars, each providing a clear signature in macroscopic observables. These solutions originate from separate regions in the parameter space, with both regions being relatively small. Twin stars in our approach generally obtain small maximum masses, while the part of the sequence corresponding to neutron stars extends to large radii, indicating that these solutions only marginally pass the astrophysical constraints. Finally, we find that all twin stars obtain sizable cores of quark matter.
  • Lipsanen, Veera (2024)
    The constant outflow of solar wind from the Sun and possible larger structures within it influence the Earth's magnetosphere. These large structures include interplanetary coronal mass ejections (ICMEs) and high speed streams (HSSs). They can contain substructures: fast enough ICMEs can have a turbulent sheath region in front of them, while a HSS can interact with the slower ambient solar wind and form a stream interaction region (SIR). Pc5 Ultra-low frequency (ULF) waves have a frequency range of 2–7 mHz and they are important in energy transfer from solar wind to the magnetosphere and they affect energetic electrons in the radiation belts. ULF waves in the magnetosphere are generated by multiple mechanisms. For example fluctuations in solar wind's dynamic pressure create waves on the dayside, Kelvin-Helmholtz instability, often caused by HSSs, on the magnetopause flanks and substorms on the nightside. This makes ULF waves MLT dependent. For this thesis a new ground-based ULF index that is MLT dependent and has a resolution of 1 min is constructed using the wavelet analysis method. The aim of this thesis is to study how this new index correlates with multiple solar wind parameters, geomagnetic indices and an already existing ULF index during substructures of four events: weak HSS and ICME and strong HSS and ICME. The ULF power is found to peak at the sheath–ejecta boundary during ICMEs and the stream interface during HSSs, primarily driven by dawn and night ULF powers. The AE index is found to correlate with ULF power during all of the events, which indicates that even the non-geoeffective event produces some kind of substorm activity. Solar wind speed is found to correlate well with the ULF power during SIRs. It is important to take into account the MLT dependence of ULF waves since their generation mechanisms are different at different parts of the magnetosphere. In addition, we found that the ULF power in the four MLT sectors can behave differently at the same moment.
  • Al-Adulrazzaq, Aula (2023)
    Dark matter direct detection experiments still have found no evidence of the dark matter WIMPs. The search has therefore been expanded for lighter dark matter candidates. Light dark matter is nearly invisible to current detectors through the elastic nuclear recoils. This thesis is meant to provide understanding on the inelastic atomic scatterings, which are one good way to detect dark matter particles with mχ ∼ GeV. In this thesis we consider spin-independent scatterings. Inelastic scatterings are based on the fact that in an atom, electrons do not follow the motion of the recoil nucleus immediately, but instead it takes time. This results in a small probability of observable ionization or excitation of the atom. This is known as the Migdal effect. We will first study the theoretical framework of dark matter-nucleus scatterings, showing how to get the event rate and how it is factorized into the astrophysical, the particle physics and the target response part. Then we will move to the inelastic processes, Migdal and Bremsstrahlung effects, deriving their event rates. In the first, we try to detect ionized electrons. The latter one, the Bremsstrahlung, is a similar process to the Migdal, but there we try to detect photons emitted from the de-excitations of atoms excited in the inelastic recoils. We will also look into the Migdal in semiconductors. Because of the smaller gap for electron excitations in crystals, we find that the rate for the Migdal effect is much higher in semiconductors than in atomic targets, thus allowing the search for even lighter dark matter particles. The rate can be expressed in terms of the energy loss function of the target material.
  • Nurminen, Niilo Waltteri (2021)
    Phase transitions in the early Universe and in condensed matter physics are active fields of research. During these transitions, objects such as topological solitons and defects are produced by the breaking of symmetry. Studying such objects more thoroughly could shed light on some of the modern problems in cosmology such as baryogenesis and explain many aspects in materials research. One example of such topological solitons are the (1+1) dimensional kinks and their respective higher dimensional domain walls. The dynamics of kink collisions are complicated and very sensitive to initial conditions. Making accurate predictions within such a system has proven to be difficult, and research has been conducted since the 70s. Especially difficult is predicting the location of resonance windows and giving a proper theoretical explanation for such a structure. Deeper understanding of these objects is interesting in its own right but can also bring insight in predicting their possibly generated cosmological signatures. In this thesis we have summarized the common field theoretic tools and methods for the analytic treatment of kinks. Homotopy theory and its applications are also covered in the context of classifying topological solitons and defects. We present our numerical simulation scheme and results on kink-antikink and kink-impurity collisions in the $\phi^4$ model. Kink-antikink pair production from a wobbling kink is also studied, in which case we found that the separation velocity of the produced kink-antikink pair is directly correlated with the excitation amplitude of the wobbling kink. Direct annihilation of the produced pair was also observed. We modify the $\phi^4$ model by adding a small linear term $\delta \phi^3$, which modifies the kinks into accelerating bubble walls. The collision dynamics and pair production of these objects are explored with the same simulation methods. We observe multiple new effects in kink-antikink collisions, such as potentially perpetual bouncing and faster bion formation in comparison to the $\phi^4$ model. We also showed that the $\delta$ term defines the preferred vacuum by inevitably annihilating any kink-antikink pair. During pair production we noticed a momentum transfer between the produced bion and the original kink and that direct annihilation seems unlikely in such processes. For wobbling kink - impurity collisions we found an asymmetric spectral wall. Future research prospects and potential expansions for our analysis are also discussed.
  • Lehtinen, Simo (2021)
    The solar corona constantly emits a flow of charged particles, called the solar wind, into interplanetary space. This flow is diverted around the Earth by the magnetic pressure of the Earth’s own geomagnetic field, shielding the Earth from the effect of this particle radiation. On occasion the Sun ejects a large amount of plasma outwards from the corona in an event called a Coronal Mass Ejection (CME). Such events can drive discontinuities in the solar wind plasma, called interplanetary shocks. Shocks can affect the Earth’s magnetosphere, compressing it inwards and generating electromagnetic waves inside it. In this thesis we will cover a study of the ultra-low frequency (ULF) wave response in the magnetosphere to CME-driven shocks. Geomagnetic pulsations are ultra-low frequency plasma waves in the magnetosphere, observable from ground-based magnetometers. The compression of the magnetosphere by interplanetary shocks generates geomagnetic pulsations in the Pc4 and Pc5 frequency ranges (2 - 22 mHz). These waves play an important role in magnetospheric dynamics and the acceleration and depletion of high energy electrons in the radiation belts. We consider 39 interplanetary shock events driven by CMEs, and analyse ground-based magnetometer data from stations located near local noon at the time of the shock arrival. Solar wind measurements are used to categorise interplanetary shocks based on their Mach number and the dynamic pressure differential as main indicators of shock strength. The importance of these parameters in determining the strength of the wave response in the geomagnetic field is then studied using wavelet analysis and superposed epoch analysis. Stronger shocks are found to result in larger increases in wave activity, especially in the Pc4 range. Ground stations at higher latitudes observe higher wavepower, but there is an interesting anomaly in the Pc4 range at stations magnetically connected to regions near the plasmapause, which show an enhanced wavepower response. We quantify the decay time of the wave activity and find that it is around 20 hours for Pc5 waves and 7 hours for Pc4 waves.
  • Mozejko, Arik (2023)
    Dark matter (DM) is introduced and explored in a holistic perspective. Topics include observational evidence, various DM properties, potential candidates, and the tenets of indirect versus direct DM detection. Then an emphasis is placed on understanding the cryogenic detection of weakly interacting massive particles, with explicit connection to phonon-based detection of DM. The importance of improving methods of DM direct detection are emphasised, with specifically the usage of molecular dynamics simulations as an avenue of studying defect creation in cryogenic detector materials. Previous investigations into this area are reviewed and expanded upon through novel experimentation into how defect properties vary when changing thermal motion of the crystal lattice. This experimentation is conducted via the usage of molecular dynamics simulations on sapphire (Al2O3) as a DM direct detection material, and it is found that while atomic velocity does not impact the overall emergent defect structure, it does have an impact on the energy lost in these defects. Changing the temperature of the lattice produces the expected results, generating greater variance in both defect band structure as well as average energy loss.
  • Gonzalez Ateca, Marcos (2020)
    The distribution of matter in space is not homogeneous. Large structures such as galaxy groups, clusters or big empty spaces called voids can be observed at large scales in the Universe. The large scale structure of the Universe will depend on both the cosmological parameters and the dynamics of galaxy formation and evolution. One of the main observables that allow us to quantify this structure is the two-point correlation function, with which we can trace different galaxy properties such as luminosity, stellar mass and also, it enables us to track its evolution with redshift. In galaxy surveys, we do not obtain the location of galaxies in real space. We obtain our data in what it is called redshift space. This redshift space can be defined as a distortion of the real space generated by the redshift introduced by the peculiar velocities of galaxies and from the Hubble expansion of the Universe. Therefore, the distribution of galaxies in redshift space will look different from the one obtained in real space. These differences between both spaces are small but not negligible, and they depend strictly on the cosmology. In this work, we will assume a ΛCDM cosmology. Therefore, in order to find the different 1-dimensional or 2-dimensional correlations functions, we will use the most updated version of the code provided by the Euclid consortium, which belongs officially to the ESA Euclid mission. Moreover, we will also need different galaxy catalogues. These catalogues have already been simulated and they are called Minerva mocks, which are a set of 300 different cosmological mocks produced with N-body simulations. Finally, as there is a well-defined relation between real and redshift space, one could also assume that there is a relation between the two-point correlation functions in both real and redshift space. In this project, we will prove that the real-space one-dimensional two-point correlation function, which is the physically meaningful one, can be derived from the two-dimensional two-point correlation function in redshift space following a geometrical procedure independent of approximations. This method, in theory, should work for all distance scales.
  • Gibson, Natalie (2023)
    The search for a profound connection between gravity and quantum mechanics has been a longstanding goal in theoretical physics. One such connection is known as the holographic principle, which suggests that the dynamics within a given region of spacetime can be fully described on its boundary surface. This concept led to the realization that string theory provides a lower-dimensional description that encapsulates essential aspects of spacetime. While the "AdS/CFT correspondence" exemplifies the success of this holographic theory, it was discovered soon after that the Universe has a positive cosmological constant, Λ. This immediately sparked interest in a potential correspondence centered around de Sitter (dS) space, which is also characterized by a positive cosmological constant. This thesis comprehensively explores the de Sitter/Conformal Field Theory (dS/CFT) correspondence from various perspectives, along with the unique challenges posed by the distinct nature of dS space. The original dS/CFT duality proposes that a two-dimensional Conformal Field Theory resides on the boundary of three-dimensional asymptotic dS space. However, the definition and interpretation of physical observables within the dS/CFT framework remain open questions. Therefore, the discussions in this thesis not only cover the original dS/CFT conjecture, but also extend into more recent advancements in the field. These advancements include a higher-spin dS/CFT duality, the relationship between string theory and dS space, and the intriguing proposal of an "elliptical" dS space. While the dS/CFT correspondence is still far from being well-defined, there have been extensive efforts devoted to shedding light on its intricate framework and exploring its potential applications. As the Universe may be evolving towards an approximately de Sitter phase, understanding the dS/CFT correspondence offers a unique opportunity for gaining fresh insights into the link between gravity and quantum field theory.
  • Häkkinen, Jenni (2024)
    Gravitational waves from cosmological phase transitions are a promising probe of the early universe. Many theories beyond the Standard Model predict the early universe to have undergone a cosmological first-order phase transition at the electroweak scale. This transition would have produced gravitational waves potentially detectable with the future space-based detector Laser Interferometer Space Antenna (LISA). We study the gravitational wave power spectrum generated by sound waves, which are a dominant source of gravitational waves from first-order phase transitions. We compare two methods for calculating the sound wave power spectrum: a simulation-motivated broken power-law fit of the shape of the spectrum, and a wider theoretical framework called the Sound Shell Model, which includes hydrodynamic calculations of the phase transition. We present an implementation of the Sound Shell Model into the PTPlot tool, which is currently based on the broken power-law fit. With PTPlot, we calculate the signal-to-noise ratios of LISA for the sound wave power spectrum of each method. The signal-to-noise ratio allows us to estimate the detectability of gravitational wave signals with LISA. We analyse how the detectability of certain particle physics models changes between the two different methods. Our results show that the Sound Shell Model has a potentially significant impact on the signal-to-noise ratio predictions, but it does not uniformly improve or worsen the detectability of the gravitational wave signals compared to the broken power law. The code implementation is overall successful and lays the foundation for an updated release of PTPlot and future work within this topic.
  • Takala, Saara (2024)
    Ultra-low frequency (ULF) waves in the Pc4-Pc5, 2 – 25 mHz range have been observed to accelerate trapped 1 – 10 MeV electrons in the Earth’s radiation belts. This acceleration can lead to particle losses and injections that occur on timescales comparable to the particle drift periods. Current models rely on diffusion equations written in terms of Fokker-Planck equations and are not suitable for describing fast temporal variations in the distribution function. This thesis is a study of fast transport of equatorially trapped electrons in the radiation belts. We look at solutions for the time evolution of the linear part of the perturbed distribution function using both analytical and numerical methods. Based on this work we build a simple model of fast transport in the radiation belts using a spectral PDE framework called Dedalus. The resulting program is a computationally inexpensive, simple approach to modelling drift-periodic signatures on fast timescales. In this study we investigate the behavior of the distribution function in three systems: a simple system without a wave term, and systems with a single non-resonant and resonant ULF wave. The wave solutions are evaluated with magnetic field perturbations of different magnitudes. The Earth’s magnetic field is modelled with the Mead field. The numerical solution of the perturbed differential equation is studied for relativistic equatorially trapped electrons. Phase-mixing is found to happen regardless of field fluctuations or resonance. The non-resonant wave solution shows time-delayed, spatially localized structures in the equatorial plane forming in the presence of large magnetic field fluctuations. These transients are also seen in the analytical solution and provide a new theoretical explanation for the ubiquitous observation of drift echoes in the inner and outer radiation belts of the Earth (Li et al., 2024).
  • Suni, Jonas (2021)
    Magnetosheath jets are a class of structures in the Earth's magnetosheath usually defined by an enhancement of the dynamic pressure of the plasma. Magnetosheath jets have been observed by several different spacecraft over the past few decades, but their origin and formation mechanism have remained unclear. The aim of this thesis is to use data from a global simulation to investigate the origin of magnetosheath jets. We defined two different kinds of structures, magnetosheath jets and foreshock compressive structures (FCS), and collected a database of individual jets and FCSs from 4 Vlasiator global hybrid-Vlasov simulation runs, all of which simulate only the ecliptic plane. We then conducted a statistical analysis of the properties of jets and FCSs, and their occurrence rates as a function of the definition of the FCS criterion. Jets were separated into two categories: jets that form in contact with FCSs (FCS-jets), and those that do not (non-FCS-jets). We found that up to 75% of magnetosheath jets form in association with an FCS impacting the Earth's bow shock. We also found that FCS-jets penetrate deeper into the magnetosheath than non-FCS-jets. Finally, we found no conclusive explanation for the formation of non-FCS-jets. The properties of both jets and FCSs agree qualitatively and to some extent quantitatively with spacecraft observations and other simulations in the literature. The formation of jets from FCSs impacting the bow shock is similar to the proposed theory that jets are linked to Short Large-Amplitude Magnetic Structures (SLAMS). In the future, we will study magnetosheath jets and FCSs in polar plane simulation runs as well, and ultimately in full 3D simulation runs. If made possible by new simulations, the effects of electron kinetic effects on jets and FCSs will also be studied. Comparison studies with spacecraft observations of jet formation from FCSs will also be conducted, if and when such observations are found and become available.
  • Kostamo, Iida (2023)
    In this thesis, 16 gas-free galaxy merger simulations were run in order to determine the dominant effect that causes the formation of cores, i.e. regions of ”missing light”, in the centers of massive early-type galaxies. The simulations were run on the Mahti supercomputer at the Finnish IT Centre for Science (CSC), using the simulation codes GADGET-3 and KETJU. The merging of galaxies will eventually lead to the merging of their central supermassive black holes (SMBHs). The evolution of a SMBH binary can be divided into three phases: the dynamical friction phase, the three-body interaction phase and gravitational wave (GW) emission phase. After the GW emission phase, the merged SMBH may receive a recoil velocity. To study the effects of these three phases on the formation of the core, three kinds of simulations were run. These include three GADGET runs and thirteen KETJU runs that can be divided into two groups, since eight of the runs had GW recoils enabled. When using only GADGET, the three-body interactions are not modeled due to the softening of gravity. Using KETJU allows for modeling the later evolution of the SMBH binary, including the three-body interaction phase, the GW emission phase and the resulting GW recoil. The initial conditions for the progenitor galaxies were motivated by the galaxy NGC 1600. The same stellar and dark matter profiles were used for each simulated galaxy. Two different SMBH masses were used. One KETJU merger was also run without SMBHs. For the runs with spinning SMBHs, the spin directions and magnitudes were chosen in such a way that different magnitudes of recoil velocity would be achieved. The size of the core can be determined from the brightness profile of the galaxy, assuming a constant mass-to-light ratio. The commonly used core-Sérsic profile was fit to the surface mass density profiles of the merger remnant galaxies. The best-fit parameter values were then used to estimate the sizes of the cores, such as the core radii. The amount of missing mass in the centers of the galaxies, i.e. the mass deficits, were also computed based on the core-Sérsic fits. We found that the mass deficit correlates positively with the mass of the SMBH for all the KETJU runs. Including GW recoils in the simulations was found to increase the mass deficits by roughly the equivalent of one SMBH mass compared to the KETJU runs without the recoils. The formation of cores was the weakest for the GADGET runs, since the cores were created only by the large-scale dynamics. No core was formed in the run without SMBHs, as expected. Thus, we can conclude that SMBHs are essential for the formation of cores in massive early type galaxies, and the largest cores are formed when GW recoils are included in the model.
  • Grön, Julia (2023)
    Mars is a rocky planet in the Solar System, fairly similar to the Earth. It is known for its red color and the theory that there has been liquid water and possibly life in Mars at some point in its history. Mars has been a target of study since space exploration began, with the first fly-by mission occuring in 1965. In the past, Martian climate could have supported a functioning hydrological system and Mars even could have had an ocean. This system would have been similar to the one happening in Antarctica due to cold weather. The Martian atmosphere is mostly composed of carbon dioxide, much thinner and colder compared to the Earth's atmosphere. Methane is a main point of interest due to it being a possible biosignature, a sign of life. Another interesting feature about the atmosphere are the dust devils that contibute to the climate by lifting dust into the air. The Martian soil is composed of rocks and dust containing toxic perchlorates. To colonize Mars, several requirements need to be met: the transportation and funding, designing the settlement efficiently with everything needed and determining how many people are necessary. The risk assessment to both people and equipment needs to be made and taken into account.