When looking up to the sky and trying to detect the photons of far away stars that are still in formation, we can get access to knowledge we never had before. In order to do this, we must make sure we truly understand the nature of how photons behave in accordance to each other and to our detectors.In quantum optics we can see that these photons have a very dubious character. They are both particle and wave and neither. This makes studying how they are distributed over the spectrum of the submillimeter and millimeter wave astronomy a very big challenge. We know that this very particlewave duality is what makes up the noise that is fundamental to photons. This duality also makes them obey the BoseEinstein statistics. Hence due to their particlelike behavior, we have a noise component that resembles that of a Poisson distribution. Very much like raindrops falling down from the sky. On the other hand we have the wave nature of the photons, this causes them to arrive in bunches rather than these random raindrops. It is observed that photons arriving at a detector are corrrelated.
In this thesis we will be investigating how to incorporate the noise of the bunching of the photons into the model of TiEMPO. This model simulates the signal processing of a measurement done by the wideband spectrometer named DESHIMA. Due to the fact that DESHIMA operates in a wideband frequency range, we do theoretical research that explains the fundamental theory behind calculating the photon noise over this wideband range. We show that taking the wideband integral of the photon noise is mathematically equivalent of summing the narrowband approximation for infinitely many subbands and adding them up top each other. This approach is the previous method of calculating the photon noise over the wideband. Due to this method being valid, the question of how these variances over these smaller subbands can be additive?
Since we are dealing with the detection of photons which as previously stated is a correlated signal. By modeling a simplified version of wideband photon detection, we have come to the conclusion that due to the small coherence time these photons are independent in the wideband signal. The photons in these smaller subbands of the wideband signal can also be viewed as statistically independent. If we decrease the frequency bandwidth, we increase the coherence time. Thus measuring the signal over this subband equates to having a larger uncertainty in time. Hence when a photon is detected in this subband, due to the large coherence time we have that the knowledge of when this photon arrived is mostly lost. Making the time correlations irrelevant to a measurement of an integration time this long.

Quantum Markov Semigroups (QMS) describe the evolution of a quantum system by evolving a projection or density operator in time. QMS are generated by a generator obeying the well-known Lindblad equation. However, this is a difficult equation. Therefore, the result that the Lindblad form greatly simplifies in the case of the generator commuting with the modular automorphisms group, is useful. Unfortunately, the proof only works for finite dimensional Hilbert spaces, which is why the aim of this thesis is to generalise this result to countably infinite dimensional Hilbert spaces. To this end, the Lindblad equation is derived from both a mathematical and physical perspective. Where the former relies on rigorous proof and the latter relies on approximations.   In the rigorous case the theory of unital completely positive maps is used. Furthermore, multiple topologies are considered which put less stringent conditions on the operators of interest than the norm topology. Additionally, the Haar measure is used on the unitaries of the bounded linear operators to construct the explicit Lindblad form. To derive the result by employing physical assumptions the interaction picture is used. The physical derivation starts from the Von Neumann equation and uses multiple assumptions to obtain the final Lindblad form. The most important physical assumptions are: the Born approximation, the Markov approximation and the rotating wave approximation.   Furthermore, the main result is the generalisation of the simplified Lindblad form. This simplified form holds for generators commuting with the modular automorphisms group in case the Hilbert spaces are countably infinite dimensional. However, this requires the domain of the generator to be restricted to trace class operators with the identity operator artificially added. Additionally, the generator needs to map strongly convergent sequences to weakly convergent sequences. It also needs to be self-adjoint with respect to the Hilbert-Schmidt inner product. Lastly, the generator is assumed to be self-adjoint with respect to the Gelfand-Naimark-Segal (GNS) inner product <X, Y>=Tr(σ X*Y) for σ a density operator. This last assumption implies that the generator commutes with the modular automorphisms group, which is the symmetry we are considering. Hence, the two previous assumptions are the additional requirements needed to generalise the result, besides the restriction of the domain. Therefore, it is recommended for further research to generalise the result for the domain extended to the bounded operators B(H). It should be noted that the proof heavily relies on the Hilbert space structure induced by the Hilbert-Schmidt inner product. Consequently, the generalisation for the bounded operators would probably require a different approach. Another recommendation is to try and lift the sequence and self-adjoint requirements on the generator. In addition, it is interesting to investigate which physical systems actually have the symmetry of generators commuting with the modular automorphisms group.

The goal of this thesis is to implement and experiment with a Bayesian way of estimating a (smooth) monotone regression function by applying it to climate data. The method we use is proposed by Neelon and Dunson. This method uses a piece-wise linear model for the unknown regression function and enforces the monotonicity constraint by the specification of the prior distribution of the slopes. This thesis is also aimed at providing solutions to specific problems that we encounter during the process of applying this method. We encounter two main problems: a numerical problem and a boundary problem. The numerical problem concerns a fraction of very small numbers, which we can solve using an asymptotic approximation of the Mills Ratio. The boundary problem appears in the form of a steep upward-sloping curve at the left boundary. We provide a pragmatic solution to this boundary effect, in which we use an extension of the data set, to obtain a better curve estimate at this boundary.

In this bachelor thesis we use a stochastic model to aspire to explain biodiversity patterns in different ecosystems with selection advantage. The stochastic model we use is an extension of the mean­-field voter model where we include a selection factor. In the model individuals with two different types of alleles in two different ecosystems are considered. The model is a stochastic Markov process that describes interactions of individuals with each­other over time. This means that the ratio of individuals with certain alleles stochastically drifts over time. The main goal of this bachelor thesis is investigate whether is it possible that individuals with two different types of alleles can coexist (a stable equi­librium) in two populations. We do this by taking the limit of this Markov process such that we can show convergence to ordinary differential equations. By studying these differential equations we ob­tain results: vector fields with equilibrium points. We conclude that, under certain conditions and with a selection advantage, coexistence of individuals with two different alleles in two different ecosystems is possible (see Section 4.4) in the form of a stable equilibrium. Furthermore we claim that the typical time to absorption, reaching an absorbing state where all the individuals have the same allele, of the two­-dimensional mean­-field voter model with selection scales exponentially with the system size.

In this report, the conversion from spin-echo signals, obtained with a low-field hand-held MRI scanner that was designed and built at the Leiden University, to images of the proton density within the sample is considered. This scanner does not make use of switchable gradient coils, but instead relies solely on the natural inhomogeneity of the field and on translations of the sample over this field for its spatial encoding. Specifically, an attempt is made to answer the question of how we can reconstruct a phantom using this kind of scanner. This is done by deriving a signal model, discretising it and writing it as a linear least squares problem. Then, we can make use of the techniques of Cojugate Gradient for Least Squares (CGLS) wih `2-regularization and Generalized Conjugate Gradient Minimal Error (GCGME) with `1-regularization for the difference between neighbouring pixels in order to solve this inverse problem. Firstly, we theoretically consider combinations for magnetic field geometry and measurement strategy for their usability for image reconstruction. After this, the obtained strategies are tested in three experiments, with two magnets and two samples. We start by doing this numerically, using a simulated phantom in combination with a measured magnetic field and the translation strategy. By doing this, we can determine if reconstruction is possible using that combination of field and strategy. Finally, the strategy is tested on real samples. Using numerical phantoms in combination with the magnetic field and translation strategy used in the measurements, we were able to correctly reconstruct the phantoms. However, the reconstruction broke down when data from real samples was considered. A variety of possible improvements is discussed. The improvement that would have the most impact would be to design a magnet that has a less uniform gradient in the z-direction, and instead has some locations in the xy-plane where the field falls slowly as function of z, but quickly in other locations.

Peer assessment (PA) is a process that involves students assessing each other’s work and providing their peers with (elaborated) feedback. This thesis investigates how an elaborated form of formative PA can be implemented in the mathematics class to promote students’ achievement in mathematics. A mixed methods investigation at a third-grade mathematics class was conducted. For this investigation, a teacher’s guide and a student guideline were developed. The teacher’s guide helps the mathematics teacher in implementing PA in the class. The students’ guideline was developed to guide students doing online PA during distance learning due to the COVID-19 pandemic. Using this student guideline can require more effort for the teacher to monitor students’ participation and more effort for students to do elaborated online PA. Students were trained with the technique “two stars and a wish” for providing peer feedback. With this feedback technique, students could provide their peers with two reinforcing feedback and one suggestive feedback. The feedback has to first be written in their peer’s work and later be elaborated when the work is being returned. As a result of this investigation, implementing PA promoted students to interact and learn from one another during online distance learning. Evidence was found that by implementing PA, students were more involved in their learning process, and students’ achievement increased. PA can be used as an effective learning method for students. However, it is preferably for PA to be implemented during physical class, especially when it is implemented the first time.

Knowledge of particle deposition is important in clinical settings or when discussing environmental effects of aerosols on humans.
Particle deposition in the human respiratory tract is determined by breathing patterns and lung morphology, as well as particle properties and deposition mechanisms.
In this study we develop a 1-dimensional model that numerically solves the general dynamic aerosol equation in the human respiratory tract.
The model can be used to calculate deposition fractions for a range of initial parameters.
We use Weibel's morphometric model to describe the lung geometry.

The model is validated by comparing it with previous numerical results, and running sensitivity tests to examine its consistency with parametric variations.
The model proved to be computationally efficient, requiring just seconds to run a simulation.
We use this to perform a number of parametric studies, most notably changing the tidal volume and the breathing rate.
For both of these, an increase in either the volume or the rate decreased the deposition fraction across the spectrum of particle sizes, apart from at the tails of the distribution.
We also examine the effect of particle density on the deposition fraction, which increases with an increasing density.
The source code is published along with this thesis, allowing anyone to perform arbitrary parametric studies of their own.

This thesis focuses on optimising the Emergency Medical Services (EMS) in the Netherlands. In the current situation, the Netherlands is divided into 24 independent regions in which ambulances operate. These regions can be merged in order to reduce response time and increase efficiency. Different models are presented in which these regions are merged in an optimal way. The thesis starts with an explanation of the current regional system within the Netherlands. Next, two models following from existing literature are discussed. These two models are used as a basis for new models that optimise the merging of regions. Then, the results are discussed, which are in line with the assumption that merging regions leads to a better coverage and a reduction in the number of ambulances needed.

Keywords:

The Abelian sandpile model was first introduced by Bak, Tang and Wiesenfeld in 1987. Since then, a lot of researchers have studied this model and similar models, all related by the concept of self-organized criticality. In this thesis, we study a variant on the classical model where dissipative and anti-dissipative vertices are incorporated in the model. These have an influence on the critical behaviour of the model. We first introduce a definition of criticality in this model and investigate which levels of dissipation are required to guarantee non-critical behaviour. In studying this variant, we encounter random walks and Green's functions.

A model is developed to simulate NaCl concentrations across a cation exchange membrane in the transient regime. The convectionless Nersnt-Planck equation is combined with the continuity equation to obtain a partial dierential equation in time and space. The boundary conditions at the edges of the cell for this equationare a constant bulk concentration. Furthermore, in presence of an applied potential, the potential gradient is assumed to be linear with dierent values for the gradient inside and outside the membrane. This dierent value inside the membrane is due to Donnan potentials that arise at the membrane/electrolyteinterfaces. On top of this, because of the Donnan exclusion, the concentrationsare not continuous on the interface. To account for this, the cell is splitinto three regions: the membrane and two electrolyte regions on either side of the membrane. The boundary conditions to tie these regions together, is that the ux should be continuous across the interface. The problem has the form of an initial value problem. The partial differential equation is discretized using a finite difference method and is integrated using the forward Euler method.The results show an accumulation of ions at the membrane interface closest to the positive electrode, while ions are depleted at the other interface. This accumulation of ions is too big however, with the Na+ concentration at one point reaching 12M after only 0.10 seconds, while the concentration was 1M atthe start. This is clearly not what one would expect to happen physically. Therewere, however, no errors found in the numerical method as it was possible tocorrectly deduce and predict concentration changes, based on the concentrations and ux in a certain time step. Therefore it appears that there is an error in the assumptions that were made when developing the model. It is believed that this error lies in the assumption that the potential gradients only has two different values. This assumption ignores the diffusion boundary layers on the outside of the membrane. In this diffusion boundary layer, the potential gradient has another, different value. Therefore, to obtain more accurate results, the cell should actually be split into five different regions: the membrane, two diffusion boundary layers and two electrolyte solutions.

In this thesis, the methods of Matrix Product States are examined. We give background in the origin of the methods and explain in detail how the method can be applied to a series of Josephson junctions. Numerical results are given for the ground state energy analysis of these series.

We have explored the concept of mobile computational fluid dynamics (CFD) solving. Using the computational
power of a smartphone we tried to generate real time simulations of relatively simple transient natural
convection problems in the laminar regime. The focus lies on finding a compromise between accuracy, speed
and stability: we want to make a real timemobile CFD solver that is as accurate as possible.
A JAVA application has been created that runs on Android. The SIMPLE algorithm has been implemented
in order to solve for the flow and heat transfer. The SIMPLE algorithm on the application was tested for
runtime performance and accuracy.
For a test problem, a real time simulation on a Nexus 5X has been realized with an accuracy of 5.4% on a
21x21 grid. Running the algorithm on a desktop gave a simulation 10 times faster than real time. Comparing
this algorithm to a version converted to MATLAB, gave similar solving speed. We concluded that the bottleneck
in the algorithm, solving matrix equations, could not be easily improved, because MATLAB solvers
perform quite optimal: a significantly faster CFD solver implementing the SIMPLE algorithm probably does
not exist.
Nevertheless, when trying to use the obtained results on water and air, we could not obtain any real time
solutions. It is up to further research to define restictions that can guarantee real time solutions for fluids.