BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Matthew Anderson

Matthew Anderson appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
gr-qcastro-ph.HEComputational ComplexityDistributed, Parallel, and Cluster Computingastro-ph.IMLogic in Computer ScienceArtificial Intelligenceastro-ph.COastro-ph.GAData Structures and AlgorithmsDiscrete Mathematicseess.SPeess.SYFormal Languages and Automata TheoryInformation TheoryMachine Learningmath.ITmath.NAMultiagent SystemsNumerical Analysisphysics.comp-phRoboticsSymbolic ComputationSystems and Control

Catalog footprint

What is connected

23works
24topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Published work

23 published item(s)

preprint2022arXiv

Design of the new CHARA instrument SILMARIL: pushing the sensitivity of a 3-beam combiner in the H- and K-bands

Optical interferometry is a powerful technique to achieve high angular resolution. However, its main issue is its lack of sensitivity, compared to other observation techniques. Efforts have been made in the previous decade to improve the sensitivity of optical interferometry, with instruments such as PIONIER and GRAVITY at VLTI, or MIRC-X and MYSTIC at CHARA. While those instruments pushed on sensitivity, their design focus was not the sensitivity but relative astrometric accuracy, imaging capability, or spectral resolution. Our goal is to build an instrument specifically designed to optimize for sensitivity. This meant focusing our design efforts on different parts of the instrument and investigating new technologies and techniques. First, we make use of the low-noise C-RED One camera using e-APD technology and provided by First Light Imaging, already used in the improvement of sensitivity in recent new instruments. We forego the use of single-mode fibers but still favor an image plane design that offers more sensitivity than a pupil plane layout. We also use a minimum number of optical elements to maximize the throughput of the design, using a long focal length cylindrical mirror. We chose to limit our design to 3 beams, to have the capability to obtain closure phases, but not dilute the incoming flux in more beam combinations. We also use in our design an edge filter to have the capability to observe H- and K-band at the same time. We use a low spectral resolution, allowing for group delay fringe tracking but maximizing the SNR of the fringes for each spectral channel. All these elements will lead to a typical limiting magnitude between 10 and 11 in both H- and K-bands.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Matrix Multiplication: Verifying Strong Uniquely Solvable Puzzles

Cohn and Umans proposed a framework for developing fast matrix multiplication algorithms based on the embedding computation in certain groups algebras. In subsequent work with Kleinberg and Szegedy, they connected this to the search for combinatorial objects called strong uniquely solvable puzzles (strong USPs). We begin a systematic computer-aided search for these objects. We develop and implement constraint-based algorithms build on reductions to $\mathrm{SAT}$ and $\mathrm{IP}$ to verify that puzzles are strong USPs, and to search for large strong USPs. We produce tight bounds on the maximum size of a strong USP for width $k \le 5$, construct puzzles of small width that are larger than previous work, and improve the upper bounds on strong USP size for $k \le 12$. Although our work only deals with puzzles of small-constant width, the strong USPs we find imply matrix multiplication algorithms that run in $O(n^ω)$ time with exponent $ω\le 2.66$. While our algorithms do not beat the fastest algorithms, our work provides evidence and, perhaps, a path to finding families of strong USPs that imply matrix multiplication algorithms that are more efficient than those currently known.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Automated Trajectory Synthesis for UAV Swarms Based on Resilient Data Collection Objectives

The use of Unmanned Aerial Vehicles (UAVs) for collecting data from remotely located sensor systems is emerging. The data can be time-sensitive and require to be transmitted to a data processing center. However, planning the trajectory of a collaborative UAV swarm depends on multi-fold constraints, such as data collection requirements, UAV maneuvering capacity, and budget limitation. Since a UAV may fail or be compromised, it is important to provide necessary resilience to such contingencies, thus ensuring data security. It is important to provide the UAVs with efficient spatio-temporal trajectories so that they can efficiently cover necessary data sources. In this work, we present Synth4UAV, a formal approach for automated synthesis of efficient trajectories for a UAV swarm by logically modeling the aerial space and data point topology, UAV moves, and associated constraints in terms of the turning and climbing angle, fuel usage, data collection point coverage, data freshness, and resiliency properties. We use efficient, logical formulas to encode and solve the complex model. The solution to the model provides the routing and maneuvering plan for each UAV, including the time to visit the points on the paths and corresponding fuel usage such that the necessary data points are visited while satisfying the resiliency requirements. We evaluate the proposed trajectory synthesizer, and the results show that the relationship among different parameters follow the requirements while the tool scales well with the problem size.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Design and Autonomous Stabilization of a Ballistically Launched Multirotor

Aircraft that can launch ballistically and convert to autonomous, free flying drones have applications in many areas such as emergency response, defense, and space exploration, where they can gather critical situational data using onboard sensors. This paper presents a ballistically launched, autonomously stabilizing multirotor prototype (SQUID, Streamlined Quick Unfolding Investigation Drone) with an onboard sensor suite, autonomy pipeline, and passive aerodynamic stability. We demonstrate autonomous transition from passive to vision based, active stabilization, confirming the ability of the multirotor to autonomously stabilize after a ballistic launch in a GPS denied environment.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Exploring Critical Collapse in the Semilinear Wave Equation using Space-Time Finite Elements

A fully implicit numerical approach based on the space-time finite element method is implemented for the semilinear wave equation in 1(space) + 1(time) and 2 + 1 dimensions to explore critical collapse and search for self-similar solutions. Previous work studied this behavior by exploring the threshold of singularity formation using time marching finite difference techniques while this work introduces an adaptive time parallel numerical method to the problem. The semilinear wave equation with a $p = 7$ term is examined in spherical symmetry. The impact of mesh refinement and the time additive Schwarz preconditioner in conjunction with Krylov Subspace Methods are examined.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Unequal mass binary neutron star mergers and multimessenger signals

We study the merger of binary neutron stars with different mass ratios adopting three different realistic, microphysical nuclear equations of state, as well as incorporating neutrino cooling effects. In particular, we concentrate on the influence of the equation of state on the gravitational wave signature and also on its role, in combination with neutrino cooling, in determining the properties of the resulting hypermassive neutron star, of the neutrinos produced, and of the ejected material. The ejecta we find are consistent with other recent studies that find that small mass ratios produce more ejecta than equal mass cases (up to some limit) and this ejecta is more neutron rich. This trend indicates the importance with future kilonovae observations of measuring the individual masses of an associated binary neutron star system, presumably from concurrent gravitational wave observations, in order to be able to extract information about the nuclear equation of state

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Identity Testing and Lower Bounds for Read-$k$ Oblivious Algebraic Branching Programs

Read-$k$ oblivious algebraic branching programs are a natural generalization of the well-studied model of read-once oblivious algebraic branching program (ROABPs). In this work, we give an exponential lower bound of $\exp(n/k^{O(k)})$ on the width of any read-$k$ oblivious ABP computing some explicit multilinear polynomial $f$ that is computed by a polynomial size depth-$3$ circuit. We also study the polynomial identity testing (PIT) problem for this model and obtain a white-box subexponential-time PIT algorithm. The algorithm runs in time $2^{\tilde{O}(n^{1-1/2^{k-1}})}$ and needs white box access only to know the order in which the variables appear in the ABP.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Relativistic Hydrodynamics with Wavelets

Methods to solve the relativistic hydrodynamic equations are a key computational kernel in a large number of astrophysics simulations and are crucial to understanding the electromagnetic signals that originate from the merger of astrophysical compact objects. Because of the many physical length scales present when simulating such mergers, these methods must be highly adaptive and capable of automatically resolving numerous localized features and instabilities that emerge throughout the computational domain across many temporal scales. While this has been historically accomplished with adaptive mesh refinement (AMR) based methods, alternatives based on wavelet bases and the wavelet transformation have recently achieved significant success in adaptive representation for advanced engineering applications. This work presents a new method for the integration of the relativistic hydrodynamic equations using iterated interpolating wavelets and introduces a highly adaptive implementation for multidimensional simulation. The wavelet coefficients provide a direct measure of the local approximation error for the solution and place collocation points that naturally adapt to the fluid flow while providing good conservation of fluid quantities. The resulting implementation, OAHU, is applied to a series of demanding one- and two-dimensional problems which explore high Lorentz factor outflows and the formation of several instabilities, including the Kelvin-Helmholtz instability and the Rayleigh-Taylor instability.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Magnetized Neutron Stars With Realistic Equations of State and Neutrino Cooling

We incorporate realistic, tabulated equations of state into fully relativistic simulations of magnetized neutron stars along with a neutrino leakage scheme which accounts for cooling via neutrino emission. Both these improvements utilize open-source code (GR1D) and tables from http://stellarcollapse.org. Our implementation makes use of a novel method for the calculation of the optical depth which simplifies its use with distributed adaptive mesh refinement. We present various tests with and without magnetization and preliminary results both from single stars and from the merger of a binary system.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

On Symmetric Circuits and Fixed-Point Logics

We study properties of relational structures such as graphs that are decided by families of Boolean circuits. Circuits that decide such properties are necessarily invariant to permutations of the elements of the input structures. We focus on families of circuits that are symmetric, i.e., circuits whose invariance is witnessed by automorphisms of the circuit induced by the permutation of the input structure. We show that the expressive power of such families is closely tied to definability in logic. In particular, we show that the queries defined on structures by uniform families of symmetric Boolean circuits with majority gates are exactly those definable in fixed-point logic with counting. This shows that inexpressibility results in the latter logic lead to lower bounds against polynomial-size families of symmetric circuits.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The Mass of the Central Black Hole in the Nearby Seyfert Galaxy NGC5273

We present the results of a reverberation-mapping program targeting NGC5273, a nearby early-type galaxy with a broad-lined active galactic nucleus. Over the course of the monitoring program, NGC5273 showed strong variability that allowed us to measure time delays in the responses of the broad optical recombination lines to changes in the continuum flux. A weighted average of these measurements results in a black hole mass determination of $M_{\rm BH} = (4.7 \pm 1.6) \times 10^6$ M$_{\odot}$. An estimate of the size of the black hole sphere of influence in NGC5273 puts it just at the limit of the resolution achievable with current ground-based large aperture telescopes. NGC5273 is therefore an important future target for a black hole mass determination from stellar dynamical modeling, especially because it is the only nearby early-type galaxy hosting an AGN with a reverberation-based mass, allowing the best comparison for the masses determined from these two techniques.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Gravitational and electromagnetic outputs from binary neutron star mergers

The late stage of an inspiraling neutron star binary gives rise to strong gravitational wave emission due to its highly dynamic, strong gravity. Moreover, interactions between the stellar magnetospheres can produce considerable electromagnetic radiation. We study this scenario using fully general relativistic, resistive magneto-hydrodynamics simulations. We show that these interactions extract kinetic energy from the system, dissipate heat, and power radiative Poynting flux, as well as develop current sheets. Our results indicate that this power can: (i) outshine pulsars in binaries, (ii) display a distinctive angular- and time-dependent pattern, and (iii) radiate within large opening angles. These properties suggest that some binary neutron star mergers are ideal candidates for multimessenger astronomy.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Independent Vector Analysis: Identification Conditions and Performance Bounds

Recently, an extension of independent component analysis (ICA) from one to multiple datasets, termed independent vector analysis (IVA), has been the subject of significant research interest. IVA has also been shown to be a generalization of Hotelling's canonical correlation analysis. In this paper, we provide the identification conditions for a general IVA formulation, which accounts for linear, nonlinear, and sample-to-sample dependencies. The identification conditions are a generalization of previous results for ICA and for IVA when samples are independently and identically distributed. Furthermore, a principal aim of IVA is the identification of dependent sources between datasets. Thus, we provide the additional conditions for when the arbitrary ordering of the sources within each dataset is common. Performance bounds in terms of the Cramer-Rao lower bound are also provided for the demixing matrices and interference to source ratio. The performance of two IVA algorithms are compared to the theoretical bounds.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Linking electromagnetic and gravitational radiation in coalescing binary neutron stars

We expand on our study of the gravitational and electromagnetic emissions from the late stage of an inspiraling neutron star binary as presented in Ref. \cite{Palenzuela:2013hu}. Interactions between the stellar magnetospheres, driven by the extreme dynamics of the merger, can yield considerable outflows. We study the gravitational and electromagnetic waves produced during the inspiral and merger of a binary neutron star system using a full relativistic, resistive MHD evolution code. We show that the interaction between the stellar magnetospheres extracts kinetic energy from the system and powers radiative Poynting flux and heat dissipation. These features depend strongly on the configuration of the initial stellar magnetic moments. Our results indicate that this power can strongly outshine pulsars in binaries and have a distinctive angular and time-dependent pattern. Our discussion provides more detail than Ref. \cite{Palenzuela:2013hu}, showing clear evidence of the different effects taking place during the inspiral. Our simulations include a few milliseconds after the actual merger and study the dynamics of the magnetic fields during the formation of the hypermassive neutron star. We also briefly discuss the possibility of observing such emissions.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Maximum Matching and Linear Programming in Fixed-Point Logic with Counting

We establish the expressibility in fixed-point logic with counting (FPC) of a number of natural polynomial-time problems. In particular, we show that the size of a maximum matching in a graph is definable in FPC. This settles an open problem first posed by Blass, Gurevich and Shelah, who asked whether the existence of perfect matchings in general graphs could be determined in the more powerful formalism of choiceless polynomial time with counting. Our result is established by showing that the ellipsoid method for solving linear programs can be implemented in FPC. This allows us to prove that linear programs can be optimised in FPC if the corresponding separation oracle problem can be defined in FPC. On the way to defining a suitable separation oracle for the maximum matching problem, we provide FPC formulas defining maximum flows and canonical minimum cuts in capacitated graphs.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Neutron Star Evolutions using Tabulated Equations of State with a New Execution Model

The addition of nuclear and neutrino physics to general relativistic fluid codes allows for a more realistic description of hot nuclear matter in neutron star and black hole systems. This additional microphysics requires that each processor have access to large tables of data, such as equations of state, and in large simulations the memory required to store these tables locally can become excessive unless an alternative execution model is used. In this work we present relativistic fluid evolutions of a neutron star obtained using a message driven multi-threaded execution model known as ParalleX. These neutron star simulations would require substantial memory overhead dedicated entirely to the equation of state table if using a more traditional execution model. We introduce a ParalleX component based on Futures for accessing large tables of data, including out-of-core sized tables, which does not require substantial memory overhead and effectively hides any increased network latency.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Adaptive Mesh Refinement for Astrophysics Applications with ParalleX

Several applications in astrophysics require adequately resolving many physical and temporal scales which vary over several orders of magnitude. Adaptive mesh refinement techniques address this problem effectively but often result in constrained strong scaling performance. The ParalleX execution model is an experimental execution model that aims to expose new forms of program parallelism and eliminate any global barriers present in a scaling-impaired application such as adaptive mesh refinement. We present two astrophysics applications using the ParalleX execution model: a tabulated equation of state component for neutron star evolutions and a cosmology model evolution. Performance and strong scaling results from both simulations are presented. The tabulated equation of state data are distributed with transparent access over the nodes of the cluster. This allows seamless overlapping of computation with the latencies introduced by the remote access to the table. Because of the expected size increases to the equation of state table, this type of table partitioning for neutron star simulations is essential while the implementation is greatly simplified by ParalleX semantics.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

An Application Driven Analysis of the ParalleX Execution Model

Exascale systems, expected to emerge by the end of the next decade, will require the exploitation of billion-way parallelism at multiple hierarchical levels in order to achieve the desired sustained performance. The task of assessing future machine performance is approached by identifying the factors which currently challenge the scalability of parallel applications. It is suggested that the root cause of these challenges is the incoherent coupling between the current enabling technologies, such as Non-Uniform Memory Access of present multicore nodes equipped with optional hardware accelerators and the decades older execution model, i.e., the Communicating Sequential Processes (CSP) model best exemplified by the message passing interface (MPI) application programming interface. A new execution model, ParalleX, is introduced as an alternative to the CSP model. In this paper, an overview of the ParalleX execution model is presented along with details about a ParalleX-compliant runtime system implementation called High Performance ParalleX (HPX). Scaling and performance results for an adaptive mesh refinement numerical relativity application developed using HPX are discussed. The performance results of this HPX-based application are compared with a counterpart MPI-based mesh refinement code. The overheads associated with HPX are explored and hardware solutions are introduced for accelerating the runtime system.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Improving the scalability of parallel N-body applications with an event driven constraint based execution model

The scalability and efficiency of graph applications are significantly constrained by conventional systems and their supporting programming models. Technology trends like multicore, manycore, and heterogeneous system architectures are introducing further challenges and possibilities for emerging application domains such as graph applications. This paper explores the space of effective parallel execution of ephemeral graphs that are dynamically generated using the Barnes-Hut algorithm to exemplify dynamic workloads. The workloads are expressed using the semantics of an Exascale computing execution model called ParalleX. For comparison, results using conventional execution model semantics are also presented. We find improved load balancing during runtime and automatic parallelism discovery improving efficiency using the advanced semantics for Exascale computing.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Mergers of Magnetized Neutron Stars with Spinning Black Holes: Disruption, Accretion and Fallback

We investigate the merger of a neutron star (of compaction ratio $0.1$) in orbit about a spinning black hole in full general relativity with a mass ratio of $5:1$, allowing for the star to have an initial magnetization of $10^{12} {\rm Gauss}$. We present the resulting gravitational waveform and analyze the fallback accretion as the star is disrupted. The evolutions suggest no significant effects from the initial magnetization. We find that only a negligible amount of matter becomes unbound; $99%$ of the neutron star material has a fallback time of 10 seconds or shorter to reach the region of the central engine and that $99.99%$ of the star will interact with the central disk and black hole within 3 hours.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Post-merger electromagnetic emissions from disks perturbed by binary black holes

We simulate the possible emission from a disk perturbed by a recoiling super-massive black hole. To this end, we study radiation transfer from the system incorporating bremsstrahlung emission from a Maxwellian plasma and absorption given by Kramer's opacity law modified to incorporate blackbody effects. We employ this model in the radiation transfer integration to compute the luminosity at several frequencies, and compare with previous bremsstrahlung luminosity estimations from a transparent limit (in which the emissivity is integrated over the computational domain and over all frequencies) and with a simple thermal emission model. We find close agreement between the radiation transfer results and the estimated bremsstrahlung luminosity from previous work for electromagnetic signals above $10^{14}$ Hz. For lower frequencies, we find a self-eclipsing behavior in the disk, resulting in a strong intensity variability connected to the orbital period of the disk.

0 citations0 reviewsNo code linkedOpen access
preprint2004arXiv

Extended Lifetime in Computational Evolution of Isolated Black Holes

Solving the 4-d Einstein equations as evolution in time requires solving equations of two types: the four elliptic initial data (constraint) equations, followed by the six second order evolution equations. Analytically the constraint equations remain solved under the action of the evolution, and one approach is to simply monitor them ({\it unconstrained} evolution). The problem of the 3-d computational simulation of even a single isolated vacuum black hole has proven to be remarkably difficult. Recently, we have become aware of two publications that describe very long term evolution, at least for single isolated black holes. An essential feature in each of these results is {\it constraint subtraction}. Additionally, each of these approaches is based on what we call "modern," hyperbolic formulations of the Einstein equations. It is generally assumed, based on computational experience, that the use of such modern formulations is essential for long-term black hole stability. We report here on comparable lifetime results based on the much simpler ("traditional") $\dot g$ - $\dot K$ formulation. We have also carried out a series of {\it constrained} 3-d evolutions of single isolated black holes. We find that constraint solution can produce substantially stabilized long-term single hole evolutions. However, we have found that for large domains, neither constraint-subtracted nor constrained $\dot g$ - $\dot K$ evolutions carried out in Cartesian coordinates admit arbitrarily long-lived simulations. The failure appears to arise from features at the inner excision boundary; the behavior does generally improve with resolution.

0 citations0 reviewsNo code linkedOpen access