Researcher profile

A. Marek

A. Marek contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
8works
0followers
9topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

8 published item(s)

preprint2018arXiv

Optimizations of the Eigensolvers in the ELPA Library

The solution of (generalized) eigenvalue problems for symmetric or Hermitian matrices is a common subtask of many numerical calculations in electronic structure theory or materials science. Solving the eigenvalue problem can easily amount to a sizeable fraction of the whole numerical calculation. For researchers in the field of computational materials science, an efficient and scalable solution of the eigenvalue problem is thus of major importance. The ELPA-library is a well-established dense direct eigenvalue solver library, which has proven to be very efficient and scalable up to very large core counts. In this paper, we describe the latest optimizations of the ELPA-library for new HPC architectures of the Intel Skylake processor family with an AVX-512 SIMD instruction set, or for HPC systems accelerated with recent GPUs. We also describe a complete redesign of the API in a modern modular way, which, apart from a much simpler and more flexible usability, leads to a new path to access system-specific performance optimizations. In order to ensure optimal performance for a particular scientific setting or a specific HPC system, the new API allows the user to influence in straightforward way the internal details of the algorithms and of performance-critical parameters used in the ELPA-library. On top of that, we introduced an autotuning functionality, which allows for finding the best settings in a self-contained automated way. In situations where many eigenvalue problems with similar settings have to be solved consecutively, the autotuning process of the ELPA-library can be done "on-the-fly". Practical applications from materials science which rely on so-called self-consistency iterations can profit from the autotuning. On some examples of scientific interest, simulated with the FHI-aims application, the advantages of the latest optimizations of the ELPA-library are demonstrated.

preprint2013arXiv

SASI Activity in Three-Dimensional Neutrino-Hydrodynamics Simulations of Supernova Cores

The relevance of the standing accretion shock instability (SASI) compared to neutrino-driven convection in three-dimensional (3D) supernova-core environments is still highly controversial. Studying a 27 Msun progenitor, we demonstrate, for the first time, that violent SASI activity can develop in 3D simulations with detailed neutrino transport despite the presence of convection. This result was obtained with the Prometheus-Vertex code with the same sophisticated neutrino treatment so far used only in 1D and 2D models. While buoyant plumes initially determine the nonradial mass motions in the postshock layer, bipolar shock sloshing with growing amplitude sets in during a phase of shock retraction and turns into a violent spiral mode whose growth is only quenched when the infall of the Si/SiO interface leads to strong shock expansion in response to a dramatic decrease of the mass accretion rate. In the phase of large-amplitude SASI sloshing and spiral motions, the postshock layer exhibits nonradial deformation dominated by the lowest-order spherical harmonics (l=1, m=0,-1,+1) in distinct contrast to the higher multipole structures associated with neutrino-driven convection. We find that the SASI amplitudes, shock asymmetry, and nonradial kinetic energy in 3D can exceed those of the corresponding 2D case during extended periods of the evolution. We also perform parametrized 3D simulations of a 25 Msun progenitor, using a simplified, gray neutrino transport scheme, an axis-free Yin-Yang grid, and different amplitudes of random seed perturbations. They confirm the importance of the SASI for another progenitor, its independence of the choice of spherical grid, and its preferred growth for fast accretion flows connected to small shock radii and compact proto-neutron stars as previously found in 2D setups.

preprint2012arXiv

A New Multi-Dimensional General Relativistic Neutrino Hydrodynamics Code for Core-Collapse Supernovae II. Relativistic Explosion Models of Core-Collapse Supernovae

We present the first two-dimensional general relativistic (GR) simulations of stellar core collapse and explosion with the CoCoNuT hydrodynamics code in combination with the VERTEX solver for energy-dependent, three-flavor neutrino transport, using the extended conformal flatness condition for approximating the spacetime metric and a ray-by-ray-plus ansatz to tackle the multi-dimensionality of the transport. For both of the investigated 11.2 and 15 solar mass progenitors we obtain successful, though seemingly marginal, neutrino-driven supernova explosions. This outcome and the time evolution of the models basically agree with results previously obtained with the PROMETHEUS hydro solver including an approximative treatment of relativistic effects by a modified Newtonian potential. However, GR models exhibit subtle differences in the neutrinospheric conditions compared to Newtonian and pseudo-Newtonian simulations. These differences lead to significantly higher luminosities and mean energies of the radiated electron neutrinos and antineutrinos and therefore to larger energy-deposition rates and heating efficiencies in the gain layer with favorable consequences for strong non-radial mass motions and ultimately for an explosion. Moreover, energy transfer to the stellar medium around the neutrinospheres through nucleon recoil in scattering reactions of heavy-lepton neutrinos also enhances the mentioned effects. Together with previous pseudo-Newtonian models the presented relativistic calculations suggest that the treatment of gravity and energy-exchanging neutrino interactions can make differences of even 50-100% in some quantities and is likely to contribute to a finally successful explosion mechanism on no minor level than hydrodynamical differences between different dimensions.

preprint2011arXiv

Neutrino Signal of Electron-Capture Supernovae from Core Collapse to Cooling

An 8.8 solar mass electron-capture supernova (SN) was simulated in spherical symmetry consistently from collapse through explosion to nearly complete deleptonization of the forming neutron star. The evolution time of about 9 s is short because of nucleon-nucleon correlations in the neutrino opacities. After a brief phase of accretion-enhanced luminosities (~200 ms), luminosity equipartition among all species becomes almost perfect and the spectra of electron antineutrinos and muon/tau antineutrinos very similar. We discuss consequences for the neutrino-driven wind as a nucleosynthesis site and for flavor oscillations of SN neutrinos.

preprint2009arXiv

Axisymmetric General Relativistic Simulations of the Accretion-Induced Collapse of White Dwarfs

(Abridged.) The accretion-induced collapse (AIC) of a white dwarf (WD) may lead to the formation of a protoneutron star and a collapse-driven supernova explosion. This process represents a path alternative to thermonuclear disruption of accreting white dwarfs in Type Ia supernovae. Neutrino and gravitational-wave (GW) observations may provide crucial information necessary to reveal a potential AIC. Motivated by the need for systematic predictions of the GW signature of AIC, we present results from an extensive set of general-relativistic AIC simulations using a microphysical finite-temperature equation of state and an approximate treatment of deleptonization during collapse. Investigating a set of 114 progenitor models in rotational equilibrium, with a wide range of rotational configurations, temperatures and central densities, we extend previous Newtonian studies and find that the GW signal has a generic shape akin to what is known as a "Type III" signal in the literature. We discuss the detectability of the emitted GWs, showing that the signal-to-noise ratio for current or next-generation interferometer detectors could be high enough to detect such events in our Galaxy. Some of our AIC models form massive quasi-Keplerian accretion disks after bounce. In rapidly differentially rotating models, the disk mass can be as large as ~0.8-Msun. Slowly and/or uniformly rotating models produce much smaller disks. Finally, we find that the postbounce cores of rapidly spinning white dwarfs can reach sufficiently rapid rotation to develop a nonaxisymmetric rotational instability.

preprint2007arXiv

Effects of Inelastic Neutrino-Nucleus Scattering on Supernova Dynamics and Radiated Neutrino Spectra

Based on the shell model for Gamow-Teller and the Random Phase Approximation for forbidden transitions, we have calculated reaction rates for inelastic neutrino-nucleus scattering (INNS) under supernova (SN) conditions, assuming a matter composition given by Nuclear Statistical Equilibrium. The rates have been incorporated into state-of-the-art stellar core-collapse simulations with detailed energy-dependent neutrino transport. While no significant effect on the SN dynamics is observed, INNS increases the neutrino opacities noticeably and strongly reduces the high-energy tail of the neutrino spectrum emitted in the neutrino burst at shock breakout. Relatedly the expected event rates for the observation of such neutrinos by earthbound detectors are reduced by up to about 60%.

preprint2005arXiv

Exploiting the neutronization burst of a galactic supernova

One of the robust features found in simulations of core-collapse supernovae (SNe) is the prompt neutronization burst, i.e. the first $\sim 25$ milliseconds after bounce when the SN emits with very high luminosity mainly $ν_e$ neutrinos. We examine the dependence of this burst on variations in the input of current SN models and find that recent improvements of the electron capture rates as well as uncertainties in the nuclear equation of state or a variation of the progenitor mass have only little effect on the signature of the neutronization peak in a megaton water Cherenkov detector for different neutrino mixing schemes. We show that exploiting the time-structure of the neutronization peak allows one to identify the case of a normal mass hierarchy and large 13-mixing angle $θ_{13}$, where the peak is absent. The robustness of the predicted total event number in the neutronization burst makes a measurement of the distance to the SN feasible with a precision of about 5%, even in the likely case that the SN is optically obscured.