Showing 1–50 of 116 results for author: Henning, P

Convergence of a Riemannian gradient method for the Gross-Pitaevskii energy functional in a rotating frame

Authors: Patrick Henning, Mahima Yadav

Abstract: This paper investigates the numerical approximation of ground states of rotating Bose-Einstein condensates. This problem requires the minimization of the Gross-Pitaevskii energy $E$ on a Riemannian manifold $\mathbb{S}$. To find a corresponding minimizer $u$, we use a generalized Riemannian gradient method that is based on the concept of Sobolev gradients in combination with an adaptively changing… ▽ More This paper investigates the numerical approximation of ground states of rotating Bose-Einstein condensates. This problem requires the minimization of the Gross-Pitaevskii energy $E$ on a Riemannian manifold $\mathbb{S}$. To find a corresponding minimizer $u$, we use a generalized Riemannian gradient method that is based on the concept of Sobolev gradients in combination with an adaptively changing metric on the manifold. By a suitable choice of the metric, global energy dissipation for the arising gradient method can be proved. The energy dissipation property in turn implies global convergence to the density $|u|^2$ of a critical point $u$ of $E$ on $\mathbb{S}$. Furthermore, we present a precise characterization of the local convergence rates in a neighborhood of each ground state $u$ and how these rates depend on the first spectral gap of $E^{\prime\prime}(u)$ restricted to the $L^2$-orthogonal complement of $u$. With this we establish the first convergence results for a Riemannian gradient method to minimize the Gross-Pitaevskii energy functional in a rotating frame. At the same, we refine previous results obtained in the case without rotation. The major complication in our new analysis is the missing isolation of minimizers, which are at most unique up to complex phase shifts. For that, we introduce an auxiliary iteration in the tangent space $T_{\mathrm{i} u} \mathbb{S}$ and apply the Ostrowski theorem to characterize the asymptotic convergence rates through a weighted eigenvalue problem. Afterwards, we link the auxiliary iteration to the original Riemannian gradient method and bound the spectrum of the weighted eigenvalue problem to obtain quantitative convergence rates. Our findings are validated in numerical experiments. △ Less

Submitted 6 June, 2024; originally announced June 2024.

MSC Class: 35Q55; 65N12; 65N25; 81Q05

Vortex-capturing multiscale spaces for the Ginzburg-Landau equation

Authors: Maria Blum, Christian Döding, Patrick Henning

Abstract: This paper considers minimizers of the Ginzburg-Landau energy functional in particular multiscale spaces which are based on finite elements. The spaces are constructed by localized orthogonal decomposition techniques and their usage for solving the Ginzburg-Landau equation was first suggested in [Dörich, Henning, SINUM 2024]. In this work we further explore their approximation properties and give… ▽ More This paper considers minimizers of the Ginzburg-Landau energy functional in particular multiscale spaces which are based on finite elements. The spaces are constructed by localized orthogonal decomposition techniques and their usage for solving the Ginzburg-Landau equation was first suggested in [Dörich, Henning, SINUM 2024]. In this work we further explore their approximation properties and give an analytical explanation for why vortex structures of energy minimizers can be captured more accurately in these spaces. We quantify the necessary mesh resolution in terms of the Ginzburg-Landau parameter $κ$ and a stabilization parameter $β\ge 0$ that is used in the construction of the multiscale spaces. Furthermore, we analyze how $κ$ affects the necessary locality of the multiscale basis functions and we prove that the choice $β=0$ yields typically the highest accuracy. Our findings are supported by numerical experiments. △ Less

Submitted 27 May, 2024; v1 submitted 23 May, 2024; originally announced May 2024.

MSC Class: 65N12; 65N15; 65N30; 35Q56

MHONGOOSE -- A MeerKAT Nearby Galaxy HI Survey

Authors: W. J. G. de Blok, J. Healy, F. M. Maccagni, D. J. Pisano, A. Bosma, J. English, T. Jarrett, A. Marasco, G. R. Meurer, S. Veronese, F. Bigiel, L. Chemin, F. Fraternali, B. W. Holwerda, P. Kamphuis, H. R. Klöckner, D. Kleiner, A. K. Leroy, M. Mogotsi, K. A. Oman, E. Schinnerer, L. Verdes-Montenegro, T. Westmeier, O. I. Wong, N. Zabel , et al. (35 additional authors not shown)

Abstract: The MHONGOOSE (MeerKAT HI Observations of Nearby Galactic Objects: Observing Southern Emitters) survey maps the distribution and kinematics of the neutral atomic hydrogen (HI) gas in and around 30 nearby star-forming spiral and dwarf galaxies to extremely low HI column densities. The HI column density sensitivity (3 sigma over 16 km/s) ranges from ~ 5 x 10^{17} cm^{-2} at 90'' resolution to ~4 x 1… ▽ More The MHONGOOSE (MeerKAT HI Observations of Nearby Galactic Objects: Observing Southern Emitters) survey maps the distribution and kinematics of the neutral atomic hydrogen (HI) gas in and around 30 nearby star-forming spiral and dwarf galaxies to extremely low HI column densities. The HI column density sensitivity (3 sigma over 16 km/s) ranges from ~ 5 x 10^{17} cm^{-2} at 90'' resolution to ~4 x 10^{19} cm^{-2} at the highest resolution of 7''. The HI mass sensitivity (3 sigma over 50 km/s) is ~5.5 X 10^5 M_sun at a distance of 10 Mpc (the median distance of the sample galaxies). The velocity resolution of the data is 1.4 km/s. One of the main science goals of the survey is the detection of cold, accreting gas in the outskirts of the sample galaxies. The sample was selected to cover a range in HI masses, from 10^7 M_sun to almost 10^{11} M_sun, to optimally sample possible accretion scenarios and environments. The distance to the sample galaxies ranges from 3 to 23 Mpc. In this paper, we present the sample selection, survey design, and observation and reduction procedures. We compare the integrated HI fluxes based on the MeerKAT data with those derived from single-dish measurement and find good agreement, indicating that our MeerKAT observations are recovering all flux. We present HI moment maps of the entire sample based on the first ten percent of the survey data, and find that a comparison of the zeroth- and second-moment values shows a clear separation between the physical properties of the HI in areas with star formation and areas without, related to the formation of a cold neutral medium. Finally, we give an overview of the HI-detected companion and satellite galaxies in the 30 fields, five of which have not previously been catalogued. We find a clear relation between the number of companion galaxies and the mass of the main target galaxy. △ Less

Submitted 6 June, 2024; v1 submitted 2 April, 2024; originally announced April 2024.

Comments: Accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 688, A109 (2024)

Dynamics of electronic phase separation at the laser-induced insulator-metal transition in (La$_{0.6}$Pr$_{0.4}$)$_{0.7}$Ca$_{0.3}$MnO$_3$

Authors: Tim Titze, Maximilian Staabs, Pia Henning, Karen Stroh, Stefan Mathias, Vasily Moshnyaga, Daniel Steil

Abstract: Ultrafast optical excitations allow creating new metastable and hidden states in quantum materials. However, the fundamental material properties required to support new emergent order are largely unknown. Here we show for two colossal magnetoresistive (CMR) manganites that electronic phase separation (EPS) strongly favors non-thermal behavior and exploit this to stabilize an optically-induced cond… ▽ More Ultrafast optical excitations allow creating new metastable and hidden states in quantum materials. However, the fundamental material properties required to support new emergent order are largely unknown. Here we show for two colossal magnetoresistive (CMR) manganites that electronic phase separation (EPS) strongly favors non-thermal behavior and exploit this to stabilize an optically-induced conducting state. Our results shed light on the role of EPS in optical control of CMR manganites and provide guidance for the design of materials that can exhibit non-equilibrium states of matter. △ Less

Submitted 8 December, 2023; originally announced December 2023.

Comments: Main Text: 6 pages, 4 figures Supplemental: 5 pages, 6 figures

Error bounds for discrete minimizers of the Ginzburg-Landau energy in the high-$κ$ regime

Authors: Benjamin Dörich, Patrick Henning

Abstract: In this work, we study discrete minimizers of the Ginzburg-Landau energy in finite element spaces. Special focus is given to the influence of the Ginzburg-Landau parameter $κ$. This parameter is of physical interest as large values can trigger the appearance of vortex lattices. Since the vortices have to be resolved on sufficiently fine computational meshes, it is important to translate the size o… ▽ More In this work, we study discrete minimizers of the Ginzburg-Landau energy in finite element spaces. Special focus is given to the influence of the Ginzburg-Landau parameter $κ$. This parameter is of physical interest as large values can trigger the appearance of vortex lattices. Since the vortices have to be resolved on sufficiently fine computational meshes, it is important to translate the size of $κ$ into a mesh resolution condition, which can be done through error estimates that are explicit with respect to $κ$ and the spatial mesh width $h$. For that, we first work in an abstract framework for a general class of discrete spaces, where we present convergence results in a problem-adapted $κ$-weighted norm. Afterwards we apply our findings to Lagrangian finite elements and a particular generalized finite element construction. In numerical experiments we further explore the asymptotic optimality of our derived $L^2$- and $H^1$-error estimates with respect to $κ$ and $h$. Preasymptotic effects are observed for large mesh sizes $h$. △ Less

Submitted 19 February, 2024; v1 submitted 24 March, 2023; originally announced March 2023.

MSC Class: 65N12; 65N15; 65N30; 35Q56

On discrete ground states of rotating Bose-Einstein condensates

Authors: Patrick Henning, Mahima Yadav

Abstract: The ground states of Bose-Einstein condensates in a rotating frame can be described as constrained minimizers of the Gross-Pitaevskii energy functional with an angular momentum term. In this paper we consider the corresponding discrete minimization problem in Lagrange finite element spaces of arbitrary polynomial order and we investigate the approximation properties of discrete ground states. In p… ▽ More The ground states of Bose-Einstein condensates in a rotating frame can be described as constrained minimizers of the Gross-Pitaevskii energy functional with an angular momentum term. In this paper we consider the corresponding discrete minimization problem in Lagrange finite element spaces of arbitrary polynomial order and we investigate the approximation properties of discrete ground states. In particular, we prove a priori error estimates of optimal order in the $L^2$- and $H^1$-norm, as well as for the ground state energy and the corresponding chemical potential. A central issue in the analysis of the problem is the missing uniqueness of ground states, which is mainly caused by the invariance of the energy functional under complex phase shifts. Our error analysis is therefore based on an Euler-Lagrange functional that we restrict to certain tangent spaces in which we have local uniqueness of ground states. This gives rise to an error decomposition that is ultimately used to derive the desired a priori error estimates. We also present numerical experiments to illustrate various aspects of the problem structure. △ Less

Submitted 25 March, 2024; v1 submitted 1 March, 2023; originally announced March 2023.

MSC Class: 35Q55; 65N15; 65N30

A two level approach for simulating Bose-Einstein condensates by localized orthogonal decomposition

Authors: Christian Döding, Patrick Henning, Johan Wärnegård

Abstract: In this work, we consider the numerical computation of ground states and dynamics of single-component Bose-Einstein condensates (BECs). The corresponding models are spatially discretized with a multiscale finite element approach known as Localized Orthogonal Decomposition (LOD). Despite the outstanding approximation properties of such a discretization in the context of BECs, taking full advantage… ▽ More In this work, we consider the numerical computation of ground states and dynamics of single-component Bose-Einstein condensates (BECs). The corresponding models are spatially discretized with a multiscale finite element approach known as Localized Orthogonal Decomposition (LOD). Despite the outstanding approximation properties of such a discretization in the context of BECs, taking full advantage of it without creating severe computational bottlenecks can be tricky. In this paper, we therefore present two fully-discrete numerical approaches that are formulated in such a way that they take special account of the structure of the LOD spaces. One approach is devoted to the computation of ground states and another one for the computation of dynamics. A central focus of this paper is also the discussion of implementation aspects that are very important for the practical realization of the methods. In particular, we discuss the use of suitable data structures that keep the memory costs economical. The paper concludes with various numerical experiments in 1d, 2d and 3d that investigate convergence rates and approximation properties of the methods and which demonstrate their performance and computational efficiency, also in comparison to spectral and standard finite element approaches. △ Less

Submitted 23 May, 2024; v1 submitted 14 December, 2022; originally announced December 2022.

MSC Class: 35Q55; 65M60; 65M15; 81Q05

Uniform $L^\infty$-bounds for energy-conserving higher-order time integrators for the Gross-Pitaevskii equation with rotation

Authors: Christian Döding, Patrick Henning

Abstract: In this paper, we consider an energy-conserving continuous Galerkin discretization of the Gross-Pitaevskii equation with a magnetic trapping potential and a stirring potential for angular momentum rotation. The discretization is based on finite elements in space and time and allows for arbitrary polynomial orders. It was first analyzed in [O. Karakashian, C. Makridakis; SIAM J. Numer. Anal. 36(6):… ▽ More In this paper, we consider an energy-conserving continuous Galerkin discretization of the Gross-Pitaevskii equation with a magnetic trapping potential and a stirring potential for angular momentum rotation. The discretization is based on finite elements in space and time and allows for arbitrary polynomial orders. It was first analyzed in [O. Karakashian, C. Makridakis; SIAM J. Numer. Anal. 36(6):1779-1807, 1999] in the absence of potential terms and corresponding a priori error estimates were derived in 2D. In this work we revisit the approach in the generalized setting of the Gross-Pitaevskii equation with rotation and we prove uniform $L^\infty$-bounds for the corresponding numerical approximations in 2D and 3D without coupling conditions between the spatial mesh size and the time step size. With this result at hand, we are in particular able to extend the previous error estimates to the 3D setting while avoiding artificial CFL conditions. △ Less

Submitted 21 September, 2023; v1 submitted 4 October, 2022; originally announced October 2022.

MSC Class: 35Q55; 65M60; 65M15; 81Q05

LADUMA: Discovery of a luminous OH megamaser at $z > 0.5$

Authors: Marcin Glowacki, Jordan D. Collier, Amir Kazemi-Moridani, Bradley Frank, Hayley Roberts, Jeremy Darling, Hans-Rainer Klöckner, Nathan Adams, Andrew J. Baker, Matthew Bershady, Tariq Blecher, Sarah-Louise Blyth, Rebecca Bowler, Barbara Catinella, Laurent Chemin, Steven M. Crawford, Catherine Cress, Romeel Davé, Roger Deane, Erwin de Blok, Jacinta Delhaize, Kenneth Duncan, Ed Elson, Sean February, Eric Gawiser , et al. (43 additional authors not shown)

Abstract: In the local Universe, OH megamasers (OHMs) are detected almost exclusively in infrared-luminous galaxies, with a prevalence that increases with IR luminosity, suggesting that they trace gas-rich galaxy mergers. Given the proximity of the rest frequencies of OH and the hyperfine transition of neutral atomic hydrogen (HI), radio surveys to probe the cosmic evolution of HI in galaxies also offer exc… ▽ More In the local Universe, OH megamasers (OHMs) are detected almost exclusively in infrared-luminous galaxies, with a prevalence that increases with IR luminosity, suggesting that they trace gas-rich galaxy mergers. Given the proximity of the rest frequencies of OH and the hyperfine transition of neutral atomic hydrogen (HI), radio surveys to probe the cosmic evolution of HI in galaxies also offer exciting prospects for exploiting OHMs to probe the cosmic history of gas-rich mergers. Using observations for the Looking At the Distant Universe with the MeerKAT Array (LADUMA) deep HI survey, we report the first untargeted detection of an OHM at $z > 0.5$, LADUMA J033046.20$-$275518.1 (nicknamed "Nkalakatha"). The host system, WISEA J033046.26$-$275518.3, is an infrared-luminous radio galaxy whose optical redshift $z \approx 0.52$ confirms the MeerKAT emission line detection as OH at a redshift $z_{\rm OH} = 0.5225 \pm 0.0001$ rather than HI at lower redshift. The detected spectral line has 18.4$σ$ peak significance, a width of $459 \pm 59\,{\rm km\,s^{-1}}$, and an integrated luminosity of $(6.31 \pm 0.18\,{\rm [statistical]}\,\pm 0.31\,{\rm [systematic]}) \times 10^3\,L_\odot$, placing it among the most luminous OHMs known. The galaxy's far-infrared luminosity $L_{\rm FIR} = (1.576 \pm 0.013) \times 10^{12}\,L_\odot$ marks it as an ultra-luminous infrared galaxy; its ratio of OH and infrared luminosities is similar to those for lower-redshift OHMs. A comparison between optical and OH redshifts offers a slight indication of an OH outflow. This detection represents the first step towards a systematic exploitation of OHMs as a tracer of galaxy growth at high redshifts. △ Less

Submitted 5 April, 2022; originally announced April 2022.

Comments: 10 pages, 4 figures. Accepted to ApJ Letters

The dependency of spectral gaps on the convergence of the inverse iteration for a nonlinear eigenvector problem

Authors: Patrick Henning

Abstract: In this paper we consider the generalized inverse iteration for computing ground states of the Gross-Pitaevskii eigenvector problem (GPE). For that we prove explicit linear convergence rates that depend on the maximum eigenvalue in magnitude of a weighted linear eigenvalue problem. Furthermore, we show that this eigenvalue can be bounded by the first spectral gap of a linearized Gross-Pitaevskii o… ▽ More In this paper we consider the generalized inverse iteration for computing ground states of the Gross-Pitaevskii eigenvector problem (GPE). For that we prove explicit linear convergence rates that depend on the maximum eigenvalue in magnitude of a weighted linear eigenvalue problem. Furthermore, we show that this eigenvalue can be bounded by the first spectral gap of a linearized Gross-Pitaevskii operator, recovering the same rates as for linear eigenvector problems. With this we establish the first local convergence result for the basic inverse iteration for the GPE without damping. We also show how our findings directly generalize to extended inverse iterations, such as the Gradient Flow Discrete Normalized (GFDN) proposed in [W. Bao, Q. Du, SIAM J. Sci. Comput., 25 (2004)] or the damped inverse iteration suggested in [P. Henning, D. Peterseim, SIAM J. Numer. Anal., 53 (2020)]. Our analysis also reveals why the inverse iteration for the GPE does not react favourably to spectral shifts. This empirical observation can now be explained with a blow-up of a weighting function that crucially contributes to the convergence rates. Our findings are illustrated by numerical experiments. △ Less

Submitted 8 March, 2024; v1 submitted 15 February, 2022; originally announced February 2022.

MSC Class: 35Q55; 65N12; 65N25; 65N30; 81Q05

On optimal convergence rates for discrete minimizers of the Gross-Pitaevskii energy in LOD spaces

Authors: Patrick Henning, Anna Persson

Abstract: In this paper we revisit a two-level discretization based on the Localized Orthogonal Decomposition (LOD). It was originally proposed in [P.Henning, A.Målqvist, D.Peterseim. SIAM J. Numer. Anal.52-4:1525-1550, 2014] to compute ground states of Bose-Einstein condensates by finding discrete minimizers of the Gross-Pitaevskii energy functional. The established convergence rates for the method appeare… ▽ More In this paper we revisit a two-level discretization based on the Localized Orthogonal Decomposition (LOD). It was originally proposed in [P.Henning, A.Målqvist, D.Peterseim. SIAM J. Numer. Anal.52-4:1525-1550, 2014] to compute ground states of Bose-Einstein condensates by finding discrete minimizers of the Gross-Pitaevskii energy functional. The established convergence rates for the method appeared however suboptimal compared to numerical observations and a proof of optimal rates in this setting remained open. In this paper we shall close this gap by proving optimal order error estimates for the $L^2$- and $H^1$-error between the exact ground state and discrete minimizers, as well as error estimates for the ground state energy and the ground state eigenvalue. In particular, the achieved convergence rates for the energy and the eigenvalue are of $6$th order with respect to the mesh size on which the discrete LOD space is based, without making any additional regularity assumptions. These high rates justify the use of very coarse meshes, which significantly reduces the computational effort for finding accurate approximations of ground states. In addition, we include numerical experiments that confirm the optimality of the new theoretical convergence rates, both for smooth and discontinuous potentials. △ Less

Submitted 14 April, 2023; v1 submitted 15 December, 2021; originally announced December 2021.

Gender Systematics in the NRAO Proposal Review System

Authors: Gareth Hunt, Frederic R. Schwab, P. A. Henning, Dana S. Balser

Abstract: Several recent investigations indicate the existence of gender-related systematic trends in the peer review of proposals for observations on astronomical facilities. This includes the National Radio Astronomy Observatory (NRAO) where there is evidence of a gender imbalance in the rank of proposals with male principal investigators (PIs) favored over female PIs. Since semester 2017A (17A), the NRAO… ▽ More Several recent investigations indicate the existence of gender-related systematic trends in the peer review of proposals for observations on astronomical facilities. This includes the National Radio Astronomy Observatory (NRAO) where there is evidence of a gender imbalance in the rank of proposals with male principal investigators (PIs) favored over female PIs. Since semester 2017A (17A), the NRAO has taken the following steps: (1) inform science review panels (SRPs) and the telescope time allocation committee (TAC) about the gender imbalance; and (2) increase the female representation on SRPs and the TAC to reflect the community demographics. Here we analyze SRP normalized rank-ordered scores, or linear ranks, by PI gender for NRAO observing proposals from semesters 12A-21A. We use bootstrap resampling to generate modeled distributions and the Anderson-Darling (AD) test to evaluate the probability that the linear rank distributions for male and female PIs are drawn from the same parent sample. We find that between semesters 12A-17A that male PIs are favored over female PIs (AD p-value 0.0084), whereas between semesters 17B-21A female PIs are favored over male PIs, but at a lower significance (AD p-value 0.11). Therefore the gender imbalance is currently being ameliorated, but this imbalance may have been reversed. Regardless, we plan to adopt a dual-anonymous approach to proposal review to reduce the possibility of bias to occur. △ Less

Submitted 24 October, 2021; originally announced October 2021.

Comments: Accepted for publication in the PASP

Coherent phonon transport and minimum of thermal conductivity in LaMnO$_3$/SrMnO$_3$ superlattices

Authors: Dennis Meyer, Daniel Metternich, Pia Henning, Jan Philipp Bange, Robert Gruhl, Vitaly Bruchmann-Bamberg, Vasily Moshnyaga, Henning Ulrichs

Abstract: In this Letter, we present investigations of thermal transport in (LaMnO$_3$)$_n$/(SrMnO$_3$)$_n$ superlattices (SLs) with SL periods $Λ=2n$ between $2$ and $12$ unit cells. The analysis of the experimental data clearly reveals a minimum in the thermal conductivity $κ$ at $Λ=6$ u.c. Furthermore, we theoretically estimate the phonon mean free path to be about $3Λ$ in our SLs. These findings clearly… ▽ More In this Letter, we present investigations of thermal transport in (LaMnO$_3$)$_n$/(SrMnO$_3$)$_n$ superlattices (SLs) with SL periods $Λ=2n$ between $2$ and $12$ unit cells. The analysis of the experimental data clearly reveals a minimum in the thermal conductivity $κ$ at $Λ=6$ u.c. Furthermore, we theoretically estimate the phonon mean free path to be about $3Λ$ in our SLs. These findings clearly show that, in a wide temperature range of $\pm 100\,$K around room temperature, thermal transport in our manganite SLs is coherent. In addition, we demonstrate that antiferromagnetic magnetic ordering can control the degree of coherency. This opens the challenging perspective for active control of coherent transport in correlated matter. △ Less

Submitted 12 August, 2021; originally announced August 2021.

EB-DEVS: A Formal Framework for Modeling and Simulation of Emergent Behavior in Dynamic Complex Systems

Authors: Daniel J. Foguelman, Philipp Henning, Adelinde Uhrmacher, Rodrigo Castro

Abstract: Emergent behavior is a key feature defining a system under study as a complex system. Simulation has been recognized as the only way to deal with the study of the emergency of properties (at a macroscopic level) among groups of system components (at a microscopic level), for the manifestations of emergent structures cannot be deduced from analysing components in isolation. A systems-oriented gener… ▽ More Emergent behavior is a key feature defining a system under study as a complex system. Simulation has been recognized as the only way to deal with the study of the emergency of properties (at a macroscopic level) among groups of system components (at a microscopic level), for the manifestations of emergent structures cannot be deduced from analysing components in isolation. A systems-oriented generalisation must consider the presence of feedback loops (micro components react to macro properties), interaction among components of different classes (modular composition) and layered interaction of subsystems operating at different spatio-temporal scales (hierarchical organisation). In this work we introduce Emergent Behavior-DEVS (EB-DEVS) a Modeling and Simulation (M&S) formalism that permits reasoning about complex systems where emergent behavior is placed at the forefront of the analysis activity. EB-DEVS builds on the DEVS formalism, adding upward/downward communication channels to well-established capabilities for modular and hierarchical M&S of heterogeneous multi-formalism systems. EB-DEVS takes a minimalist stance on expressiveness, introducing a small set of extensions on Classic DEVS that can cope with emergent behavior, and making both formalisms interoperable (the modeler decides which subsystems deserve to be expressed via micro-macro dynamics). We present three case studies: flocks of birds with learning, population epidemics with vaccination and sub-cellular dynamics with homeostasis, through which we showcase how EB-DEVS performs by placing emergent properties at the center of the M&S process. △ Less

Submitted 13 October, 2020; v1 submitted 10 October, 2020; originally announced October 2020.

Comments: 38 page document with: original content 25 pages, references 3 pages, appendices 10 pages

Journal ref: Journal of Computational Science Volume 53, July 2021, 101387

Superconvergence of time invariants for the Gross-Pitaevskii equation

Authors: Patrick Henning, Johan Wärnegård

Abstract: This paper considers the numerical treatment of the time-dependent Gross-Pitaevskii equation. In order to conserve the time invariants of the equation as accurately as possible, we propose a Crank-Nicolson-type time discretization that is combined with a suitable generalized finite element discretization in space. The space discretization is based on the technique of Localized Orthogonal Decomposi… ▽ More This paper considers the numerical treatment of the time-dependent Gross-Pitaevskii equation. In order to conserve the time invariants of the equation as accurately as possible, we propose a Crank-Nicolson-type time discretization that is combined with a suitable generalized finite element discretization in space. The space discretization is based on the technique of Localized Orthogonal Decompositions (LOD) and allows to capture the time invariants with an accuracy of order $\mathcal{O}(H^6)$ with respect to the chosen mesh size $H$. This accuracy is preserved due to the conservation properties of the time stepping method. Furthermore, we prove that the resulting scheme approximates the exact solution in the $L^{\infty}(L^2)$-norm with order $\mathcal{O}(τ^2 + H^4)$, where $τ$ denotes the step size. The computational efficiency of the method is demonstrated in numerical experiments for a benchmark problem with known exact solution. △ Less

Submitted 19 October, 2021; v1 submitted 17 August, 2020; originally announced August 2020.

MSC Class: 35Q55; 65M60; 65M15; 81Q05

Localization and delocalization of ground states of Bose-Einstein condensates under disorder

Authors: Robert Altmann, Patrick Henning, Daniel Peterseim

Abstract: This paper studies the localization behaviour of Bose-Einstein condensates in disorder potentials, modeled by a Gross-Pitaevskii eigenvalue problem on a bounded interval. In the regime of weak particle interaction, we are able to quantify exponential localization of the ground state, depending on statistical parameters and the strength of the potential. Numerical studies further show delocalizatio… ▽ More This paper studies the localization behaviour of Bose-Einstein condensates in disorder potentials, modeled by a Gross-Pitaevskii eigenvalue problem on a bounded interval. In the regime of weak particle interaction, we are able to quantify exponential localization of the ground state, depending on statistical parameters and the strength of the potential. Numerical studies further show delocalization if we leave the identified parameter range, which is in agreement with experimental data. These mathematical and numerical findings allow the prediction of physically relevant regimes where localization of ground states may be observed experimentally. △ Less

Submitted 8 October, 2021; v1 submitted 1 June, 2020; originally announced June 2020.

Comments: accepted for publication in SIAM J. Appl. Math

MSC Class: 47H40; 81Q10; 65N25

WALLABY -- An SKA Pathfinder HI Survey

Authors: B. S. Koribalski, L. Staveley-Smith, T. Westmeier, P. Serra, K. Spekkens, O. I. Wong, C. D. P. Lagos, D. Obreschkow, E. V. Ryan-Weber, M. Zwaan, V. Kilborn, G. Bekiaris, K. Bekki, F. Bigiel, A. Boselli, A. Bosma, B. Catinella, G. Chauhan, M. E. Cluver, M. Colless, H. M. Courtois, R. A. Crain, W. J. G. de Blok, H. Dénes, A. R. Duffy , et al. (45 additional authors not shown)

Abstract: The Widefield ASKAP L-band Legacy All-sky Blind surveY (WALLABY) is a next-generation survey of neutral hydrogen (HI) in the Local Universe. It uses the widefield, high-resolution capability of the Australian Square Kilometer Array Pathfinder (ASKAP), a radio interferometer consisting of 36 x 12-m dishes equipped with Phased-Array Feeds (PAFs), located in an extremely radio-quiet zone in Western A… ▽ More The Widefield ASKAP L-band Legacy All-sky Blind surveY (WALLABY) is a next-generation survey of neutral hydrogen (HI) in the Local Universe. It uses the widefield, high-resolution capability of the Australian Square Kilometer Array Pathfinder (ASKAP), a radio interferometer consisting of 36 x 12-m dishes equipped with Phased-Array Feeds (PAFs), located in an extremely radio-quiet zone in Western Australia. WALLABY aims to survey three-quarters of the sky (-90 degr < Dec < +30 degr) to a redshift of z < 0.26, and generate spectral line image cubes at ~30 arcsec resolution and ~1.6 mJy/beam per 4 km/s channel sensitivity. ASKAP's instantaneous field of view at 1.4 GHz, delivered by the PAF's 36 beams, is about 30 sq deg. At an integrated signal-to-noise ratio of five, WALLABY is expected to detect over half a million galaxies with a mean redshift of z ~ 0.05 (~200 Mpc). The scientific goals of WALLABY include: (a) a census of gas-rich galaxies in the vicinity of the Local Group; (b) a study of the HI properties of galaxies, groups and clusters, in particular the influence of the environment on galaxy evolution; and (c) the refinement of cosmological parameters using the spatial and redshift distribution of low-bias gas-rich galaxies. For context we provide an overview of previous large-scale HI surveys. Combined with existing and new multi-wavelength sky surveys, WALLABY will enable an exciting new generation of panchromatic studies of the Local Universe. - First results from the WALLABY pilot survey are revealed, with initial data products publicly available in the CSIRO ASKAP Science Data Archive (CASDA). △ Less

Submitted 7 July, 2020; v1 submitted 17 February, 2020; originally announced February 2020.

Comments: Accepted for publication in ApSS (38 pages, 14 figures), see also https://www.atnf.csiro.au/research/WALLABY/ - Contact email: [email protected]

Journal ref: Ap&SS 365, 118 (2020) - https://rdcu.be/b5Bfg

Shadow Lagrangian dynamics for superfluidity

Authors: Patrick Henning, Anders M. N. Niklasson

Abstract: Motivated by a similar approach for Born-Oppenheimer molecular dynamics, this paper proposes an extended "shadow" Lagrangian density for quantum states of superfluids. The extended Lagrangian contains an additional field variable that is forced to follow the wave function of the quantum state through a rapidly oscillating extended harmonic oscillator. By considering the adiabatic limit for large f… ▽ More Motivated by a similar approach for Born-Oppenheimer molecular dynamics, this paper proposes an extended "shadow" Lagrangian density for quantum states of superfluids. The extended Lagrangian contains an additional field variable that is forced to follow the wave function of the quantum state through a rapidly oscillating extended harmonic oscillator. By considering the adiabatic limit for large frequencies of the harmonic oscillator, we can derive the two equations of motions, a Schrödinger-type equation for the quantum state and a wave equation for the extended field variable. The equations are coupled in a nonlinear way, but each equation individually is linear with respect to the variable that it defines. The computational advantage of this new system is that it can be easily discretized using linear time stepping methods, where we propose to use a Crank-Nicolson-type approach for the Schrödinger equation and an extended leapfrog scheme for the wave equation. Furthermore, the difference between the quantum state and the extended field variable defines a consistency error that should go to zero if the frequency tends to infinity. By coupling the time-step size in our discretization to the frequency of the harmonic oscillator we can extract an easily computable consistency error indicator that can be used to estimate the numerical error without any additional costs. The findings are illustrated in numerical experiments. △ Less

Submitted 11 January, 2021; v1 submitted 21 December, 2019; originally announced December 2019.

Computational homogenization of time-harmonic Maxwell's equations

Authors: Patrick Henning, Anna Persson

Abstract: In this paper we consider a numerical homogenization technique for curl-curl-problems that is based on the framework of the Localized Orthogonal Decomposition and which was proposed in [D. Gallistl, P. Henning, B. Verfürth. SIAM J. Numer. Anal. 56-3:1570-1596, 2018] for problems with essential boundary conditions. The findings of the aforementioned work establish quantitative homogenization result… ▽ More In this paper we consider a numerical homogenization technique for curl-curl-problems that is based on the framework of the Localized Orthogonal Decomposition and which was proposed in [D. Gallistl, P. Henning, B. Verfürth. SIAM J. Numer. Anal. 56-3:1570-1596, 2018] for problems with essential boundary conditions. The findings of the aforementioned work establish quantitative homogenization results for the time-harmonic Maxwell's equations that hold beyond assumptions of periodicity, however, a practical realization of the approach was left open. In this paper, we transfer the findings from essential boundary conditions to natural boundary conditions and we demonstrate that the approach yields a computable numerical method. We also investigate how boundary values of the source term can effect the computational complexity and accuracy. Our findings will be supported by various numerical experiments, both in $2D$ and $3D$. △ Less

Submitted 3 March, 2020; v1 submitted 16 October, 2019; originally announced October 2019.

MSC Class: 35Q61; 65N30; 65N12; 78M10

The Arecibo L-band Feed Array Zone of Avoidance (ALFAZOA) Shallow Survey

Authors: Monica Sanchez-Barrantes, Patricia A. Henning, Travis McIntyre, Emmanuel Momjian, Robert Minchin, Jessica L. Rosenberg, Stephen Schneider, Lister Staveley-Smith, Wim van Driel, Mpati Ramatsoku, Zhon Butcher, Ensieh Vaez

Abstract: The Arecibo L-band Feed Array Zone of Avoidance (ALFAZOA) Shallow Survey is a blind HI survey of the extragalactic sky behind the northern Milky Way conducted with the ALFA receiver on the 305m Arecibo Radio Telescope. ALFAZOA Shallow covered 900 square degrees at full sensitivity from 30° ${\leq} l {\leq} $75° and |b|$ {\leq}$ 10° and an additional 460 square degrees at limited sensitivity at lat… ▽ More The Arecibo L-band Feed Array Zone of Avoidance (ALFAZOA) Shallow Survey is a blind HI survey of the extragalactic sky behind the northern Milky Way conducted with the ALFA receiver on the 305m Arecibo Radio Telescope. ALFAZOA Shallow covered 900 square degrees at full sensitivity from 30° ${\leq} l {\leq} $75° and |b|$ {\leq}$ 10° and an additional 460 square degrees at limited sensitivity at latitudes up to 20°. It has an rms sensitivity of 5-7 mJy and a velocity resolution of 9 - 20.6 km s$^{-1}$, and detected 403 galaxies out to a recessional velocity of 12,000 km s$^{-1}$, with an angular resolution of 3.4' and a positional accuracy between 0.2' and 1.7'. The survey is complete above an integrated line flux $F_{HI}$ = 2.0 Jy km s$^{-1}$ for half the survey, and above $F_{HI}$ = 2.8 Jy km s$^{-1}$ for the other half. Forty-three percent of the ALFAZOA HI detections have at least one possible optical/NIR counterpart in the literature, and an additional 16% have counterparts that only included previous HI measurements. There are fewer counterparts in regions of high extinction and for galaxies with lower HI mass. Comparing the results of the survey to the predictions of Erdogdu et al. (2006), and using their nomenclature, ALFAZOA confirms the position and extent in the ZOA of the C7, C$ζ$, Pegasus, Corona Borealis and Delphinus structures, but not of the Cygnus void. Two new structures are identified, both connected to the C7 overdensity; one extends to 35°, and the other crosses the ZOA. △ Less

Submitted 11 October, 2019; originally announced October 2019.

Comments: 19 pages, 18 figures, accepted for publication in the Astronomical Journal

The $J$-method for the Gross-Pitaevskii eigenvalue problem

Authors: Robert Altmann, Patrick Henning, Daniel Peterseim

Abstract: This paper studies the $J$-method of [E. Jarlebring, S. Kvaal, W. Michiels. SIAM J. Sci. Comput. 36-4:A1978-A2001, 2014] for nonlinear eigenvector problems in a general Hilbert space framework. This is the basis for variational discretization techniques and a mesh-independent numerical analysis. A simple modification of the method mimics an energy-decreasing discrete gradient flow. In the case of… ▽ More This paper studies the $J$-method of [E. Jarlebring, S. Kvaal, W. Michiels. SIAM J. Sci. Comput. 36-4:A1978-A2001, 2014] for nonlinear eigenvector problems in a general Hilbert space framework. This is the basis for variational discretization techniques and a mesh-independent numerical analysis. A simple modification of the method mimics an energy-decreasing discrete gradient flow. In the case of the Gross-Pitaevskii eigenvalue problem, we prove global convergence towards an eigenfunction for a damped version of the $J$-method. More importantly, when the iterations are sufficiently close to an eigenfunction, the damping can be switched off and we recover a local linear convergence rate previously known from the discrete setting. This quantitative convergence analysis is closely connected to the~$J$-method's unique feature of sensitivity with respect to spectral shifts. Contrary to classical gradient flows, this allows both the selective approximation of excited states as well as the amplification of convergence beyond linear rates in the spirit of the Rayleigh quotient iteration for linear eigenvalue problems. These advantageous convergence properties are demonstrated in a series of numerical experiments involving exponentially localized states under disorder potentials and vortex lattices in rotating traps. △ Less

Submitted 9 December, 2020; v1 submitted 1 August, 2019; originally announced August 2019.

MSC Class: 35Q55; 65N12; 65N25; 81Q05

A note on optimal $H^1$-error estimates for Crank-Nicolson approximations to the nonlinear Schrödinger equation

Authors: Patrick Henning, Johan Wärnegård

Abstract: In this paper we consider a mass- and energy--conserving Crank-Nicolson time discretization for a general class of nonlinear Schrödinger equations. This scheme, which enjoys popularity in the physics community due to its conservation properties, was already subject to several analytical and numerical studies. However, a proof of optimal $L^{\infty}(H^1)$-error estimates is still open, both in the… ▽ More In this paper we consider a mass- and energy--conserving Crank-Nicolson time discretization for a general class of nonlinear Schrödinger equations. This scheme, which enjoys popularity in the physics community due to its conservation properties, was already subject to several analytical and numerical studies. However, a proof of optimal $L^{\infty}(H^1)$-error estimates is still open, both in the semi-discrete Hilbert space setting, as well as in fully-discrete finite element settings. This paper aims at closing this gap in the literature. We also suggest a fixed point iteration to solve the arising nonlinear system of equations that makes the method easy to implement and efficient. This is illustrated by numerical experiments. △ Less

Submitted 7 July, 2020; v1 submitted 5 July, 2019; originally announced July 2019.

MSC Class: 35Q55; 65M60; 65M15; 81Q05

The FAIR Funder pilot programme to make it easy for funders to require and for grantees to produce FAIR Data

Authors: P. Wittenburg, H. Pergl Sustkova, A. Montesanti, S. M. Bloemers, S. H. de Waard, M. A. Musen, J. B. Graybeal, K. M. Hettne, A. Jacobsen, R. Pergl, R. W. W. Hooft, C. Staiger, C. W. G. van Gelder, S. L. Knijnenburg, A. C. van Arkel, B. Meerman, M. D. Wilkinson, S-A Sansone, P. Rocca-Serra, P. McQuilton, A. N. Gonzalez-Beltran, G. J. C. Aben, P. Henning, S. Alencar, C. Ribeiro , et al. (35 additional authors not shown)

Abstract: There is a growing acknowledgement in the scientific community of the importance of making experimental data machine findable, accessible, interoperable, and reusable (FAIR). Recognizing that high quality metadata are essential to make datasets FAIR, members of the GO FAIR Initiative and the Research Data Alliance (RDA) have initiated a series of workshops to encourage the creation of Metadata for… ▽ More There is a growing acknowledgement in the scientific community of the importance of making experimental data machine findable, accessible, interoperable, and reusable (FAIR). Recognizing that high quality metadata are essential to make datasets FAIR, members of the GO FAIR Initiative and the Research Data Alliance (RDA) have initiated a series of workshops to encourage the creation of Metadata for Machines (M4M), enabling any self-identified stakeholder to define and promote the reuse of standardized, comprehensive machine-actionable metadata. The funders of scientific research recognize that they have an important role to play in ensuring that experimental results are FAIR, and that high quality metadata and careful planning for FAIR data stewardship are central to these goals. We describe the outcome of a recent M4M workshop that has led to a pilot programme involving two national science funders, the Health Research Board of Ireland (HRB) and the Netherlands Organisation for Health Research and Development (ZonMW). These funding organizations will explore new technologies to define at the time that a request for proposals is issued the minimal set of machine-actionable metadata that they would like investigators to use to annotate their datasets, to enable investigators to create such metadata to help make their data FAIR, and to develop data-stewardship plans that ensure that experimental data will be managed appropriately abiding by the FAIR principles. The FAIR Funders design envisions a data-management workflow having seven essential stages, where solution providers are openly invited to participate. The initial pilot programme will launch using existing computer-based tools of those who attended the M4M Workshop. △ Less

Submitted 6 March, 2019; v1 submitted 26 February, 2019; originally announced February 2019.

Comments: This is a pre-print of the FAIR Funders pilot, an outcome of the first Metadata for Machines workshop, see: https://www.go-fair.org/resources/go-fair-workshop-series/metadata-for-machines-workshops/. Corresponding author: E. A Schultes, ORCID 0000-0001-8888-635X

Sobolev gradient flow for the Gross-Pitaevskii eigenvalue problem: global convergence and computational efficiency

Authors: Patrick Henning, Daniel Peterseim

Abstract: We propose a new normalized Sobolev gradient flow for the Gross-Pitaevskii eigenvalue problem based on an energy inner product that depends on time through the density of the flow itself. The gradient flow is well-defined and converges to an eigenfunction. For ground states we can quantify the convergence speed as exponentially fast where the rate depends on spectral gaps of a linearized operator.… ▽ More We propose a new normalized Sobolev gradient flow for the Gross-Pitaevskii eigenvalue problem based on an energy inner product that depends on time through the density of the flow itself. The gradient flow is well-defined and converges to an eigenfunction. For ground states we can quantify the convergence speed as exponentially fast where the rate depends on spectral gaps of a linearized operator. The forward Euler time discretization of the flow yields a numerical method which generalizes the inverse iteration for the nonlinear eigenvalue problem. For sufficiently small time steps, the method reduces the energy in every step and converges globally in $H^1$ to an eigenfunction. In particular, for any nonnegative starting value, the ground state is obtained. A series of numerical experiments demonstrates the computational efficiency of the method and its competitiveness with established discretizations arising from other gradient flows for this problem. △ Less

Submitted 2 April, 2020; v1 submitted 3 December, 2018; originally announced December 2018.

MSC Class: 35Q55; 65N12; 65N25; 65N30; 81Q05

CHILES: HI morphology and galaxy environment at z=0.12 and z=0.17

Authors: Kelley M. Hess, Nicholas M. Luber, Ximena Fernández, Hansung B. Gim, J. H. van Gorkom, Emmanuel Momjian, Julia Gross, Martin Meyer, Attila Popping, Luke J. M. Davies, Lucas Hunt, Kathryn Kreckel, Danielle Lucero, D. J. Pisano, Monica Sanchez-Barrantes, Min S. Yun, Richard Dodson, Kevin Vinsen, Andreas Wicenec, Chen Wu, Matthew A. Bershady, Aeree Chung, Julie D. Davis, Jennifer Donovan Meyer, Patricia Henning , et al. (5 additional authors not shown)

Abstract: We present a study of 16 HI-detected galaxies found in 178 hours of observations from Epoch 1 of the COSMOS HI Large Extragalactic Survey (CHILES). We focus on two redshift ranges between 0.108 <= z <= 0.127 and 0.162 <= z <= 0.183 which are among the worst affected by radio frequency interference (RFI). While this represents only 10% of the total frequency coverage and 18% of the total expected t… ▽ More We present a study of 16 HI-detected galaxies found in 178 hours of observations from Epoch 1 of the COSMOS HI Large Extragalactic Survey (CHILES). We focus on two redshift ranges between 0.108 <= z <= 0.127 and 0.162 <= z <= 0.183 which are among the worst affected by radio frequency interference (RFI). While this represents only 10% of the total frequency coverage and 18% of the total expected time on source compared to what will be the full CHILES survey, we demonstrate that our data reduction pipeline recovers high quality data even in regions severely impacted by RFI. We report on our in-depth testing of an automated spectral line source finder to produce HI total intensity maps which we present side-by-side with significance maps to evaluate the reliability of the morphology recovered by the source finder. We recommend that this become a common place manner of presenting data from upcoming HI surveys of resolved objects. We use the COSMOS 20k group catalogue, and we extract filamentary structure using the topological DisPerSE algorithm to evaluate the \hi\ morphology in the context of both local and large-scale environments and we discuss the shortcomings of both methods. Many of the detections show disturbed HI morphologies suggesting they have undergone a recent interaction which is not evident from deep optical imaging alone. Overall, the sample showcases the broad range of ways in which galaxies interact with their environment. This is a first look at the population of galaxies and their local and large-scale environments observed in HI by CHILES at redshifts beyond the z=0.1 Universe. △ Less

Submitted 28 November, 2018; originally announced November 2018.

Comments: 23 pages, 12 figures, 1 interactive 3D figure, accepted to MNRAS

The Nancay HI Zone of Avoidance survey of 2MASS bright galaxies

Authors: Renée C. Kraan-Korteweg, Wim van Driel, Anja C. Schröder, Mpati Ramatsoku, Patricia A. Henning

Abstract: To complement the 2MASS Redshift Survey (2MRS) and the 2MASS Tully-Fisher survey (2MTF) a search for 21cm HI line emission of 2MASS bright galaxy candidates has been pursued along the dust-obscured plane of the Milky Way with the 100m Nancay Radio Telescope. For our sample selection we adopted an isophotal extinction-corrected K-band magnitude limit of $K_s^o = 11.25$mag, corresponding to the firs… ▽ More To complement the 2MASS Redshift Survey (2MRS) and the 2MASS Tully-Fisher survey (2MTF) a search for 21cm HI line emission of 2MASS bright galaxy candidates has been pursued along the dust-obscured plane of the Milky Way with the 100m Nancay Radio Telescope. For our sample selection we adopted an isophotal extinction-corrected K-band magnitude limit of $K_s^o = 11.25$mag, corresponding to the first 2MRS data release and 2MTF, for which the 2MASX completeness level remains fairly constant deep into the Zone of Avoidance (ZoA). About one thousand galaxies without prior redshift measurement accessible from Nancay (Dec > -40°) were observed to an rms noise level of ~3 mJy for the velocity range -250 to 10'600 km/s. This resulted in 220 clear and 12 marginal detections of the target sample. Only few detections have redshifts above 8000 km/s due to recurring radio frequency interference (RFI). A further 29 detections and 6 marginals have their origin in non-target galaxies in the telescope beam. The newly detected galaxies are on average considerably more \HI-rich (mostly $10^9 - 10^{10}$M$_\odot$) compared to systematic (blind) HI surveys. The HI detections reveal various new filaments crossing the mostly uncharted northern ZoA (e.g. at $\ell \sim 90°, 130°, 160°$), whilst consolidating galaxy agglomerations in Monoceros and Puppis ($\ell \sim 220°, 240°$). Considerably new insight has been gained about the extent of the Perseus-Pisces Supercluster through the confirmation of a ridge ($\ell \sim 160°$) encompassing the 3C129 cluster that links Perseus-Pisces to Lynx, and the continuation of the second Perseus-Pisces arm ($\ell \sim 90°$) across the ZoA. △ Less

Submitted 17 August, 2018; originally announced August 2018.

Comments: 19 pages, 11 figures, 7 tables; accepted for publication in MNRAS; supplementary material available at http://www.acgc.uct.ac.za/~kraan/NRT_ZOA/Supplement.pdf

Numerical comparison of mass-conservative schemes for the Gross-Pitaevskii equation

Authors: Patrick Henning, Johan Wärnegård

Abstract: In this paper we present a numerical comparison of various mass-conservative discretizations for the time-dependent Gross-Pitaevskii equation. We have three main objectives. First, we want to clarify how purely mass-conservative methods perform compared to methods that are additionally energy-conservative or symplectic. Second, we shall compare the accuracy of energy-conservative and symplectic me… ▽ More In this paper we present a numerical comparison of various mass-conservative discretizations for the time-dependent Gross-Pitaevskii equation. We have three main objectives. First, we want to clarify how purely mass-conservative methods perform compared to methods that are additionally energy-conservative or symplectic. Second, we shall compare the accuracy of energy-conservative and symplectic methods among each other. Third, we will investigate if a linearized energy-conserving method suffers from a loss of accuracy compared to an approach which requires to solve a full nonlinear problem in each time-step. In order to obtain a representative comparison, our numerical experiments cover different physically relevant test cases, such as traveling solitons, stationary multi-solitons, Bose-Einstein condensates in an optical lattice and vortex pattern in a rapidly rotating superfluid. We shall also consider a computationally severe test case involving a pseudo Mott insulator. Our space discretization is based on finite elements throughout the paper. We will also give special attention to long time behavior and possible coupling conditions between time-step sizes and mesh sizes. The main observation of this paper is that mass conservation alone will not lead to a competitive method in complex settings. Furthermore, energy-conserving and symplectic methods are both reliable and accurate, yet, the energy-conservative schemes achieve a visibly higher accuracy in our test cases. Finally, the scheme that performs best throughout our experiments is an energy-conserving relaxation scheme with linear time-stepping proposed by C. Besse (SINUM,42(3):934--952,2004). △ Less

Submitted 26 June, 2019; v1 submitted 27 April, 2018; originally announced April 2018.

MSC Class: 35Q55; 65M60; 65Y20; 65Z99; 65P10; 81Q05

Quantitative Anderson localization of Schrödinger eigenstates under disorder potentials

Authors: Robert Altmann, Patrick Henning, Daniel Peterseim

Abstract: This paper concerns spectral properties of linear Schrödinger operators under oscillatory high-amplitude potentials on bounded domains. Depending on the degree of disorder, we prove the existence of spectral gaps amongst the lowermost eigenvalues and the emergence of exponentially localized states. We quantify the rate of decay in terms of geometric parameters that characterize the potential. The… ▽ More This paper concerns spectral properties of linear Schrödinger operators under oscillatory high-amplitude potentials on bounded domains. Depending on the degree of disorder, we prove the existence of spectral gaps amongst the lowermost eigenvalues and the emergence of exponentially localized states. We quantify the rate of decay in terms of geometric parameters that characterize the potential. The proofs are based on the convergence theory of iterative solvers for eigenvalue problems and their optimal local preconditioning by domain decomposition. △ Less

Submitted 9 February, 2020; v1 submitted 27 March, 2018; originally announced March 2018.

Comments: accepted for publication in M3AS

MSC Class: 35P15; 81Q10; 65N25; 65N55; 47B80

An Overview of the MHONGOOSE Survey: Observing Nearby Galaxies with MeerKAT

Authors: W. J. G. de Blok, E. A. K. Adams, P. Amram, E. Athanassoula, I. Bagetakos, C. Balkowski, M. A. Bershady, R. Beswick, F. Bigiel, S. -L. Blyth, A. Bosma, R. S. Booth, A. Bouchard, E. Brinks, C. Carignan, L. Chemin, F. Combes, J. Conway, E. C. Elson, J. English, B. Epinat, B. S. Frank, J. Fiege, F. Fraternali, J. S. Gallagher , et al. (43 additional authors not shown)

Abstract: MHONGOOSE is a deep survey of the neutral hydrogen distribution in a representative sample of 30 nearby disk and dwarf galaxies with HI masses from 10^6 to ~10^{11} M_sun, and luminosities from M_R ~ -12 to M_R ~ -22. The sample is selected to uniformly cover the available range in log(M_HI). Our extremely deep observations, down to HI column density limits of well below 10^{18} cm^{-2} - or a few… ▽ More MHONGOOSE is a deep survey of the neutral hydrogen distribution in a representative sample of 30 nearby disk and dwarf galaxies with HI masses from 10^6 to ~10^{11} M_sun, and luminosities from M_R ~ -12 to M_R ~ -22. The sample is selected to uniformly cover the available range in log(M_HI). Our extremely deep observations, down to HI column density limits of well below 10^{18} cm^{-2} - or a few hundred times fainter than the typical HI disks in galaxies - will directly detect the effects of cold accretion from the intergalactic medium and the links with the cosmic web. These observations will be the first ever to probe the very low-column density neutral gas in galaxies at these high resolutions. Combination with data at other wavelengths, most of it already available, will enable accurate modelling of the properties and evolution of the mass components in these galaxies and link these with the effects of environment, dark matter distribution, and other fundamental properties such as halo mass and angular momentum. MHONGOOSE can already start addressing some of the SKA-1 science goals and will provide a comprehensive inventory of the processes driving the transformation and evolution of galaxies in the nearby universe at high resolution and over 5 orders of magnitude in column density. It will be a Nearby Galaxies Legacy Survey that will be unsurpassed until the advent of the SKA, and can serve as a highly visible, lasting statement of MeerKAT's capabilities. △ Less

Submitted 25 September, 2017; originally announced September 2017.

Comments: To be published in "MeerKAT Science: On the Pathway to the SKA". Proceedings of Science. Workshop held 25-27 May, 2016 Stellenbosch, South Africa

Numerical homogenization of H(curl)-problems

Authors: Dietmar Gallistl, Patrick Henning, Barbara Verfürth

Abstract: If an elliptic differential operator associated with an $\mathbf{H}(\mathrm{curl})$-problem involves rough (rapidly varying) coefficients, then solutions to the corresponding $\mathbf{H}(\mathrm{curl})$-problem admit typically very low regularity, which leads to arbitrarily bad convergence rates for conventional numerical schemes. The goal of this paper is to show that the missing regularity can b… ▽ More If an elliptic differential operator associated with an $\mathbf{H}(\mathrm{curl})$-problem involves rough (rapidly varying) coefficients, then solutions to the corresponding $\mathbf{H}(\mathrm{curl})$-problem admit typically very low regularity, which leads to arbitrarily bad convergence rates for conventional numerical schemes. The goal of this paper is to show that the missing regularity can be compensated through a corrector operator. More precisely, we consider the lowest order Nédélec finite element space and show the existence of a linear corrector operator with four central properties: it is computable, $\mathbf{H}(\mathrm{curl})$-stable, quasi-local and allows for a correction of coarse finite element functions so that first-order estimates (in terms of the coarse mesh-size) in the $\mathbf{H}(\mathrm{curl})$ norm are obtained provided the right-hand side belongs to $\mathbf{H}(\mathrm{div})$. With these four properties, a practical application is to construct generalized finite element spaces which can be straightforwardly used in a Galerkin method. In particular, this characterizes a homogenized solution and a first order corrector, including corresponding quantitative error estimates without the requirement of scale separation. △ Less

Submitted 9 June, 2017; originally announced June 2017.

MSC Class: 35Q61; 65N12; 65N15; 65N30; 78M10

Crank-Nicolson Galerkin approximations to nonlinear Schrödinger equations with rough potentials

Authors: Patrick Henning, Daniel Peterseim

Abstract: This paper analyses the numerical solution of a class of non-linear Schrödinger equations by Galerkin finite elements in space and a mass- and energy conserving variant of the Crank-Nicolson method due to Sanz-Serna in time. The novel aspects of the analysis are the incorporation of rough, discontinuous potentials in the context of weak and strong disorder, the consideration of some general class… ▽ More This paper analyses the numerical solution of a class of non-linear Schrödinger equations by Galerkin finite elements in space and a mass- and energy conserving variant of the Crank-Nicolson method due to Sanz-Serna in time. The novel aspects of the analysis are the incorporation of rough, discontinuous potentials in the context of weak and strong disorder, the consideration of some general class of non-linearities, and the proof of convergence with rates in $L^{\infty}(L^2)$ under moderate regularity assumptions that are compatible with discontinuous potentials. For sufficiently smooth potentials, the rates are optimal without any coupling condition between the time step size and the spatial mesh width. △ Less

Submitted 17 June, 2017; v1 submitted 7 August, 2016; originally announced August 2016.

Highest Redshift Image of Neutral Hydrogen in Emission: A CHILES Detection of a Starbursting Galaxy at z=0.376

Authors: Ximena Fernández, Hansung B. Gim, J. H. van Gorkom, Min S. Yun, Emmanuel Momjian, Attila Popping, Laura Chomiuk, Kelley M. Hess, Lucas Hunt, Kathryn Kreckel, Danielle Lucero, Natasha Maddox, Tom Oosterloo, D. J. Pisano, M. A. W. Verheijen, Christopher A. Hales, Aeree Chung, Richard Dodson, Kumar Golap, Julia Gross, Patricia Henning, John Hibbard, Yara L. Jaffé, Jennifer Donovan Meyer, Martin Meyer , et al. (10 additional authors not shown)

Abstract: Our current understanding of galaxy evolution still has many uncertainties associated with the details of accretion, processing, and removal of gas across cosmic time. The next generation of radio telescopes will image the neutral hydrogen (HI) in galaxies over large volumes at high redshifts, which will provide key insights into these processes. We are conducting the COSMOS HI Large Extragalactic… ▽ More Our current understanding of galaxy evolution still has many uncertainties associated with the details of accretion, processing, and removal of gas across cosmic time. The next generation of radio telescopes will image the neutral hydrogen (HI) in galaxies over large volumes at high redshifts, which will provide key insights into these processes. We are conducting the COSMOS HI Large Extragalactic Survey (CHILES) with the Karl G. Jansky Very Large Array, which is the first survey to simultaneously observe HI from z=0 to z~0.5. Here, we report the highest redshift HI 21-cm detection in emission to date of the luminous infrared galaxy (LIRG) COSMOS J100054.83+023126.2 at z=0.376 with the first 178 hours of CHILES data. The total HI mass is $(2.9\pm1.0)\times10^{10}~M_\odot$, and the spatial distribution is asymmetric and extends beyond the galaxy. While optically the galaxy looks undisturbed, the HI distribution suggests an interaction with candidate a candidate companion. In addition, we present follow-up Large Millimeter Telescope CO observations that show it is rich in molecular hydrogen, with a range of possible masses of $(1.8-9.9)\times10^{10}~M_\odot$. This is the first study of the HI and CO in emission for a single galaxy beyond z~0.2. △ Less

Submitted 31 May, 2016; originally announced June 2016.

Comments: 7 pages, 5 figures, accepted for publication in ApJL

Multiscale methods for wave problems in heterogeneous media

Authors: Assyr Abdulle, Patrick Henning

Abstract: In this paper we give a survey on various multiscale methods for the numerical solution of second order hyperbolic equations in highly heterogeneous media. We concentrate on the wave equation and distinguish between two classes of applications. First we discuss numerical methods for the wave equation in heterogeneous media without scale separation. Such a setting is for instance encountered in the… ▽ More In this paper we give a survey on various multiscale methods for the numerical solution of second order hyperbolic equations in highly heterogeneous media. We concentrate on the wave equation and distinguish between two classes of applications. First we discuss numerical methods for the wave equation in heterogeneous media without scale separation. Such a setting is for instance encountered in the geosciences, where natural structures often exhibit a continuum of different scales, that all need to be resolved numerically to get meaningful approximations. Approaches tailored for these settings typically involve the construction of generalized finite element spaces, where the basis functions incorporate information about the data variations. In the second part of the paper, we discuss numerical methods for the case of structured media with scale separation. This setting is for instance encountered in engineering sciences, where materials are often artificially designed. If this is the case, the structure and the scale separation can be explicitly exploited to compute appropriate homogenized/upscaled wave models that only exhibit a single coarse scale and that can be hence solved at significantly reduced computational costs. △ Less

Submitted 25 May, 2016; originally announced May 2016.

The WSRT ZoA Perseus-Pisces Filament wide-field HI imaging survey I. HI catalogue and atlas

Authors: M. Ramatsoku, M. A. W Verheijen, R. C. Kraan-Korteweg, G. I. G. Józsa, A. C. Schröder, T. H. Jarrett, E. C Elson, W. van Driel, W. J. G. de Blok, P. A. Henning

Abstract: We present results of a blind 21cm HI-line imaging survey of a galaxy overdensity located behind the Milky Way at $\ell,b$ $\approx$ 160 deg, 0.5 deg. The overdensity corresponds to a Zone-of-Avoidance crossing of the Perseus-Pisces Supercluster filament. Although it is known that this filament contains an X-ray galaxy cluster (3C129) hosting two strong radio galaxies, little is known about galaxi… ▽ More We present results of a blind 21cm HI-line imaging survey of a galaxy overdensity located behind the Milky Way at $\ell,b$ $\approx$ 160 deg, 0.5 deg. The overdensity corresponds to a Zone-of-Avoidance crossing of the Perseus-Pisces Supercluster filament. Although it is known that this filament contains an X-ray galaxy cluster (3C129) hosting two strong radio galaxies, little is known about galaxies associated with this potentially rich cluster because of the high Galactic dust extinction. We mapped a sky area of $\sim$9.6 sq.deg using the Westerbork Synthesis Radio Telescope in a hexagonal mosaic of 35 pointings observed for 12 hours each, in the radial velocity range $cz = 2400 - 16600$ km/s. The survey has a sensitivity of 0.36 mJy/beam rms at a velocity resolution of 16.5 km/s. We detected 211 galaxies, 62% of which have a near-infrared counterpart in the UKIDSS Galactic Plane Survey. We present a catalogue of the HI properties and an HI atlas containing total intensity maps, position-velocity diagrams, global HI profiles and UKIDSS counterpart images. For the resolved galaxies we also present HI velocity fields and radial HI surface density profiles. A brief analysis of the structures outlined by these galaxies finds that 87 of them lie at the distance of the Perseus-Pisces Supercluster ($cz \sim 4000 - 8000$ km/s) and seem to form part of the 3C129 cluster. Further 72 detections trace an overdensity at a velocity of $cz \approx$ 10000 km/s and seem to coincide with a structure predicted from mass density reconstructions in the first 2MASS Redshift Survey. △ Less

Submitted 9 May, 2016; originally announced May 2016.

Comments: 28 pages, 18 figures, 6 Tables, accepted for publication in MNRAS

A multiscale method for linear elasticity reducing Poisson locking

Authors: Patrick Henning, Anna Persson

Abstract: We propose a generalized finite element method for linear elasticity equations with highly varying and oscillating coefficients. The method is formulated in the framework of localized orthogonal decomposition techniques introduced by Målqvist and Peterseim (Math. Comp., 83(290): 2583--2603, 2014). Assuming only $L_\infty$-coefficients we prove linear convergence in the $H^1$-norm, also for materia… ▽ More We propose a generalized finite element method for linear elasticity equations with highly varying and oscillating coefficients. The method is formulated in the framework of localized orthogonal decomposition techniques introduced by Målqvist and Peterseim (Math. Comp., 83(290): 2583--2603, 2014). Assuming only $L_\infty$-coefficients we prove linear convergence in the $H^1$-norm, also for materials with large Lamé parameter $λ$. The theoretical a priori error estimate is confirmed by numerical examples. △ Less

Submitted 31 March, 2016; originally announced March 2016.

The Parkes HI Zone of Avoidance Survey

Authors: L. Staveley-Smith, R. C. Kraan-Korteweg, A. C. Schröder, P. A. Henning, B. S. Koribalski, I. M. Stewart, G. Heald

Abstract: A blind HI survey of the extragalactic sky behind the southern Milky Way has been conducted with the multibeam receiver on the 64-m Parkes radio telescope. The survey covers the Galactic longitude range 212 < l < 36 and Galactic latitudes |b| < 5, and yields 883 galaxies to a recessional velocity of 12,000 km/s. The survey covers the sky within the HIPASS area to greater sensitivity, finding lower… ▽ More A blind HI survey of the extragalactic sky behind the southern Milky Way has been conducted with the multibeam receiver on the 64-m Parkes radio telescope. The survey covers the Galactic longitude range 212 < l < 36 and Galactic latitudes |b| < 5, and yields 883 galaxies to a recessional velocity of 12,000 km/s. The survey covers the sky within the HIPASS area to greater sensitivity, finding lower HI-mass galaxies at all distances, and probing more completely the large-scale structures at and beyond the distance of the Great Attractor. Fifty-one percent of the HI detections have an optical/NIR counterpart in the literature. A further 27% have new counterparts found in existing, or newly obtained, optical/NIR images. The counterpart rate drops in regions of high foreground stellar crowding and extinction, and for low-HI mass objects. Only 8% of all counterparts have a previous optical redshift measurement. A notable new galaxy is HIZOA J1353-58, a possible companion to the Circinus galaxy. Merging this catalog with the similarly-conducted northern extension (Donley et al. 2005), large-scale structures are delineated, including those within the Puppis and Great Attractor regions, and the Local Void. Several newly-identified structures are revealed here for the first time. Three new galaxy concentrations (NW1, NW2 and NW3) are key in confirming the diagonal crossing of the Great Attractor Wall between the Norma cluster and the CIZA J1324.7-5736 cluster. Further contributors to the general mass overdensity in that area are two new clusters (CW1 and CW2) in the nearer Centaurus Wall, one of which forms part of the striking 180 deg (100/h Mpc) long filament that dominates the southern sky at velocities of ~3000 km/s, and the suggestion of a further Wall at the Great Attractor distance at slightly higher longitudes. △ Less

Submitted 10 February, 2016; v1 submitted 9 February, 2016; originally announced February 2016.

Comments: Published in Astronomical Journal 9 February 2016 (accepted 26 September 2015); 42 pages, 7 tables, 18 figures, main figures data tables only available in the on-line version of journal

Journal ref: Astronomical Journal, 151, 52 (2016)

Efficient implementation of the Localized Orthogonal Decomposition method

Authors: Christian Engwer, Patrick Henning, Axel Målqvist, Daniel Peterseim

Abstract: In this paper we present algorithms for an efficient implementation of the Localized Orthogonal Decomposition method (LOD). The LOD is a multiscale method for the numerical simulation of partial differential equations with a continuum of inseparable scales. We show how the method can be implemented in a fairly standard Finite Element framework and discuss its realization for different types of pro… ▽ More In this paper we present algorithms for an efficient implementation of the Localized Orthogonal Decomposition method (LOD). The LOD is a multiscale method for the numerical simulation of partial differential equations with a continuum of inseparable scales. We show how the method can be implemented in a fairly standard Finite Element framework and discuss its realization for different types of problems, such as linear elliptic problems with rough coefficients and linear eigenvalue problems. △ Less

Submitted 20 February, 2019; v1 submitted 4 February, 2016; originally announced February 2016.

A new Heterogeneous Multiscale Method for time-harmonic Maxwell's equations based on divergence-regularization

Authors: Patrick Henning, Mario Ohlberger, Barbara Verfürth

Abstract: In this paper, we suggest a new heterogeneous multiscale method (HMM) for the time-harmonic Maxwell equations in locally periodic media. The method is constructed by using a divergence-regularization in one of the cell problems. This allows us to introduce fine-scale correctors that are not subject to a cumbersome divergence-free constraint and which can hence easily be implemented. To analyze the… ▽ More In this paper, we suggest a new heterogeneous multiscale method (HMM) for the time-harmonic Maxwell equations in locally periodic media. The method is constructed by using a divergence-regularization in one of the cell problems. This allows us to introduce fine-scale correctors that are not subject to a cumbersome divergence-free constraint and which can hence easily be implemented. To analyze the method, we first revisit classical homogenization theory for time-harmonic Maxwell equations and derive a new homogenization result that makes use of the divergence-regularization in the two-scale homogenized equation. We then show that the HMM is equivalent to a discretization of this equation. In particular, writing both problems in a fully coupled two-scale formulation is the crucial starting point for a corresponding numerical analysis of the method. With this approach we are able to prove rigorous a priori error estimates in the $\mathbf{H}(\mbox{curl})$- and the $H^{-1}$-norm and we derive reliable and efficient localized residual-based a posteriori error estimates. △ Less

Submitted 14 September, 2015; v1 submitted 10 September, 2015; originally announced September 2015.

Comments: 31 pages, fixed typos in this version

MSC Class: 35B27; 65N15; 65N30; 78M40

The ALFA ZOA Deep Survey: First Results

Authors: Travis P. McIntyre, Patricia A. Henning, Robert F. Minchin, Emmanuel Momjian, Zhon Butcher

Abstract: The Arecibo L-Band Feed Array Zone of Avoidance (ALFA ZOA) Deep Survey is the deepest and most sensitive blind Hi survey undertaken in the ZOA. ALFA ZOA Deep will cover about 300 square degrees of sky behind the Galactic plane in both the inner (30 deg < l < 75 deg; b < |2 deg|) and outer (175 deg < l < 207 deg; -2 deg < b < +1 deg) Galaxy, using the Arecibo Radio Telescope. First results from the… ▽ More The Arecibo L-Band Feed Array Zone of Avoidance (ALFA ZOA) Deep Survey is the deepest and most sensitive blind Hi survey undertaken in the ZOA. ALFA ZOA Deep will cover about 300 square degrees of sky behind the Galactic plane in both the inner (30 deg < l < 75 deg; b < |2 deg|) and outer (175 deg < l < 207 deg; -2 deg < b < +1 deg) Galaxy, using the Arecibo Radio Telescope. First results from the survey have found 61 galaxies within a 15 square degree area centered on l = 192 deg and b = -2 deg. The survey reached its expected sensitivity of rms = 1 mJy at 9 km/s channel resolution, and is shown to be complete above integrated flux, F_HI = 0.5 Jy km/s. The positional accuracy of the survey is 28 arcsec and detections are found out to a recessional velocity of nearly 19,000 km/s. The survey confirms the extent of the Orion and Abell 539 clusters behind the plane of the Milky Way and discovers expansive voids, at 10,000 km/s and 18,000 km/s. 26 detections (43%) have a counterpart in the literature, but only two of these have known redshift. Counterparts are 20% less common beyond v_hel = 10,000 km/s and 33% less common at extinctions higher than AB = 3.5 mag. ALFA ZOA Deep survey is able to probe large scale structure beyond redshifts that even the most modern wide-angle surveys have been able to detect in the Zone of Avoidance at any wavelength. △ Less

Submitted 23 June, 2015; originally announced June 2015.

Comments: 49 pages, Accepted by AJ, AASTeX5.2

Monitoring the Sky with the Prototype All-Sky Imager on the LWA1

Authors: K. S. Obenberger, G. B. Taylor, J. M. Hartman, T. E. Clarke, J. Dowell, A. Dubois, D. Dubois, P. A. Henning, J. Lazio, S. Michalak, F. K. Schinzel

Abstract: We present a description of the Prototype All-Sky Imager (PASI), a backend correlator and imager of the first station of the Long Wavelength Array (LWA1). PASI cross-correlates a live stream of 260 dual-polarization dipole antennas of the LWA1, creates all-sky images, and uploads them to the LWA-TV website in near real-time. PASI has recorded over 13,000 hours of all-sky images at frequencies betw… ▽ More We present a description of the Prototype All-Sky Imager (PASI), a backend correlator and imager of the first station of the Long Wavelength Array (LWA1). PASI cross-correlates a live stream of 260 dual-polarization dipole antennas of the LWA1, creates all-sky images, and uploads them to the LWA-TV website in near real-time. PASI has recorded over 13,000 hours of all-sky images at frequencies between 10 and 88 MHz creating opportunities for new research and discoveries. We also report rate density and pulse energy density limits on transients at 38, 52, and 74 MHz, for pulse widths of 5 s. We limit transients at those frequencies with pulse energy densities of $>2.7\times 10^{-23}$, $>1.1\times 10^{-23}$, and $>2.8\times 10^{-23}$ J m$^{-2}$ Hz$^{-1}$ to have rate densities $<1.2\times10^{-4}$, $<5.6\times10^{-4}$, and $<7.2\times10^{-4}$ yr$^{-1}$ deg$^{-2}$ △ Less

Submitted 17 March, 2015; originally announced March 2015.

Comments: 27 pages, 10 Figures, 1 Table

The Finite Element Method for the time-dependent Gross-Pitaevskii equation with angular momentum rotation

Authors: Patrick Henning, Axel Målqvist

Abstract: We consider the time-dependent Gross-Pitaevskii equation describing the dynamics of rotating Bose-Einstein condensates and its discretization with the finite element method. We analyze a mass conserving Crank-Nicolson-type discretization and prove corresponding a priori error estimates with respect to the maximum norm in time and the $L^2$- and energy-norm in space. The estimates show that we obta… ▽ More We consider the time-dependent Gross-Pitaevskii equation describing the dynamics of rotating Bose-Einstein condensates and its discretization with the finite element method. We analyze a mass conserving Crank-Nicolson-type discretization and prove corresponding a priori error estimates with respect to the maximum norm in time and the $L^2$- and energy-norm in space. The estimates show that we obtain optimal convergence rates under the assumption of additional regularity for the solution to the Gross-Pitaevskii equation. We demonstrate the performance of the method in numerical experiments. △ Less

Submitted 7 June, 2016; v1 submitted 17 February, 2015; originally announced February 2015.

Multiscale mixed finite elements

Authors: Fredrik Hellman, Patrick Henning, Axel Målqvist

Abstract: In this work, we propose a mixed finite element method for solving elliptic multiscale problems based on a localized orthogonal decomposition (LOD) of Raviart-Thomas finite element spaces. It requires to solve local problems in small patches around the elements of a coarse grid. These computations can be perfectly parallelized and are cheap to perform. Using the results of these patch problems, we… ▽ More In this work, we propose a mixed finite element method for solving elliptic multiscale problems based on a localized orthogonal decomposition (LOD) of Raviart-Thomas finite element spaces. It requires to solve local problems in small patches around the elements of a coarse grid. These computations can be perfectly parallelized and are cheap to perform. Using the results of these patch problems, we construct a low dimensional multiscale mixed finite element space with very high approximation properties. This space can be used for solving the original saddle point problem in an efficient way. We prove convergence of our approach, independent of structural assumptions or scale separation. Finally, we demonstrate the applicability of our method by presenting a variety of numerical experiments, including a comparison with an MsFEM approach. △ Less

Submitted 20 June, 2016; v1 submitted 22 January, 2015; originally announced January 2015.

MSC Class: 35J15; 35M10; 65N12; 65N30; 76S05

A reduced basis localized orthogonal decomposition

Authors: Assyr Abdulle, Patrick Henning

Abstract: In this work we combine the framework of the Reduced Basis method (RB) with the framework of the Localized Orthogonal Decomposition (LOD) in order to solve parametrized elliptic multiscale problems. The idea of the LOD is to split a high dimensional Finite Element space into a low dimensional space with comparably good approximation properties and a remainder space with negligible information. The… ▽ More In this work we combine the framework of the Reduced Basis method (RB) with the framework of the Localized Orthogonal Decomposition (LOD) in order to solve parametrized elliptic multiscale problems. The idea of the LOD is to split a high dimensional Finite Element space into a low dimensional space with comparably good approximation properties and a remainder space with negligible information. The low dimensional space is spanned by locally supported basis functions associated with the node of a coarse mesh obtained by solving decoupled local problems. However, for parameter dependent multiscale problems, the local basis has to be computed repeatedly for each choice of the parameter. To overcome this issue, we propose an RB approach to compute in an "offline" stage LOD for suitable representative parameters. The online solution of the multiscale problems can then be obtained in a coarse space (thanks to the LOD decomposition) and for an arbitrary value of the parameters (thanks to a suitable "interpolation" of the selected RB). The online RB-LOD has a basis with local support and leads to sparse systems. Applications of the strategy to both linear and nonlinear problems are given. △ Less

Submitted 12 May, 2015; v1 submitted 13 October, 2014; originally announced October 2014.

MSC Class: 65N30; 65M60; 74Q05; 74Q15

Journal ref: Journal of Computational Physics, Vol. 295, pp 379--401, 2015

Localized orthogonal decomposition method for the wave equation with a continuum of scales

Authors: Assyr Abdulle, Patrick Henning

Abstract: This paper is devoted to numerical approximations for the wave equation with a multiscale character. Our approach is formulated in the framework of the Localized Orthogonal Decomposition (LOD) interpreted as a numerical homogenization with an $L^2$-projection. We derive explicit convergence rates of the method in the $L^{\infty}(L^2)$-, $W^{1,\infty}(L^2)$- and $L^{\infty}(H^1)$-norms without any… ▽ More This paper is devoted to numerical approximations for the wave equation with a multiscale character. Our approach is formulated in the framework of the Localized Orthogonal Decomposition (LOD) interpreted as a numerical homogenization with an $L^2$-projection. We derive explicit convergence rates of the method in the $L^{\infty}(L^2)$-, $W^{1,\infty}(L^2)$- and $L^{\infty}(H^1)$-norms without any assumptions on higher order space regularity or scale-separation. The order of the convergence rates depends on further graded assumptions on the initial data. We also prove the convergence of the method in the framework of G-convergence without any structural assumptions on the initial data, i.e. without assuming that it is well-prepared. This rigorously justifies the method. Finally, the performance of the method is demonstrated in numerical experiments. △ Less

Submitted 22 September, 2015; v1 submitted 24 June, 2014; originally announced June 2014.

Detection of Radio Emission from Fireballs

Authors: K. S. Obenberger, G. B. Taylor, J. M. Hartman, J. Dowell, S. W. Ellingson, J. F. Helmboldt, P. A. Henning, M. Kavic, F. K. Schinzel, J. H. Simonetti, K. Stovall, T. L. Wilson

Abstract: We present the findings from the Prototype All-Sky Imager (PASI), a backend correlator of the first station of the Long Wavelength Array (LWA1), which has recorded over 11,000 hours of all-sky images at frequencies between 25 and 75 MHz. In a search of this data for radio transients, we have found 49 long (10s of seconds) duration transients. Ten of these transients correlate both spatially and te… ▽ More We present the findings from the Prototype All-Sky Imager (PASI), a backend correlator of the first station of the Long Wavelength Array (LWA1), which has recorded over 11,000 hours of all-sky images at frequencies between 25 and 75 MHz. In a search of this data for radio transients, we have found 49 long (10s of seconds) duration transients. Ten of these transients correlate both spatially and temporally with large meteors (fireballs), and their signatures suggest that fireballs emit a previously undiscovered low frequency, non-thermal pulse. This emission provides a new probe into the physics of meteors and identifies a new form of naturally occurring radio transient foreground. △ Less

Submitted 26 May, 2014; originally announced May 2014.

Comments: 14 Pages, 5 Figures

Journal ref: ApJ, 788, L26 2014

On Multiscale Methods in Petrov-Galerkin formulation

Authors: Daniel Elfverson, Victor Ginting, Patrick Henning

Abstract: In this work we investigate the advantages of multiscale methods in Petrov-Galerkin (PG) formulation in a general framework. The framework is based on a localized orthogonal decomposition of a high dimensional solution space into a low dimensional multiscale space with good approximation properties and a high dimensional remainder space{, which only contains negligible fine scale information}. The… ▽ More In this work we investigate the advantages of multiscale methods in Petrov-Galerkin (PG) formulation in a general framework. The framework is based on a localized orthogonal decomposition of a high dimensional solution space into a low dimensional multiscale space with good approximation properties and a high dimensional remainder space{, which only contains negligible fine scale information}. The multiscale space can then be used to obtain accurate Galerkin approximations. As a model problem we consider the Poisson equation. We prove that a Petrov-Galerkin formulation does not suffer from a significant loss of accuracy, and still preserve the convergence order of the original multiscale method. We also prove inf-sup stability of a PG Continuous and a Discontinuous Galerkin Finite Element multiscale method. Furthermore, we demonstrate that the Petrov-Galerkin method can decrease the computational complexity significantly, allowing for more efficient solution algorithms. As another application of the framework, we show how the Petrov-Galerkin framework can be used to construct a locally mass conservative solver for two-phase flow simulation that employs the Buckley-Leverett equation. To achieve this, we couple a PG Discontinuous Galerkin Finite Element method with an upwind scheme for a hyperbolic conservation law. △ Less

Submitted 30 December, 2014; v1 submitted 22 May, 2014; originally announced May 2014.

Limits on GRB Prompt Radio Emission Using the LWA1

Authors: K. S. Obenberger, J. M. Hartman, G. B. Taylor, J. Craig, J. Dowell, J. F. Helmboldt, P. A. Henning, F. K. Schinzel, T. L. Wilson

Abstract: As a backend to the first station of the Long Wavelength Array (LWA1) the Prototype All Sky Imager (PASI) has been imaging the sky $>$ -26$^{\circ}$ declination during 34 Gamma Ray Bursts (GRBs) between January 2012 and May 2013. Using this data we were able to put the most stringent limits to date on prompt low frequency emission from GRBs. While our limits depend on the zenith angle of the obser… ▽ More As a backend to the first station of the Long Wavelength Array (LWA1) the Prototype All Sky Imager (PASI) has been imaging the sky $>$ -26$^{\circ}$ declination during 34 Gamma Ray Bursts (GRBs) between January 2012 and May 2013. Using this data we were able to put the most stringent limits to date on prompt low frequency emission from GRBs. While our limits depend on the zenith angle of the observed GRB, we estimate a 1$σ$ RMS sensitivity of 68, 65 and 70 Jy for 5 second integrations at 37.9, 52.0, and 74.0 MHz at zenith. These limits are relevant for pulses $\geq$ 5 s and are limited by dispersion smearing. For pulses of length 5 s we are limited to dispersion measures ($DM$s) $\leq$ 220, 570, and 1,600 pc cm$^{-3}$ for the frequencies above. For pulses lasting longer than 5s, the $DM$ limits increase linearly with the duration of the pulse. We also report two interesting transients, which are, as of yet, of unknown origin, and are not coincident with any known GRBs. For general transients, we give rate density limits of $\leq$ $7.5\times10^{-3}$, $2.9\times10^{-2}$, and $1.4\times10^{-2}$ yr$^{-1}$ deg$^{-2}$ with pulse energy densities $>1.3\times 10^{-22}$, $1.1\times 10^{-22}$, and $1.4\times 10^{-22}$ J m$^{-2}$ Hz$^{-1}$ and pulse widths of 5 s at the frequencies given above. △ Less

Submitted 14 March, 2014; originally announced March 2014.

Comments: 11 pages, 6 figure, 3 tables, accepted for publication in ApJ

Multiscale Partition of Unity

Authors: Daniel Peterseim, Patrick Henning, Philipp Morgenstern

Abstract: We introduce a new Partition of Unity Method for the numerical homogenization of elliptic partial differential equations with arbitrarily rough coefficients. We do not restrict to a particular ansatz space or the existence of a finite element mesh. The method modifies a given partition of unity such that optimal convergence is achieved independent of oscillation or discontinuities of the diffusion… ▽ More We introduce a new Partition of Unity Method for the numerical homogenization of elliptic partial differential equations with arbitrarily rough coefficients. We do not restrict to a particular ansatz space or the existence of a finite element mesh. The method modifies a given partition of unity such that optimal convergence is achieved independent of oscillation or discontinuities of the diffusion coefficient. The modification is based on an orthogonal decomposition of the solution space while preserving the partition of unity property. This precomputation involves the solution of independent problems on local subdomains of selectable size. We deduce quantitative error estimates for the method that account for the chosen amount of localization. Numerical experiments illustrate the high approximation properties even for 'cheap' parameter choices. △ Less

Submitted 10 October, 2014; v1 submitted 20 December, 2013; originally announced December 2013.

Comments: Proceedings for Seventh International Workshop on Meshfree Methods for Partial Differential Equations, 18 pages, 3 figures

Report number: 1315 MSC Class: 65N99; 65B99; 74Q05 ACM Class: G.1.8

Detection and Flux Density Measurements of the Millisecond Pulsar J2145-0750 below 100 MHz

Authors: J. Dowell, P. S. Ray, G. B. Taylor, J. N. Blythe, T. Clarke, J. Craig, S. W. Ellingson, J. F. Helmboldt, P. A. Henning, T. J. W. Lazio, F. Schinzel, K. Stovall, C. N. Wolfe

Abstract: We present flux density measurements and pulse profiles for the millisecond pulsar PSR J2145-0750 spanning 37 to 81 MHz using data obtained from the first station of the Long Wavelength Array. These measurements represent the lowest frequency detection of pulsed emission from a millisecond pulsar to date. We find that the pulse profile is similar to that observed at 102 MHz. We also find that the… ▽ More We present flux density measurements and pulse profiles for the millisecond pulsar PSR J2145-0750 spanning 37 to 81 MHz using data obtained from the first station of the Long Wavelength Array. These measurements represent the lowest frequency detection of pulsed emission from a millisecond pulsar to date. We find that the pulse profile is similar to that observed at 102 MHz. We also find that the flux density spectrum between ~40 MHz to 5 GHz is suggestive of a break and may be better fit by a model that includes spectral curvature with a rollover around 730 MHz rather than a single power law. △ Less

Submitted 11 September, 2013; v1 submitted 2 September, 2013; originally announced September 2013.

Comments: accepted to the Astrophysical Journal Letters; 14 pages, 4 figures

Journal ref: ApJL, 775, L28 (2013)

Localized orthogonal decomposition techniques for boundary value problems

Authors: Patrick Henning, Axel Målqvist

Abstract: In this paper we propose a Local Orthogonal Decomposition method (LOD) for elliptic partial differential equations with inhomogeneous Dirichlet- and Neumann boundary conditions. For this purpose, we present new boundary correctors which preserve the common convergence rates of the LOD, even if the boundary condition has a rapidly oscillating fine scale structure. We prove a corresponding a-priori… ▽ More In this paper we propose a Local Orthogonal Decomposition method (LOD) for elliptic partial differential equations with inhomogeneous Dirichlet- and Neumann boundary conditions. For this purpose, we present new boundary correctors which preserve the common convergence rates of the LOD, even if the boundary condition has a rapidly oscillating fine scale structure. We prove a corresponding a-priori error estimate and present numerical experiments. We also demonstrate numerically that the method is reliable with respect to thin conductivity channels in the diffusion matrix. Accurate results are obtained without resolving these channels by the coarse grid and without using patches that contain the channels. △ Less

Submitted 17 July, 2014; v1 submitted 15 August, 2013; originally announced August 2013.

MSC Class: 35J15; 65N12; 65N30