Simulating waves in the solar atmosphere with MHD

This thesis deals with modeling and understanding the processes in the atmosphere of the Sun. The solar climate has implications for Earth's climate and effects many of our systems, for instance the Global Positioning System (GPS). In order to protect those systems, it is of vital importance to understand the underlying processes in the Sun that lead to solar flares and mass ejections. We model the energy transfer in the solar atmosphere by the equations of magneto-hydrodynamics (MHD) with a gravity source term, and suitable initial and boundary conditions. Since there are no analytical results available for these equations, one has to simulate the solutions. This is usually done by finite volume methods. Standard finite volume methods, however, are found to be unstable in multiple dimensions in space. Therefore, we developed and implemented new robust and accurate methods for simulating the solutions of the MHD equations. We apply those methods to a two dimensional model with appropriate steady states. The key to obtain a robust scheme is to balance the gravity source with the numerical flux and the Godunov-Powell source. The resulting high order well-balanced schemes are tested on realistic configurations and are found to resolve the complex physical phenomena quite well

Efficient encoding of the weighted MAX k-CUT on a quantum computer using QAOA

The weighted MAX k-CUT problem consists of finding a k-partition of a given weighted undirected graph G(V,E) such that the sum of the weights of the crossing edges is maximized. The problem is of particular interest as it has a multitude of practical applications. We present a formulation of the weighted MAX k-CUT suitable for running the quantum approximate optimization algorithm (QAOA) on noisy intermediate scale quantum (NISQ)-devices to get approximate solutions. The new formulation uses a binary encoding that requires only |V|log_2(k) qubits. The contributions of this paper are as follows: i) A novel decomposition of the phase separation operator based on the binary encoding into basis gates is provided for the MAX k-CUT problem for k >2. ii) Numerical simulations on a suite of test cases comparing different encodings are performed. iii) An analysis of the resources (number of qubits, CX gates) of the different encodings is presented. iv) Formulations and simulations are extended to the case of weighted graphs. For small k and with further improvements when k is not a power of two, our algorithm is a possible candidate to show quantum advantage on NISQ devices.Comment: 14 page

Approximate Riemann Solvers and Robust High-Order Finite Volume Schemes for Multi-Dimensional Ideal MHD Equations

We design stable and high-order accurate finite volume schemes for the ideal MHD equations in multi-dimensions. We obtain excellent numerical stability due to some new elements in the algorithm. The schemes are based on three- and five-wave approximate Riemann solvers of the HLL-type, with the novelty that we allow a varying normal magnetic field. This is achieved by considering the semi-conservative Godunov-Powell form of the MHD equations. We show that it is important to discretize the Godunov-Powell source term in the right way, and that the HLL-type solvers naturally provide a stable upwind discretization. Second-order versions of the ENO- and WENO-type reconstructions are proposed, together with precise modifications necessary to preserve positive pressure and density. Extending the discrete source term to second order while maintaining stability requires non-standard techniques, which we present. The first- and second-order schemes are tested on a suite of numerical experiments demonstrating impressive numerical resolution as well as stability, even on very fine meshe

Constraint Preserving Mixers for the Quantum Approximate Optimization Algorithm

The quantum approximate optimization algorithm/quantum alternating operator ansatz (QAOA) is a heuristic to find approximate solutions of combinatorial optimization problems. Most of the literature is limited to quadratic problems without constraints. However, many practically relevant optimization problems do have (hard) constraints that need to be fulfilled. In this article, we present a framework for constructing mixing operators that restrict the evolution to a subspace of the full Hilbert space given by these constraints. We generalize the “XY”-mixer designed to preserve the subspace of “one-hot” states to the general case of subspaces given by a number of computational basis states. We expose the underlying mathematical structure which reveals more of how mixers work and how one can minimize their cost in terms of the number of CX gates, particularly when Trotterization is taken into account. Our analysis also leads to valid Trotterizations for an “XY”-mixer with fewer CX gates than is known to date. In view of practical implementations, we also describe algorithms for efficient decomposition into basis gates. Several examples of more general cases are presented and analyzed.publishedVersio

Beiträge zur Geschichte des Landkreises Regensburg 25

Heimat gestern und heute II - 41 Arbeitsproben von 12 Autore

Beiträge zur Geschichte des Landkreises Regensburg 22

Heimat gestern und heute I - 48 Arbeitsproben von 14 Autore

Production of He-4 and (4) in Pb-Pb collisions at root(NN)-N-S=2.76 TeV at the LHC

Results on the production of He-4 and (4) nuclei in Pb-Pb collisions at root(NN)-N-S = 2.76 TeV in the rapidity range vertical bar y vertical bar <1, using the ALICE detector, are presented in this paper. The rapidity densities corresponding to 0-10% central events are found to be dN/dy4(He) = (0.8 +/- 0.4 (stat) +/- 0.3 (syst)) x 10(-6) and dN/dy4 = (1.1 +/- 0.4 (stat) +/- 0.2 (syst)) x 10(-6), respectively. This is in agreement with the statistical thermal model expectation assuming the same chemical freeze-out temperature (T-chem = 156 MeV) as for light hadrons. The measured ratio of (4)/He-4 is 1.4 +/- 0.8 (stat) +/- 0.5 (syst). (C) 2018 Published by Elsevier B.V.Peer reviewe

Azimuthal anisotropy of charged jet production in root s(NN)=2.76 TeV Pb-Pb collisions

We present measurements of the azimuthal dependence of charged jet production in central and semi-central root s(NN) = 2.76 TeV Pb-Pb collisions with respect to the second harmonic event plane, quantified as nu(ch)(2) (jet). Jet finding is performed employing the anti-k(T) algorithm with a resolution parameter R = 0.2 using charged tracks from the ALICE tracking system. The contribution of the azimuthal anisotropy of the underlying event is taken into account event-by-event. The remaining (statistical) region-to-region fluctuations are removed on an ensemble basis by unfolding the jet spectra for different event plane orientations independently. Significant non-zero nu(ch)(2) (jet) is observed in semi-central collisions (30-50% centrality) for 20 <p(T)(ch) (jet) <90 GeV/c. The azimuthal dependence of the charged jet production is similar to the dependence observed for jets comprising both charged and neutral fragments, and compatible with measurements of the nu(2) of single charged particles at high p(T). Good agreement between the data and predictions from JEWEL, an event generator simulating parton shower evolution in the presence of a dense QCD medium, is found in semi-central collisions. (C) 2015 CERN for the benefit of the ALICE Collaboration. Published by Elsevier B.V. This is an open access article under the CC BY license (http://creativecommons.org/licenses/by/4.0/).Peer reviewe

Pseudorapidity and transverse-momentum distributions of charged particles in proton-proton collisions at root s=13 TeV

The pseudorapidity (eta) and transverse-momentum (p(T)) distributions of charged particles produced in proton-proton collisions are measured at the centre-of-mass energy root s = 13 TeV. The pseudorapidity distribution in vertical bar eta vertical bar <1.8 is reported for inelastic events and for events with at least one charged particle in vertical bar eta vertical bar <1. The pseudorapidity density of charged particles produced in the pseudorapidity region vertical bar eta vertical bar <0.5 is 5.31 +/- 0.18 and 6.46 +/- 0.19 for the two event classes, respectively. The transverse-momentum distribution of charged particles is measured in the range 0.15 <p(T) <20 GeV/c and vertical bar eta vertical bar <0.8 for events with at least one charged particle in vertical bar eta vertical bar <1. The evolution of the transverse momentum spectra of charged particles is also investigated as a function of event multiplicity. The results are compared with calculations from PYTHIA and EPOS Monte Carlo generators. (C) 2015 CERN for the benefit of the ALICE Collaboration. Published by Elsevier B.V. This is an open access article under the CC BY license (http://creativecommons.org/licenses/by/4.0/).Peer reviewe